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Abstract 



The successes and fine-tuning problems of the Hot Big Bang theory of the Universe are 
briefly reviewed. Cosmological inflation alleviates those problems substantially and give 
rise to the primordial curvature perturbation with the properties observed in the Cosmic 
Microwave Background. It is shown how application of the quantum field theory in the 
exponentially expanding Universe leads to the conversion of quantum fluctuations into 
the classical field perturbation. The 5N formalism is reviewed and applied to calculate 
the primordial curvature perturbation £ for three examples: single field inflation, the 
end-of-inflation and the curvaton scenarios. 

The 5N formalism is extended to include the perturbation of the vector field. The 
latter is quantized in de Sitter space-time and it is found that in general the particle 
production process of the vector field is anisotropic. This anisotropy is parametrized by 
introducing two parameters p and q, which are determined by the conformal invariance 
breaking mechanism. If any of them are non-zero, generated £ is statistically anisotropic. 
Then the power spectrum of £ and the non-linearity parameter /nl have an angular 
modulation. 

This formalism is applied for two vector curvaton models and the end-of-inflation sce- 
nario. It is found that for p / 0, the magnitude of /nl an d the direction of its angular 
modulation is correlated with the anisotropy in the spectrum. If \p\ > 1, the anisotropic 
part of /nl is dominant over the isotropic one. These are distinct observational signa- 
tures; their detection would be a smoking gun for a vector field contribution to (. 

In the first curvaton model the vector field is non-minimally coupled to gravity and 
in the second one it has a time varying kinetic function and mass. In the former, only 
statistically anisotropic ( can be generated, while in the latter, isotropic ( may be realized 
too. Parameter spaces for these vector curvaton scenarios are large enough for them to be 
realized in the particle physics models. In the end-of-inflation scenario /nl have similar 
properties to the vector curvaton scenario with additional anisotropic term. 
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1. The Hot Big Bang and Inflationary 
Cosmology 



In the first Chapter of this thesis we start by reviewing briefly the Hot Big Bang (HBB) 
model of the Universe. It is one of the greatest achievements of the last century in un- 
derstanding the structure and evolution of the Universe from the first second until today 
13.7 x 10 9 years later. Predictions of the HBB model are in very impressive agreement 
with the observed distribution of the large scale structure and with the abundance mea- 
surements of the light elements. However, in addition to the dark matter and dark energy 
problems, the HBB model suffers from the need to fine tune initial conditions. The latter 
motivates us to look at the earlier stage of the evolution of the Universe. Currently the 
most popular and most predictive paradigm for this epoch is the inflationary scenario, 
which is introduced in section 1.5. 

1.1. Kinematics of HBB 

The HBB model relies on a hypothesis called the Cosmological Principle which states 
that the Universe is spatially homogeneous and isotropic on sufficiently large scales. 

The physical model of the Universe, is divided into two parts. One part describes 
the large scale behavior of the system and possesses high degree of symmetries so that 
mathematical models become simple and equations relatively easy to calculate. This 
is a background model. The second part deals with the deviations from the simplistic 
description of the background. These deviations are considered to be small compared to 
the background values. They don't influence the large scale behavior of the system: only 
the region much smaller than the scale on which background is defined. 

The Cosmological Principle is a hypothesis about the properties of the background 
distribution of matter in the Universe. The background is defined as the smeared-out 
distribution of matter, with smearing performed on large enough scales so that the dis- 
tribution appears smooth. However, a priory it is not clear that such scales do exist. It 
might be that probing larger and larger cosmological scales, we constantly discover new 
structures. This would happen if galaxies are distributed hierarchically at all distances 
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Figure 1.1.: Density fluctuations as the function of the size of the smoothing scale. The 
thick line represents a model with a scale invariant power spectrum and 
cosmological parameters Q m = 0.28, H = 72kms _1 Mpc~ 1 and r2f,/Sl m = 
0.16, r = 0.17, were £l m and is the dark matter and baryon density 
parameters respectively, H is the Hubble parameter in physical units and t 
is the optical depth [1]. 

as in the fractal Universe. In such a Universe probing larger and larger distances we find 
galaxies, clusters of galaxies, clusters of clusters of galaxies and so on. However, in the 
real Universe there is a scale at which the hierarchical structure stops and the Universe 
may be considered smooth. From the Figure 1.1 one can see that the Universe looks 
smoother and smoother if we probe it on larger scales. At around few hundreds Mpcs, 
which correspond to the size of largest superclusters, perturbations becomes smaller 
than the background value. At these scales separation of the matter distribution into the 
smooth background value and small perturbations is well justified. 

The isotropy hypothesis of the Cosmological Principle is supported by observations too. 
The strongest evidence comes from the measurements of the temperature irregularities of 
the Cosmic Microwave Background (CMB). The COBE satellite was the first to measure 
these irregularities [2] . They showed that the anisotropy in the temperature distribution 
is only of order AT/T ~ 10 -5 . In addition the evidence for the isotropy of the Universe 
is further supported by the galaxy redshift surveys, measurements of peculiar velocities 
of galaxies, distribution of radio galaxies, X-ray background and the Lyman-a forest [3]. 
Another assumption of the cosmological principle, the homogeneity of the Universe, is an 
inevitable conclusion if we assume the validity of the Copernican Principle. This principle 
states that our location in the Universe is not central or somehow special. Combined 



2 



1. The Hot Big Bang and Inflationary Cosmology 



with the evidence of the isotropy the outcome of the Copernican Principle is that the 
Universe is isotropic around every point. It can be shown that the last statement leads 
to the conclusion of spatial homogeneity. 

Accepting the validity of the Cosmological Principle we can find the metric for the ho- 
mogeneous and isotropic Universe. This can be done using only geometric considerations 
[4] giving the proper time interval as 



ds 2 = dt 2 - a z (t) 



d -' J +x 2 {de 2 + sm 2 



(1.1) 



1-Kx 2 

This metric is expressed in spherical coordinates and is called the Friedmann-Robertson- 
Walker (FRW) metric, t in Eq. (1.1) is the coordinate time and spatial coordinates 
l(t) = a (t) x are decomposed into the comoving coordinates x which are constant in 
time and the time dependent scale factor a (t) which parametrizes the evolution of the 
Universe, i.e. its expansion or contraction. In this metric, K parametrizes the curvature 
of space-time: if K < 0, the Universe is spatially open, if K > it is closed and if K = 
it is flat. As will be seen later, the inflationary paradigm predicts K ~ 0, which is in a 
very good agreement with observations. Therefore, in Chapters 2 and 3 we consider only 
the flat Universe in order to dispense with the unnecessary complications related with 
the curvature term. Furthermore, instead of using the spherical coordinate system in 
Eq. (1.1) in many situations it will be more convenient to use the Cartesian coordinate 
system. Then the flat {K = 0) FRW metric in Eq. (1.1) takes a simple form 

ds 2 = dt 2 -a 2 (t)dx l dx j . (1.2) 

1.2. Dynamics of the HBB 

From Eq. (1.1) we have seen that the evolution of the isotropic and homogeneous Uni- 
verse may be described by only one parameter, the scale factor a(t). To determine the 
dynamics of a (t) we have to specify the energy content of the Universe. In most sit- 
uations it may be well approximated by an ideal fluid whose energy-momentum tensor 
is 

T Ufl = (p + p)u fl u u +pg flu , (1.3) 

where is the velocity four-vector of the fluid, g^ v is the metric and p, p are the energy 
density and pressure of the fluid respectively. 

Using the FRW metric in Eq. (1.1) and the energy momentum conservation law 
V U T^ U = 0, where V y is the covariant derivative, we find 
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p = -3H(p + p), (1.4) 
where H is the Hubble parameter defined as 

H^- (1.5) 
a 

and the dot denotes the derivative with respect to the coordinate time t. Eq. (1.4) can 
also be rewritten as 

a^- = 3(p + p). (1.6) 
da 

For the perfect fluid the pressure is uniquely related to the energy density which is 
conveniently parametrized by the equation of state. Assuming a barotropic fluid we have 

p = wp, (1.7) 

where w is called the barotropic parameter. For different kinds of perfect fluid w will have 
different values, for example, for non-relativistic pressureless matter (sometimes called 
'dust') w = 0, for radiation (relativistic particles) w = 1/3 or w = —1 for the vacuum 
energy. Using the equation of state, from the continuity equation (1.6) it is easy to find 
the evolution of the energy density of the perfect fluid by integration: 

P = P °\ao) ' (L8) 

where '0' denotes initial values. Hence, it is clear that the energy density scales as p oc a~ 3 
for the pressureless matter, p oc a~ 4 for the relativistic matter and p = constant for the 
vacuum energy. 

To find how the content of the Universe determines the time evolution of the scale 
factor a (t) a theory of gravity must be assumed. For the purpose of this thesis it will 
be enough to consider only Einstein's theory of General Relativity (GR) with the field 
equation 

Rfiu ~ -^9^uR = mpiTpv, (1.9) 

where and R are the Ricci tensor and scalar respectively. Using temporary compo- 
nent of Einstein's field equation we can find the Friedmann equation. With the spatially 
homogeneous and isotropic metric in Eq. (1.1) it becomes 
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Furthermore, the acceleration of the Universe is obtained using the spatial components 
of Eq. (1.9). Together with Eq. (1.10) they give 

5 = -^?, (i.ii) 

which can be expressed in terms of the Hubble parameter as 

H + H 2 = _P±^P. (1.12) 
6m,p[ 

It is often useful to introduce the parameter O, which is related to the curvature of 
the space time [5]: 

0-1 = ^,. P.13) 

In the Einstein gravity this quantity measures the energy density of the Universe p relative 
to the energy density of the flat Universe p c for given H, called the critical energy density. 
This can be seen from the Friedmann equation (1.10): for zero curvature K = 0, the 
critical energy density is 

p c = 3m P1 iI 2 . (1.14) 

Plugging this back into the Friedmann equation and using Eq. (1.13), we find that in the 
Einstein's gravity 

o = A (i.i5) 

Pc 

If the energy density of the Universe is critical, = 1, the Friedmann equation becomes 

H 2 = -Ap. (1.16) 

3TOpj 

In accord with the comment above the equation (1.2) for the most of this thesis we 
consider only = 1 and the Friedmann equation of the form in Eq. (1.16). 

The early Universe is dominated by radiation, as can be seen from the scaling laws 
below Eq. (1.8). In this era it is useful to express Eq. (1.16) in terms of the temperature 
of relativistic particle species. To do this let us remember from thermodynamics that 
the energy density p of the weakly interacting gas of particles is given in terms of the 
internal degrees of freedom g^of an d its phase space distribution function / (p) as [6] 

p= WfS E{p)Hp)d3p > (li7) 

where p is the magnitude of the momentum of the particle and E is its total energy 
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E 2 = p 2 + m 2 . For particles in thermal equilibrium the distribution function / (p) is 



exp(^) ±1 



f(p) = 7^ > (1-18) 

T 

where [i is the chemical potential. The '+' sign here corresponds to the Fermi-Dirac 
species and ' — ' sign to the Bose-Einstein species. 

In the early Universe when it is dominated by the relativistic particles, i.e. T 3> m, 
we may take the limit T>/j, Inserting Eq. (1.18) into Eq. (1.17) and integrating it we 
obtain 

Pl = ~9*(T)T\ (1.19) 

where p 7 denotes the energy density of the relativistic particles and g* are the number 
of effectively massless degrees of freedom: 

9*(T)= 9i + ~ E 9i- (1-20) 

j=bosons i=fermions 

Note that g* is a function of the temperature because in this sum we included only 
relativistic species, i.e. particles with the mass m <C T. For example, at temperatures 
T <C MeV only photons and three neutrino species are relativistic giving g* = 3.36. At 
temperatures T > 300 GeV all particles of the Standard Model are relativistic resulting 
in g* = 106.75 [6]. 

Finally inserting Eq. (1.19) into Eq. (1.16) we find that in the flat, radiation dominated 
Universe the Hubble parameter is related to the temperature as 

fo~ T 2 

H = 7rJ^ — . (1.21) 
V 90 mpi v ; 

1.3. Big Bang Nucleosynthesis 

Arguably the biggest success of the HBB theory is the explanation for the origin of 
chemical elements in the Universe. According to this theory the lightest of them were 
created during the first three minutes after the Big Bang, when the Universe content 
was in a state of a very hot plasma. This process is called the Big Bang Nucleosynthesis 
(BBN). Due to its immense importance for the modern cosmology in predicting the 
abundances of light chemical elements and being a very sensitive method to constraint 
new theories of particle physics, in this section we give a summary of BBN. 

To describe the creation of light elements in the early Universe the crucial parameter is 
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the reaction rate T of some process under consideration. For illustrative purposes let us 
consider, for example, interactions of particles. Then T would represent the interaction 
rate per particle. The crucial quantity here is the ratio T/H, where H is the Hubble 
parameter. At the epoch of BBN the Universe is dominated by the matter which satisfies 
the strong energy condition, so that H^ 1 represents the age of the Universe. In this case 
T/H < 1 means that on average less than one particle interacted throughout the history 
of the Universe. In other words, we can say that particles are decoupled. If, on the other 
hand T/H > 1, particles have interacted many times and it is safe to assume that they 
are in thermal equilibrium. During the radiation dominated stage the expansion rate of 
the Universe is proportional to the temperature squared, H oc T 2 (see Eq. (1.21)), while 
the reaction rates are typically proportional to T oc T s . In the adiabatically expanding 
Universe the temperature decreases as T oc a . Hence, we can write T/H oc a 2 ~ s , from 
which we see that if s > 2 the process which was in equilibrium at some initial time, i.e. 
T/H > 1, it will fall out of equilibrium at later times. If the process we are interested 
is the interaction of particles, then we can say that after being in equilibrium, particles 
"freeze-out" when T/H becomes smaller than one, i.e. the number density is not affected 
by interactions. Then BBN can be roughly divided into three stages depending on which 
processes are in thermal equilibrium. 

When the temperature of the Universe was around 10 MeV, which corresponds to the 
age of 10 -2 s, the ratio of neutrons and protons is controlled by the weak interactions: 

n< — > p + e~ + v, n + v e < — > p + e~ , n + e + < — > p + v e . (1-22) 

where v e and v e are the electron neutrino and antineutrino, and e~ , e + are the electron 
and the positron. If the rate of these interactions are much more rapid than the expansion 
of the Universe, i.e. T nir ^ p /H > 1, the species involved in these interactions are in a 
thermal equilibrium, which means that the neutron to proton ratio evolves according to 

- = e-«/ T , (1.23) 
p 

where Q = m n — m p = 1.29 MeV is the mass difference of neutrons and protons. As we 
can see for the energies high above MeV, the number of neutrons and protons are almost 
the same, (n/p) ~ 1. 

The reaction rates for the processes in Eq. (1.22) can be calculated using Fermi theory 
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T (keV) 

4499 2107 100B 529 265 134 66 




Us) 



Figure 1.2.: Evolution of the n/p ratio. The solid red line represents the true variation 
while the dashed blue one represents an equilibrium evolution. BBN starts 
at T ~ 0.1 MeV which results in the steep decline of the red line at these 
energies. (Figure adapted from Ref. [7]) 



for the weak interactions, which gives [6] 

(~2G 2 F T 5 T»Q, m e , 

where r n is the neutron halflife and Gp = 1.1664 x 10 -5 GeV -2 is the Fermi con- 
stant. When the reaction rate T n< ^ p falls bellow the Hubble expansion rate, i.e. when 
rVi-t+p/H < 1, processes in Eq. (1.22) depart from the equilibrium and the number of neu- 
trons and protons "freezes-out". The approximate temperature of the freeze-out can be 
calculated using Eq. (1.21) and considering that T > m e : 



H I 1/6 0.8 MeV 



1.25) 



In the Standard Model of particle physics with the three (almost) massless neutrino 
species the number of relativistic degrees of freedom is g* = 10.75. Thus, the freeze-out 
temperature is found to be 

T fr ~lMeV, (1.26) 
which corresponds to about Is. The ratio n/p at this moment can be calculated from 
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Eq. (1.23) 

n 1 , 
- « -. 1.27 
p o 

However, this is not a true "freeze-out" because the n/p ratio is not constant but decreases 
slowly. This happens because of occasional weak interactions among neutrons, protons, 
and z/ e , v e eventually dominated by the free neutron j3 decay. However, this decrease 
is much slower than the equilibrium value given in Eq. (1.23) (see Figure 1.2). 

When the nucleosynthesis starts at about 0.1 MeV (corresponding to about 180 s), the 
neutron to proton ration had been decreased to 

I - \. ( , 28 , 

The main product of BBN is the helium-4, He. The production of heavier elements is 
very subdominant because there are no stable nuclei with the mass number 5 or 8 and 
hence no elements form through reactions such as n + 4 He, p + 4 He or 4 He + 4 He. In 
addition reactions such as T+ 4 He < — > 7+ 7 Li and 3 He+ He < — > 7+ 7 Be are suppressed 
because of the large Coulomb barriers. The formation of 4 He in principle could proceed 
directly through the four body collision. But the very low number densities of neutrons 
and protons renders this type of reactions negligible. Hence, the element formation must 
start with the production of deuterium through the two-body collision: 

p + ni — ^D + 7 . (1.29) 

Although the binding energy of the deuterium is Ad = 2.23 MeV, the formation of this 
element becomes effective only at much smaller temperatures. This is because of a large 
number of energetic photons which destroy deuterium. So D nuclei can start forming 
without being immediately photo-dissociated only when the number of such photons per 
baryon falls below unity, which occurs at the temperature T < 0.1 MeV [8]. Therefore, 
this period is called the deuterium bottleneck. But once deuterium starts forming, the 
whole set of reactions sets in producing other heavier elements. 

The final number density of 4 He depends on the whole nuclear network only very 
weakly. And it is a very good approximation to assume that all neutrons which didn't (3 
decay will end up being bound into the 4 He atoms. Hence, the helium mass fraction Y p 
can be calculated very easily just by power counting: 



2n 2 (n/p) 
n + p 1 + (n/p) 



Y p ~ -p!±- = 1 7';^ ~ 0.25. (1.30) 
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Number 


Reaction 


Number 


Reaction 


1 




9 


3 He (T, 7 ) 7 Be 


2 


p(n,i)d 


10 


T(T, 7 ) Y Li 


3 


D (p, 7) 3 He 


11 


7 Be (n,p) 7 Li 


4 


D (D, n) 3 He 


12 


7 Li (p, T) 4 He 


5 


D(D,p)T 


13 


4 He (D, 7 ) e Li 


6 


3 He (n,p)T 


14 


6 Li (p, T) 3 He 


7 


T (D, n) 4 He 


15 


7 Be (n, T) 4 He 


8 3 He(D,p) 4 He 


16 


7 Be (D,p) 2 4 He 



Table 1.1.: The most relevant reactions of BBN. Here, numbers of reactions correspond 
to the ones in Figure 1.3 (adapted from [8]). 
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The essential parameter for the processes of BBN is the number density of baryons. 
To quantify it, one usually uses the ratio of baryons over photons defined as 

mo = io lonB . (i.3i) 

At temperatures somewhat below T < 0.3 MeV, all the positrons have annihilated with 
the electrons and hence the number of baryons and photons in a comoving volume does 
not change. Therefore, 7710 must stay constant from BBN through recombination until 
today. And one can relate this value to the energy density parameter for the baryons 0,b 
today [9]: 

273A5n B h 2 /2.725K\ 3 /6.708 x lQ-^MeV' 2 " 



*° = 1 - 0.007Y, {-J^J \T G J' (02) 

where Tcmb is the photon temperature today and G is Newton's gravitational constant. 

The most relevant reactions for the BBN are summarized in Table 1.1 and Figure 1.3. 
The precise final abundances of each element, including 4 He, are calculated numerically, 
solving a system of coupled kinetic equations for each element as well as Einstein equa- 
tions, including the covariant conservation of total energy momentum tensor and con- 
servation of baryon number and electric charge. Results are given in Figure 1.4 together 
with observationally inferred values for some elements. 



1.4. The Problem of Initial Conditions of the Hot Big Bang 

The HBB cosmology is very successful in explaining the structure and evolution of the 
Universe after 1 s. However, in order to agree with observations the initial conditions of 
the HBB model have to be fine tuned. In this section we review briefly the problem of 
this fine tuning. 

1.4.1. The Flatness Problem 

Current observations agree very well with the density parameter O of the Universe being 
very close to one, i.e. the Universe is spatially flat. However, in the phase diagram the 
value Q = 1 is the unstable fixed point. In other words, any initially tiny departure from 
flatness will become larger and larger as the Universe evolves. This can be easily seen by 
using the Friedmann equation (1.6) and the definition of the f2 in Eq. (1.13): 
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Figure 1.4.: Abundances of light elements from the standard model of BBN. The bands 
represent 95% CL, the boxes represent observed values (smaller - ±2<r only 
statistical errors; larger - ±2a statistical and systematic errors), the narrow 
vertical column represents t]\q value inferred from CMB observations and the 
wider column indicates the BBN concordance range (both at 95% CL) [10]. 
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a 



Figure 1.5.: The evolution of the density parameter of the Universe filled with an ideal 
fluid whose barotropic parameter is w > —1/3. From the figure it is clear that 
any initial departure from critical density will grow in time, hence the spa- 
tially flat Universe is an unstable fixed point in the phase diagram (adapted 
from Ref. [11]). 



Assuming p corresponds to the energy density of a perfect fluid and using Eq. (1.8) we 
arrive at 

1 

n= l-[(O ini -l)M ni ]y3-+i' (° 4 ) 

where Q[ n i denotes the initial value of the density parameter and y = a/dQ. From this 
equation it is already clear that if w > — 1/3 any initial departure from the flat Universe 
with Qjnj 7^ 1 will grow in time. For example if initially the Universe is open, f2i n ; < 1, the 
energy density at any later time y > 1 will decrease monotonically towards zero Q — > 0, 
i.e. towards the empty Milne Universe. On the other hand, if initially the Universe is 
closed, fiini > 1, its energy density will increase rapidly and reach a singularity in a 
finite time. This behavior of the density parameter is illustrated in the phase diagram 
in Figure 1.5. Therefore, for the present Universe to be flat, f^o ~ 1, its initial energy 
density had to be extremely close to the critical value. For example, in order to reproduce 
the present Universe, the energy density at the time of BBN had to be 

I^bbn - 1\ < 1(T 16 . (1.35) 

It is extremely unlikely for ^bbn to be so close to unity by accident. 

The flatness problem of the HBB is sometimes rephrased as an age problem. This can 
be seen from Eq. (1.34) and Figure 1.5. If initially the Universe is closed Oj n i > 1, very 
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soon after its birth it recollapses without having time to form any galaxies and stars. If, 
on the other hand, initially the Universe is open f2; n i < 1, it soon becomes the empty 
Milne Universe before the formation of any structure. In both cases the Universe does 
not have time to become the one we observe today. 

1.4.2. The Horizon Problem 

In section 1.1 we have demonstrated that the isotropy and (with mild assumptions) 
homogeneity of the Universe at present as well as at the Last Scattering Surface (LSS) 
is an observationally established fact, which justifies the use of FRW metric. However, 
in the HBB model this fact is highly non-trivial due to the finite age of the Universe. 
We expect the space-time regions to be homogeneous and isotropic on scales which could 
have been in a causal contact, i.e. which could have "communicated" with each other. 
But because the maximum velocity of the signal is finite (equal to the speed of light in 
vacuum) and because the age of the Universe is finite too, there is only a limited distance 
at which two regions could have had a causal contact. 

To illustrate the horizon problem in HBB let us consider the epoch of the last scattering. 
From observations of the CMB we know that the fractional temperature variations at 
that time was of order 10 -5 . The distance of maximal causal contact at LSS is 

Causal ~ CtLSS, (1.36) 

where c is the speed of light and ^lss is the age of the Universe at LSS. At this epoch 
the size of the present horizon, given by Ztoday ~ citodayj was aLss/ a today times smaller 

^LSS ~ Cttoday— (1-37) 
Q-today 

Taking into account that the Universe was matter dominated at the last scattering and 
remained so until very recently, i.e. a cx we may compare £lss with the size of the 
causality length 

^LSS ^ / Qtoday \ ^ 3 
Causal V °LSS / 

Hydrogen recombined at the temperature Tlss ~ 1-6 X 10 5 K. Assuming adiabatic ex- 
pansion for the Universe T oc a , and using Tcmb ~ 2.7 K at a = atoday we find 



1.38) 



^causal / 



3 

~ 10 s . (1.39) 



^LSS \ 
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Figure 1.6.: Cosmic Microwave Background temperature variation map observed by the 
WMAP satellite [12]. 

Therefore, at the time of recombination the observable Universe consisted of at least 
10 5 causally disconnected regions with the fractional temperature variation of only 10~ 5 . 
No physical process could have caused such extreme smoothness in so many causally 
disconnected regions. This constitutes the horizon problem of HBB cosmology. 

1.4.3. The Origin of Primordial Perturbations 

As we have seen in Figure 1.1 the Universe can be considered isotropic and homogeneous 
only on smoothing scales larger that a few hundreds of Mpcs. On smaller scales it 
is highly inhomogeneous due to the presence of structures such as stars, galaxies and 
galaxy clusters. It is already established that this structure in the Universe formed due to 
gravitational instability, when slightly denser regions collapsed onto themselves forming 
a complicated web distribution of galaxies. However, for this process to be initiated the 
existence of some primordial seed density perturbations must be postulated. Indeed, the 
first observational proof of such perturbations was provided by the COBE satellite [2] (see 
Figure 1.6 for the high resolution CMB map from WMAP measurements). Unfortunately, 
the properties of these seed perturbations cannot be explained within the framework of 
HBB cosmology. 

Causality constraints require that seed perturbations could have formed only due to 
processes inside the causal horizon, which in the HBB model is monotonically decreasing 
as we go back in time. However, it was already realized in 60's and 70's that random dis- 
placements and movements of particles inside the horizon cannot produce the necessary 
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Figure 1.7.: The angular power spectrum of the CMB temperature-temperature (TT) 
and temperature-polarization (TE) anisotropics. Solid lines represent the 
best fit ACDM model [13]. 



16 



1. The Hot Big Bang and Inflationary Cosmology 



perturbation spectrum. Zel'dovich [14, 15] and Peebles [16] were the first to realize that 
such processes would produce matter density perturbations with the power spectrum 
P (k) oc k 4 . Such perturbations would result in an excessive overproduction of black 
holes on small length scales. Indeed, assuming a smooth power law spectrum of primor- 
dial perturbations, P (k) oc k n , from the CMB observations it was found that n ~ 0.96 
[17], very different from the n = 4 case. 

In addition, with more precise measurements of primordial perturbations other prop- 
erties became clear which cannot be explained by the standard HBB cosmology [18]. As 
seen in the upper graph of Figure 1.7 the CMB temperature power spectrum features so 
called "acoustic peaks". But more importantly this pattern is caused by adiabatic density 
perturbations. Such perturbations in the baryon-photon fluid upon horizon entry start 
oscillating only with excited cosine modes and with the same phase. Therefore, this is a 
strong indication that seed perturbations are present on scales larger than the horizon, 
i.e. they could not be created by causal processes during HBB. 

Although no viable model exists, in principle one could construct a model in which 
causal processes mimic the pattern of adiabatic acoustic peaks [19, 20, 21]. However, 
even stronger proof for the superhorizon origin of primordial perturbations is provided 
by the temperature-polarization cross correlation function [22]. The polarization signal 
is not affected on its path from LSS towards us. Therefore, by measuring polarization 
perturbations we can be certain to be probing the era of recombination. But as clearly 
seen in the lower graph of Figure 1.7 on angular scales 50 < I < 200 the temperature 
and polarization anticorrelates. Since these low multipoles represent superhorizon scales 
at LSS it is certain that primordial density perturbations were already present before 
entering the horizon. 

From this discussion, one can see that the problem of seed perturbations in the HBB 
cosmology in essence is a restatement of the horizon problem. The superhorizon origin 
of primordial perturbations is the strongest support for the inflationary scenario. 

1.5. Inflation 

1.5.1. The Accelerated Expansion 

The initial condition problems of standard HBB cosmology named in section 1.4 may be 
substantially alleviated if we postulate an accelerated expansion of the Universe at it's 
earliest stages. This epoch is called inflation. When we say "accelerated expansion" we 
mean that the distance between any two comoving points in the Universe is increasing 
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with a positive acceleration. In other words the scale factor in Eqs. (1.1) or (1.2) obeys 

a > 0. (1.40) 

This might be considered as the definition of the era when gravity is repulsive. 

Instead of Eq. (1.40) we may rewrite it in the form which gives more physical interpre- 
tation. As was discussed in section 1.2, H^ 1 defines the Hubble length. Then (aH)^ 1 is 
the comoving Hubble length. From Eq. (1.40) we find that 

1 (aH)- 1 < 0. (1.41) 

It shows that the comoving Hubble length during inflation decreases. Therefore, two 
points which initially were inside the Hubble radius (their comoving distance smaller 
than (aH)^ 1 ) at some moment goes outside this radius. This moment is called "the 
horizon exit". 

The condition in Eq. (1.40) can be written even in another form, which will be very 
useful in later sections. Substituting the derivative of the Hubble parameter H into 
Eq. (1.40) after some calculations we find —H < H 2 . When 



H 



<^H 2 , (1.42) 



the expansion is almost exponential. Even more so, if H — > 0, it is exactly exponential, 
i.e. a oc exp (Ht), and we call this de Sitter expansion. 

By postulating the early phase of accelerated expansion (Eq. (1.40)) of the early Uni- 
verse, the fine tuning problem of initial conditions of the HBB model discussed in sec- 
tion 1.4 are substantially alleviated. The crucial condition for this is Eq. (1.41). 

How the period of accelerated expansion solves the flatness problem can be seen by 
inserting Eq. (1.41) into Eq. (1.13). Because during inflation the comoving horizon is 
decreasing, a 2 H 2 grows with time and |fl — 1| is driven towards zero. Therefore, Q, = 1 
instead of being an unstable fixed point in the HBB model, becomes an attractor during 
the inflationary stage. To illustrate this let us use Eq. (1.34). As we will see shortly 
in the next section, the accelerated expansion of the Universe may be achieved if it is 
dominated by the vacuum energy for which p = —p, i.e. w = —1. Substituting the value 
w = — 1 into Eq. (1.34) we find that the density parameter approaches Q — > 1 and the 
phase diagram in the Figure 1.5 changes into the Figure 1.8. 

The decreasing comoving horizon size during inflation has an effect that initially any 
two causally connected points within the Hubble radius at some later moment leaves the 
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Figure 1.8.: The evolution of the density parameter for the Universe dominated by the 
vacuum energy with w = —1. Instead of being an unstable point in HBB 
model, Q = 1 becomes an attractor (adapted from Ref. [11]). 

horizon. This makes the causally connected regions to be larger than the horizon size 
and the postulate of the inflationary period solves a second problem of the HBB model. 

However, the most important achievement of inflation is that it explains the origin 
of superhorizon seed perturbations. According to the inflationary paradigm during ac- 
celerated expansion of the Universe vacuum quantum fluctuations are converted into 
classical perturbations. Details of this process with scalar field quantum fluctuations will 
be discussed in Chapter 2 and extended to the quantum fluctuations of vector fields in 
Chapter 3. But before that let us discuss what may cause an accelerated expansion of 
the Universe. 

1.5.2. The Scalar Field Driven Inflation 

There are several reasons why the Universe could have expanded exponentially. It might 
be that, at the relevant energy scales, the Einstein gravity is not a viable theory of Nature 
and it must be modified. Modification is such that it gives almost exponential expansion 
of space-time. The very first proposed model of inflation was due to this kind of modified 
gravity theory [23]. 

Another possibility is that inflation happens at the energy scales where Einstein gravity 
is still a viable theory of nature. Then Eq. (1.40) puts constraints on properties of 
matter which may be responsible for the inflationary expansion. This can be found using 
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Eq. (1.11): 

p + 3p<0. (1.43) 

Because the energy density p must always be positive, it follows that the Universe un- 
dergoes accelerated expansion if the pressure is negative enough, p < — p/3. The lower 
bound for pressure is determined by the dominant energy condition, which requires that 
p > —p. Should this bound be violated the propagation of energy outside the lightcone 
becomes possible and one cannot guarantee the stability of the vacuum [24]. Taking 
the extreme case p = —p and from Eq. (1.12) we find that H — > and the Universe is 
expanding exponentially (de Sitter Universe). 

The equation of state p ~ — p may be realized in the framework of GR if the Universe 
is assumed to be dominated by classical scalar fields. For simplicity we will assume 
that only one such field is relevant, which is then called the inflaton. The most general 
Lagrangian which is consistent with the GR for a single field inflation is given by [25] 

C = P(X,<ft), (1.44) 

where (ft is the scalar field, X = ^d^cftd^cft and P is some function. Inflationary models 
which study the evolution of the Universe under the influence of such fields are called 
k-inflation. But to make essential properties of inflationary models more transparent 
let us concentrate on a particular case where the field is canonically normalized. Then 
Eq. (1.44) becomes 

C= l -d^<p-V{4>), (1.45) 

where V ((ft) is the potential. The equation of motion of the field is obtained by requiring 
that the variation of the action with respect to the field vanishes. For the homogeneous 
component this gives 

eft + 3H(j> + = 0, (1.46) 

where V& = dV ((ft) /deft is the derivative of the potential with respect to the field and 
because we are interested in the homogeneous part of the field we have neglected gradient 
terms. 

The energy-momentum tensor of a scalar field may be obtained using the variation of 
the action with respect to g^ u [26] 

V = -2^p7+0i«X. (1-47) 



20 



1. The Hot Big Bang and Inflationary Cosmology 



For the scalar field with the Lagrangian of the form in Eq. (1.45) it becomes 

"1 



~ V (0) 



1.48) 



From this equation one can notice that the energy-momentum tensor for the homogeneous 
scalar field (homogenized by inflation) becomes as that of the perfect fluid so that the 
energy density p and pressure p can be defined as 



P = ^-V^). 
Using the first of these relations the Friedmann equation in Eq. (1.10) becomes 

3m|!# 2 = + V. 
Differentiating it, we derive another useful expression 



Smpjiir = — 4> 2 , 



1.49) 
1.50) 

1.51) 
1.52) 



were we have used the equation of motion in Eq. (1.46) as well. 

From expressions (1.49) and (1.50) it is clear that the condition p ~ — p is satisfied if 
the kinetic energy of the field is negligible compared to the potential one, i.e. 



H 



H 2 



< 1 ^ <t> 2 <.V ((p) . 



1.53) 



This requirement is called "the slow-roll condition" and is fulfilled if the potential of the 
field is sufficiently shallow. These conditions may be conveniently rewritten in terms 
of slow-roll parameters e and r\ in the following way. Because the field is slowly rolling 
we might also expect that the second derivative is also small, <p/H -C (j). Then, in the 
equation of motion (1.46), the first term is negligible 



3iT0 ~ —Vs- 



1.54) 
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The first slow-roll parameter becomes 




(1.55) 



On the other hand, differentiating Eq. (1.54) and using conditions in Eq. (1.53) we define 
the second slow-roll parameter 

V = ™pi^ « 1- (1-56) 

Conditions in Eq. (1.55) and (1.56) are called flatness conditions for the shape of the 
potential of the scalar field. The quasi exponential expansion during inflation lasts as 
long as those conditions are satisfied. When any of the parameters e of n becomes of 
order one, inflation ends. 

It is convenient to define a number of e-folds until the end of inflation 

N = i n ^*5°d) = f tcnd Hdl = [ 2d0 = m p , 2 I ^-d<p. (1.57) 

If the initial time in this equation is chosen to be when cosmological scales leaves the 
horizon, we may calculate how many numbers of e-folds of inflation is needed to solve the 
horizon and flatness problems of the HBB model. If the energy scale of inflation and the 
reheating temperature in Eq. (1.62) is at the supersymmetry energy scales, then N ~ 60 
[5]. In some models, with very low reheating temperature it can go down to N ~ 40. 
If we assume the validity of Einstein's gravity after inflation and fields with canonical 
kinetic terms the maximum number of e-folds is N ~ 70. 

1.5.3. The End of Inflation and Reheating 

Inflation ends when the slow-roll parameters defined in Eqs. (1.55) and (1.56) become 
of order one. Soon after this happens, the inflaton field fast-rolls towards its VEV and 
starts oscillating around the minimum of its effective potential. Expanding this potential 
around the minimum, the leading term in the series is V {(/)) = \m 2 {(/) - ((f))) 2 , where (0) 
is the VEV and m is the mass of the inflaton field. Inserting this into Eq. (1.46) we see 
that for m S> H the field acts as the underdamped harmonic oscillator with the frequency 
oj = m, much larger than the Hubble time. Therefore, in accord with the equation of 
motion of the harmonic oscillator we may write (f) 2 = 2V ((f)), were the average values 
are defined over one Hubble time. Inserting this into Eq. (1.50) we find that the average 
pressure of the oscillating scalar field is p = and the energy density from Eq. (1.49) is 
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p = (f> 2 . Using this, Eq. (1.46) may be rewritten as 

p + 3Hp = 0. (1.58) 

Taking into account that the Universe is dominated by the oscillating inflaton field with 
the zero pressure from this equation we obtain 

p = and p oc a~ 3 . (1.59) 

Therefore, the oscillating inflaton field acts as the pressureless matter and the Universe 
evolves as dominated by the non-relativistic dust particles (inflatons) [27]. 

Because the field is oscillating it might be interpreted as the collection of massive 
inflaton particles with zero momentum. Before inflaton oscillations the temperature of the 
Universe is effectively zero. However, for the successful BBN, discussed in section 1.3, the 
Universe must be radiation dominated with the temperature above lOMeV. Therefore, 
to recover the successes of the HBB cosmology, the energy stored in the inflaton field 
must be released to effectively massless particles. This process is known as 'reheating'. 
The first proposals for the mechanism to reheat the Universe were based on the single- 
body decays [28, 29]. During inflation such decays may be neglected because the field 
is not oscillating and cannot be interpreted as a collection of particles. But during the 
phase of coherent oscillations inflaton particles may decay into other scalar particles x 
or fermions ip through the terms in the Lagrangian such as gcftx 2 an d hcfitpip, where g 
is the coupling constant with the dimension of mass and h is a dimensionless coupling 
constant. Due to these couplings the equation (1.58) must include an additional friction 
term T which parametrizes the inflaton decay into these particles 

p+{3H + T)p = 0, (1.60) 

where T = r^_> xx + T^^. When the mass of the inflaton is much larger than those of 
X and ip, i.e. m 3> m x ,m^, the decay rates are known to be [29, 30] 

When H > T the number of produced particles is very small (see section 1.3) and 
they do not influence the dynamics of the Universe. However, these particles may still 
thermalise and their temperature becomes much larger than the temperature at reheating 
(given in Eq. (1.62)) [31]. At time i re h, when the Hubble parameter becomes H ~ T, 
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the decay processes become significant and practically all inflaton energy is transferred 
to the newly created particles. The temperature of the Universe at this moment may be 
calculated using the flat Friedman equation in Eq. (1.16) and assuming that new particles 
are relativistic, then from Eq. (1.19) we get 



where = 10 2 — 10 3 [32] is the number of effective relativistic degrees of freedom defined 
in Eq. (1.20). 

The mechanism of reheating described above is based on perturbative particle decay. 
However, in some inflationary scenarios the energy transfer from the inflaton field may be 
preceded by another, much more efficient process. To distinguish it from the conventional 
reheating, it is called 'preheating'. In the first such proposal, the parametric preheating, 
the inflaton field decays into relativistic particles of other fields very rapidly in short, 
explosive bursts due to the parametric resonance effects [33, 34]. At the second stage, 
these particles decay into relativistic species which finally thermalise. It should be noted, 
however, that it is not possible to transfer the total energy stored in the inflaton field by 
this process. When the amplitude of inflaton oscillations decreases below some critical 
value, the parametric resonance becomes inefficient. The residual oscillating inflaton 
field must decay through the perturbative reheating processes described above. If these 
processes are not efficient enough, due to the scaling law in Eq. (1.59), the residual 
oscillating inflaton field comes to dominate the relativistic decay products of preheating. 
In this situation the transfer of the inflaton energy into radiation is still dominated by 
the perturbative reheating processes. 



^reh - 9- 



C 1/4 v^, 
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2.1. Statistical Properties of the Curvature Perturbation 

2.1.1. Random Fields 

As it will become clear in section 2.2 the origin of cosmological perturbations is quan- 
tum mechanical. But quantum mechanical processes are non-deterministic: one can only 
predict the probability of experimental outcome. Therefore, to make quantitative de- 
scriptions of these processes one needs to use statistical methods. The same is true for 
cosmological perturbations. One cannot calculate exact values of perturbations at each 
space point, only the statistical properties may be predicted by theories and compared 
with observations. To quantify the properties of cosmological perturbations a very useful 
method is to describe them as random fields. 

Let us introduce some random field j3. It is assumed that our Universe is just one 
realization of many (hypothetical) possible universes. Then, to each of these universes 
one can assign a particular realization (3 n from the whole ensemble Depending on the 
problem to be solved, functions (3 n may parametrize, for example, the spatial distribution 
of the density, velocity or other fields. Each of the functions (3 n are realized with the 
probability p (f3 n ) dn, where n is a continuous index and p is the probability distribution 
function (PDF). 

Properties of the random field (3 are specified by the form of PDF. It is said that the 
random field f3 is statistically homogeneous if the probability p (f3 n ) of the realization f3 n 
is the same as that of realization f3 m , where (3 n (x) = /3 m (x + X) V X. In other words, 
probabilities are equal for realizations which differ only by the spatial translation. And 
(3 is said to be statistically isotropic at a point x if probabilities are equal for realizations 
which differs only by rotation, i.e. p (f3 n ) = p (/3 m ) for (3 n (x) = (3 m (IZx.), where 1Z is the 

1 As is usual in the literature the notation /3 n (x) is used to denote two things: a function itself and the 
value of that function at the point x. We will adopt the same notation here hoping that the meaning 
will be clear from the context and no confusion will arise. In addition, to denote a function itself 
(not it's value) we will use (3 n too, keeping in mind that it is a function of the spatial argument x. 
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rotation matrix, 1Z G SO (3). Analogously, (3 is parity conserving if p((3 n ) = p((3 m ) for 
/3 n (x) = /3 m (— x) V x. For the following discussion we will consider only statistically 
homogeneous and parity conserving fields. Usually in cosmology it is assumed that the 
field of the primordial density perturbation is statistically isotropic as well. But, as we 
will show in Chapter 3 this might not necessarily be so. 

Instead of working with PDF of the random field directly more convenient and obser- 
vationally more relevant quantities are iV-point correlation functions. For example the 
two-point correlation function is related to the PDF as 



Integration over n shows that it is the average over all the ensemble. In general the 
two-point correlation function does not specify the PDF uniquely, one needs to calculate 
higher order correlators which are defined analogously. 

A very powerful way to analyze correlation functions is by decomposing them into the 
eigenvectors of the translation operator. In flat space this corresponds to the decompo- 
sition into Fourier series. But to perform this decomposition it is necessary to chose the 
box of a certain size with periodic boundary conditions. In the cosmological context the 
choice of the box size is a very important issue. One requires that the box is large enough 
so that wave-vectors k could be treated as continuous and the Fourier series could be 
replaced by an integral. On the other hand, it is undesirable that the box is infinitely 
large. It might be that at very large distances the Universe becomes very anisotropic 
and inhomogeneous. This for example happens in chaotic inflationary models. Therefore, 
choosing too large a box one would have to take into account unknown physics. Usually 
it is enough for the box size to be only several orders of magnitude larger than the horizon 
of the observable Universe, so that m{HQL) ~ O (1), where L is the comoving box size 
and Hq is the Hubble parameter today. Such a box is called a minimal box [5, 35]. This 
choice is sufficient to approximate Fourier series as integrals. And we normalize Fourier 
modes such that 



Because (3 n (x) describes the distribution of real quantities in the Universe, they must 
be real functions themselves. This translates into the requirement that imaginary Fourier 
modes (3 n (k) must satisfy the reality condition j3 n (— k) = /3* (k). We note as well, that if 
the random field f3 is statistically isotropic, then f3 n (k) does not depend on the direction 
of the wave-vector k, only on it's modulus k, i.e. f3 n (k) = f3 n (k), where k = |k|. 




(2.1) 




(2.2) 
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If the random field is invariant under spatial translations, i.e. if it is statistically 
homogeneous, then the Fourier transform of the two-point correlator in Eq. (2.1) is 
determined by the reality condition 

(P n (k a ) (3* (k 2 )) = (2tt) 3 5 (kj - k 2 ) Pp (k) , (2.3) 



where Pp (k) = (|/3 n (k)|^ is called the power spectrum (remember that (...) means 
the ensemble average) . Note that the presence of the delta function in this expression is 
the result of statistical homogeneity of the random field. This relation can be rewritten 
using p n (-k) = ft* (k) as 

</3 n (kx) p n (k 2 )> = (2vr) 3 S (ki + k 2 ) Pp (ki) . (2.4) 

The power spectrum (k) is related to the two-point correlation function in the 
position space by the Wiener-Khinchin theorem. This theorem states that Pp (k) is the 
Fourier transform of the latter 

Pp (k) = j (P n (x) p n (x + r)) e- ik ' r dr. (2.5) 

It is often convenient to use another definition of the power spectrum which differs 
from the first one just by normalization 

Vp (k) = ^Pp (k) . (2.6) 

Both of these definitions have to satisfy the reality condition, i.e. Vp (— k) = Vp (k). 
For the future convenience we will parametrize the directional dependence of the power 
spectrum as [36] 



Vp^) = vr(k) 



i + //(d-k) 2 + ... 



(2.7) 



where Vp° is the average over all directions, d is some unit vector, k is the unit vector 
along k and k = |k| is the modulus of k. 

The meaning of the power spectrum Vp can be easily understood in case of statistically 
isotropic perturbations, i.e. when Vp (k) = Vp (k). Then from the inverse of Eq. (2.5) 
we find that the variance of the random field /3 is equal to 

^(x) = </3 2 (x)> = — 3 / Pp(k)d?k. (2.8) 
(Ztt) Jo 
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Since for statistically isotropic perturbations Pp (k) depends only on the modulus of k 
it is convenient to express this integral in spherical coordinates. Then the definition in 
Eq. (2.6) can be rewritten as 

<r}= P p (k) dmk. (2.9) 
Jo 

Therefore, Vp (k) corresponds to the contribution to the variance cr? per logarithmic 
interval in k. And because we assumed statistical homogeneity of /3, the variance does 
not depend on position. 

If the power spectrum Vp is scale independent then the integral in Eq. (2.9) is loga- 
rithmically divergent. Divergences for large and small k in this integral are avoided by 
introducing cutoff scales. For large k the cutoff scale R s corresponds to the smoothing 
scale and for small k (large spatial distances) i?box corresponds to the maximum size of 
the box in which we perform calculations 

o} = Vpd\nk = Vpln^. (2.10) 

box 

With the minimal box size, such that In (Rbox/Rs) is of order one, the mean-square is 
roughly of the order of the spectrum. 

If Eq. (2.5) is to be applied in the cosmological perturbation theory it requires an 
additional assumption. In practice we can observe and make measurements only of one 
Universe. Hence, the ensemble average over one Universe does not make sense and we 
cannot use this equation directly. To connect theoretical predictions with observations 
we have to assume the validity of ergodicity for our Universe. This assumption states 
that the average over the whole ensemble of universes is equivalent to the spatial average 
over one universe. To see what this means in mathematical language let us write the 
spatial average of the product of two points over the universe of realization f3 n 

J n (x) J n (x + r) = L- 3 J (3 n (x) (3 n (x + r) dx, (2.11) 

where L~ 3 is the box over which the averaging is performed. Then ergodic assumption 
states that in the limit L — > oo 

</3(x)/3(x + r)> = /3 n (x)/3 n (x + r). (2-12) 
As one can see, this assumption relates averages over the all ensemble of universes, 
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which cannot be measured, to the average over one universe, which can be measured. For 
Eq. (2.12) to be strictly valid we required an infinite box over which the measurement 
is performed. Of course this cannot be realized practically. The effect of the finite box 
introduces the so called 'cosmic variance' - when the separation between points in the 
correlators approaches the size of the box, the probability that the spatial average differs 
from the ensemble average increases. 

Until now we have considered only the two-point correlation function of Eq. (2.1) which 
is demanded by the reality condition. If the random field /3 is Gaussian, this correlator 
specifies PDF completely. Which means that the three-point correlator vanishes, while 
the four-point correlator can be expressed as the sum of two-point correlator products 
and so on. In the non-Gaussian case, the random field has a non-vanishing three-point 
correlator and the four-point correlator has additional terms which cannot be reduced to 
the product of two-point correlators. Let us limit ourselves only up to the three-point 
correlator. Although in cosmological context for some models higher order correlators 
might be as important as the three-point correlator, for the scope of this thesis the 
three-point correlator will be sufficient. It can parametrized similarly to Eq. (2.4) as 

{P (ki) (k 2 ) (k 3 )) = (2tt) 3 5 (kx + k 2 + k 3 ) Bp (kj, k 2 , k 3 ) , (2.13) 

where Bp is called the bispectrum. 

2.1.2. The Curvature Perturbation and Observational Constraints 

In the previous section we discussed random fields in general. Let us now turn to the 
discussion of the curvature perturbation £ which will be the main topic for the rest of 
this thesis. As was explained in section 1.5 the largest achievement of the inflationary 
paradigm is that it predicts the statistical properties of the curvature perturbation which 
can be compared with observations. 

Usually observational constraints on the statistical properties of £ are obtained with 
the assumption of statistical isotropy. However, one would expect that the presence of 
anisotropy at 10% level or so would not alter the results significantly. The strongest 
constraints on £ comes from the measurements of the CMB and large scale structure 
which probe the range A In A: ~ 10 [17]. The largest probable scale corresponds to the 
size of the observable Universe, k~ l ~ HZ 1 . 
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2.1.2.1. The Power Spectrum 

The shape of the power spectrum is the primary tool to contrast predictions of the 
inflationary models with observations. To quantify this shape the power spectrum is 
parametrized as 

/ u \ n(fc )-l+§n' 

V( . (fc) = p c (fc ) (_ J , (2.14) 

where ko = 0.002 Mpc -1 is the pivot scale, n is called the spectral index and parametrizes 
the scale dependence of the power spectrum and n' = dn/d In k is the running of the spec- 
tral index. Such parametrization is sufficient because according to observations n' -C n , 
thus higher derivatives are even smaller and can be neglected. Of course with such simple 
parametrization one looses sensitivity to the sharp features of the power spectrum. But 
according to some investigations (e.g. Ref. [37]) such features are not detected. The 
normalization of the power spectrum (ko) was first measured by the COBE satellite 
and most recently by the WMAP [17]. The present value is 

p c (k ) = (2.445 ± 0.096) x 10 -9 , (2.15) 

where this and later intervals are given at 68% CL. 

For the simplest, scale invariant case, called Harrison-Zel'dovich or flat power spectrum, 
n = 1 and n' = 0. However, according to current observations the spectral index is 3.1 
standard deviations away from the Harrison-Zel'dovich one. Indeed, n is smaller than 1. 
Such power spectrum is called red. With the assumption of negligible running, n' = 0, 
and no gravitational waves the spectral index is determined as 

n = 0.960 ±0.013. (2.16) 

If the running of the spectral index is allowed then this constraint is relaxed 

n = 1.017 ±0.043 and ri = -0.028 ± 0.020, (2.17) 

where gravitational wave contribution is still neglected. Letting non-negligible contri- 
bution from the gravitational waves relaxes these bound even further. However, gravi- 
tational waves are not observed yet, and as was mentioned earlier, in this thesis I will 
assume that their contribution is negligible. 

There are several reports of the detection of the angular modulation of the power 
spectrum in Refs. [38, 39]. Ref. [38] determined the modulation amplitude g defined in 
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Eq. (2.7) as 



g = 0.29 ±0.031 



(2.18) 



at 9<T confidence level. This is a definite proof of the existence of the preferred direction 
in the power spectrum. However, these two works also show that this direction is very 
close to the ecliptic poles, with the galactic coordinates (I, b) = (96,30). This is a 
very strong indication that the origin of the observed anisotropy is not cosmological but 
most probably caused by some systematic effects or comes from within the solar system. 
Although Ref. [38] have investigated the known systematic effects, including the Zodiacal 
light, but they could not find any which reproduces the observed signal. 

Given the above value of g we may place an upper limit on the anisotropy in the power 
spectrum of the cosmological origin. In this thesis we will assume that the upper bound 
on g in the primordial power spectrum is 



2.1.2.2. The Bispectrum 

Although the shape of the two-point correlator or it's counterpart in the Fourier space, 
the power spectrum, provides a very valuable information in discriminating inflationary 
scenarios and constraining physics of the early Universe, it has a limited potential. There 
are plenty of different inflationary models which predict similar power spectrum. Very 
powerful additional tools for distinguishing these models are higher order correlators. 
The Fourier transform of the three-point correlator is called the bispectrum and was 
defined in Eq. (2.13). While only two points are cross-correlated to obtain the power 
spectrum an infinitely more configurations are possible by cross-correlating three points. 
Therefore, the amount of information stored in the bispectrum is immensely richer than 
in the power spectrum, provided the curvature perturbation is non-Gaussian. 

However, as will be seen in section 2.4.1 single field, slow-roll inflationary models pre- 
dict negligible non-Gaussianity of the curvature perturbation. Observationally interesting 
non-Gaussianity can be generated only if any of the single field slow-roll assumptions or 
some combination of them are violated. These can be classified into four classes [40]: 
1) single free field, 2) canonical kinetic energy, 3) slow roll and 4) initial Bunch-Davies 
vacuum. In the first case large non-Gaussianity can be present if the curvature pertur- 
bation is generated by the different field from the one which drives inflation (two of such 
mechanisms are discussed in sections 2.4.2 and 2.4.3) or in the multifield inflation where 
the curvature perturbation is generated by many fields which drives inflation. In addition 



<?<0.3. 



(2.19) 
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if the inflaton cannot be considered as a free field, interaction terms can produce large 
non-Gaussianity as well. The second condition is violated for example in the k-inflation 
models [41]. In these class of models the speed of sound is different from the speed of light. 
The third condition might be violated if, for example, the inflaton potential has some 
sharp features which result in temporally violation of slow- roll conditions. The fourth 
assumption considers the initial fluctuations of the field. Usually it is assumed that ini- 
tial quantum fluctuations correspond to the Bunch-Davies vacuum (see section 2.2.2.4) 
which results in Gaussian statistics. If, due to some quantum gravitational effects, the 
initial quantum state does not correspond to the Bunch-Davies vacuum, field perturba- 
tions may be non-Gaussian and this non-Gaussianity will be translated into the statistical 
properties of the curvature perturbation 

Usually non- vanishing three point correlator of the curvature perturbation is paramet- 
rized by the non-linearity parameter /nl- There are several definitions of /nl in the 
literature. I will use the one which coincides with the definition used by WMAP team 



where 'c.p.' stands for 'cyclic permutations'. 2 

The simplest form of non-Gaussianity is of the local type which can be written as 



where £ g is the Gaussian part with zero mean, = 0. 

The strongest constraints on /nl comes from the measurements of the CMB sky. If the 
non-Gaussianity is of the local type in Eq. (2.21), then from WMAP5 data the constraint 
with 95% confidence level (CL) is (Ref. [17]) 



In this expression 'local' means the 'squeezed' configuration where one momentum is 
much smaller that the other two, k± ~ &2 3> k%. In the equilateral configuration with all 
three momenta of the same size, h\ = k<i = ks, the constraint on /nl is weaker 




(2.20) 




(2.21) 



C 9 (x) + ^/nl (C g 2 (x) - <C g 2 » 



9 < /nT 1 < HI- 




- 151 < f\ 



•equil 



< 253 



(2.23) 



NL 



2 The factor 6/5 comes from the fact that during matter domination, which is the case at the era of 
decoupling, the Newtonian potential $ is related to the curvature perturbation by <!> = |£. 
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at the same CL. 

The bounds on the magnitude of /nl will improve substantially in the very near future. 
If it is not detected by the Planck satellite, the constraints will reduce to |/nl| ^ 5 at 
95% CL, which is very close to the limit of an ideal experiment of |/nl| ~ 3 at 95% CL, 
limited by the cosmic variance [42]. The above bounds are given with the assumption 
that /nl is isotropic. Ref. [43] analyzed the WMAP5 data for the angular modulation 
of /nl^- However, due to the large measurement errors no conclusive statement can be 
made. 

2.2. Scalar Field Quantization 

In this section we discuss the quantization procedure of quantum field theory (QFT) in 
flat space-time (FST) and then generalize this formalism to curved space-time (CST). 
The discussion is solely about quantization of scalar fields, because they are the most 
simple ones and help to highlight the underlying principles. The extension to vector 
fields will be given in Chapter 3. 

Quantum mechanics was firstly formulated in the so called Schrodinger picture in which 
operators are time independent and state vectors evolve according to the Schrodinger 
equation. Equally well one can formulate this theory in the Heisenberg picture, where 
state vectors are constant but operators are changing with time. Quantum field theory 
can be formulated in both of these pictures as well, but this is much easier done in the 
Heisenberg picture, where operators are time dependent and satisfying field equations. 
Hence, the name quantum field theory. 

2.2.1. Quantization in Flat Space-Time 

Field quantization in FST may be presented in two ways [44]. In the first one we con- 
sider a classical field theory. Expand the field in Fourier modes and find that Fourier 
coefficients obey the equation of harmonic oscillator. With every harmonic oscillator we 
associate a position variable and the conjugate momentum in field space. The classi- 
cal harmonic oscillator is then first quantized. This is done by substituting c-numbers 
(classical numbers) of the position and momentum to the q-numbers (quantum numbers) 
and imposing commutation relations which are the result of the Heisenberg uncertainty 
principle. Then one finds that the Fourier coefficients (which are now operators) cor- 
respond to the raising and lowering operators in the Fock space, which are commonly 
called creation and annihilation operators respectively. 
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Another approach is to quantize the degrees of freedom of the classical field directly. In 
this approach one identifies the degrees of freedom of the field, finds their conjugate pairs 
of variables, changes them into q-numbers and imposes the same commutation relations 
as in the previous case. Only after quantization do we resort to the Fourier series. Again, 
we find that Fourier coefficients correspond to creation and annihilation operators. 

Results of both methods are the same. Although the first method is more intuitive and 
easier interpretable, the second method is more directly generalizable to CST. Hence, in 
this section we will take a standpoint of the second method in order to present the FST 
formalism in a way which is directly generalizable to the CST case. 

In the classical field theory equations of motion (EoM) for fields are obtained using 
the least action principle. Forming the action as 

S(4>) = J Cifad^tfx, (2.24) 

the classical field equations are calculated by requirement that the variation of the action 
should vanish 

' )S 0, (2.25) 



84>i (x) 

where x = (t, x) and <f>i (x) are classical fields. In principle these fields could be complex 
and after quantization we would find that they describe pairs of particles and antipar- 
ticles, i.e. the field would have a charge. But in context of producing the curvature 
perturbation in the Universe we are interested only in the neutral particles, which agrees 
with observations of the neutrality of the Universe. Therefore we will be interested only 
in real fields 4>i (x). 

Let us start with a free, massive, real scalar field. The relativistically invariant La- 
grangian for such field is written as 

£ = ^W-^mV, (2.26) 

where m is the mass term. The variation of the action in Eq. (2.24) with this Lagrangian 
gives the familiar Klein-Gordon field equation for the relativistic field 

[d^ + m 2 } <j> = 0. (2.27) 

The general solution for this equation can be written as the superposition of the com- 
plete set of orthonormal solutions, {u a (x)}. Where orthonormality is defined through 
the scalar product. For the Klein-Gordon equation in FST the scalar product of two 
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wave functions is 

(u m ,u n ) = i I u* n dou m dx = i J {u* m d u n - u n d u* m ) dx. (2.28) 
Then the complete set of orthonormal solutions {u a (x)} must satisfy 

(u m ,u n ) = 5 mn , (u* m ,u n ) = -b mn and (u m ,u n ) = (u* m ,u n ) = 0. (2.29) 

Indices m and n can be discrete or continuous. In the latter case Kronecker symbols 6 m n 
should be replaced by Dirac delta functions. In these expressions {u m (x)} and {u^ (x)} 
denotes a complete set of positive and negative frequency solutions respectively. Using 
these sets of solutions the general solution of Eq. (2.27) may be written as the sum of 

{u m , Urn}'- 

( X ) = [ amUm ( X ) + a \n u *m i x ) i (2.30) 

m 

where coefficients a m are given by 

a m = {(f>,u m ) (2.31) 

In the classical field theory cp (x) actually describes an infinite number of degrees of 
freedom at each space point x. One can find a conjugate momentum for each of these 
degrees of freedom by using equation 

7T = 6 -i = 4>. (2.32) 
dcp 

In this way for each spatial point x we prescribe a generalized coordinate variable, 
(f)(x), and a conjugate momentum, tt(x). Field quantization proceeds by analogy with 
quantum mechanics, which is changing c-numbers <p {x) and ir (x) into q-numbers (f) (x) 
and 7T (x) and imposing commutation relations for them 

(t,x) ,tt (t,x') =i<5 3 (x-x / ), $ (t, x) , $ (t, x') = [vr(t,x) ,tt (t,x)] = 0. (2.33) 

Because the field variable 4>{x) was promoted into the operator (j)(x), the expansion 
coefficients in Eq. (2.30) have to be operators as well, i.e. the substitution a m —> a m 
must be made. And commutation relations for these coefficients may be calculated from 
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Eqs. (2.33) [45, 46, 47]: 



0"mi Q"n 



fimn and [a m , a r 



0. 



(2.34) 



Operators d m and a\n are interpreted as the rising and lowering operators in Fock 
space, or creation and annihilation operators respectively, 

a]^ \n m , {u}) = ^n m + 1 \n m + 1, {«}) and a m |n m , {u}) = ^/n^ \n m - 1, {u}) , 

(2.35) 

where n m is the number of particles in a state m, notation of {u} in the ket reminds us 
that the definition is for particular complete set of orthonormal mode functions {u m }. In 
these equations coefficients \/n m + 1 and chosen for the correct normalization 

of the vacuum state (0, {u}\ 0, {u}) = 1, where the vacuum of this Fock space is defined 
as 

a m \0,{u}} = 0. (2.36) 
For the following discussion it will be useful to introduce an operator N such that 



(2.37) 



The meaning of this operator becomes clear when we take the expectation value 

(?1 m | N m |^m)) 

(n m \ N m \n m ) = (n m | a} m a m \n m ) = n m . (2.38) 



Hence N m can be interpreted as the number operator of m particles. 



2.2.1.1. Interpretation of a m and a^ m 

We have mentioned that operators a m and a} m are interpreted as creation and annihilation 
operators. What justifies such interpretation? To show this let us find mode functions 
{u m } explicitly. 

First of all {u m } must represent particles with positive energy, therefore these functions 
must be positive frequency solutions of Eq. (2.27), where positive frequency is defined 
along some time-like Killing vector satisfying Lie equation 

£^Um = -i^mUm, Wfn > 0. (2.39) 

To find such solution in the case of FST is a straightforward task related to the fact that 
the Poincare group is the symmetry group of Minkowski space-time. Therefore, FST 
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possesses the global time-like Killing vector £ M = (1,0,0,0). With this Killing vector 
Eq. (2.39) becomes 

-iu m u m , (2.40) 



dt 

from which it is clear that functions u m must be proportional to 



u m oc e 



—iu) m t 



(2.41) 



Let us conjecture that the full set of orthonormal mode functions have the form of 
plane waves 



.4 



k e i(k-x-£j k t) 



(2.42) 



(2^) 3/2 

where instead of the discrete indices m now we have continuous indices k, corresponding 
to the wave number of the plane wave. At the moment k is just a parameter of the 
mode function not related to the momentum. is the normalization constant which 
will be fixed later. We can easily check that these mode functions satisfy the orthonor- 
mality conditions in Eqs. (2.29). The frequency is defined using the Klein-Gordon 
equation (2.27) to be 



o; k = + 



+ m z 



(2.43) 



This equation is called relativistic dispersion relation. Note that we have chosen > 
in accordance with Eq. (2.39). 

With these mode functions the expansion of the field operator in Eq. (2.30) becomes 



4>(x) 



ake i(k.x-w k t) + a t e -i(k.x-wkt) 



dk 



(27T) 



3/2 " 



(2.44) 



And because k is the continuous index, the Kronecker delta in commutation relations of 
Eq. (2.34) must be changed into Dirac delta 



Ok, 



5 3 (k - k'l 



(2.45) 



with other commutators being zero. Using these conditions and commutation relations 
for the field operator in Eq. (2.33) we find 



(i,x),7r(i,x')] =i J ^ Wk (e tk -( x - x ')+. 



-ik(x-x') 



dk 

(2^ 



i5 3 (x - x) , (2.46) 
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which fixes the normalization constant 



(2.47) 



— 1/2 

A more physical motivation for the normalization constant being proportional to oc w k ' 
is that the expansion of the operator <f> in Eq. (2.44) with this normalization becomes 
relativistically invariant. 

To motivate the interpretation of Ok and as annihilation and creation operators 
let us exploit the space translational symmetry of the FST due to which the following 
relation must hold: 

$(t,x) = $(t,x + 6x), (2.48) 

where <5x is the infinitesimal translation vector. The process of spatial translation of the 
system may be described using a unitary transformation 



U = e iP1 , 



where 1 is the finite translation vector. Hence, we can write 



(2.49) 



>(t,x + l) = fr -1 0(i,x)*7. 



(2.50) 



For the infinitesimal translation 1 = 5x the exponent in Eq. (2.49) may be expanded 
to the first order as exp (j,P ■ 5xj = I + iP ■ Sx and the translational transformation 
Eq. (2.50) becomes 



(t, x + 6x) = \I-iP- <5xJ 4> (t, x) [I + iP ■ Sxj = (j> (t, x) + i ^ (t, x) , Pj Sx. (2.51) 
On the other hand for infinitesimal 5x we can also write 



■ (t, x + 6x) = 4> (t, x) + V x (t, x) • 5x. 



(2.52) 



Combining these two equations we find that the commutator of eft and P defines the 
gradient of the quantum field </>: 



(f>(t,x),P =-tV x ^(t,x) 



(2.53) 
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Let's concentrate now only on the positive frequency part of Eq. (2.44) 

0+ (x) = [ e 4(k - x - Wk ' ) a k dk, (2.54) 



where 

dk 

dk= — . 2.55 

(2vr) 3/2 

Inserting this expression into Eq. (2.53) we find 

ke ik - x a k dk = U+ (x) , p] . (2.56) 

Using commutation relations in Eq. (2.45) one can show that this equation is satisfied if 
the operator P is equal to 

P = j kdj^dk. (2.57) 
Acting with P on the state d^ |0) gives 

Pa{ |0, {«}) = kal |0, {«}) . (2.58) 

What does this relation mean? From Nother's theorem we know that translational 
invariance corresponds to the conservation of momentum. From which follows that the 
generator of the infinitesimal spatial translation is the operator for the total momentum, 
i.e. the operator P. From Eq. (2.58) it is clear that the state |k, {u}) = d^ |0, {u}) is the 
eigenstate of the total momentum P with the eigenvalue k. Remember, that until know 
k was just the index for the mode function. From the last relation k can be interpreted 
as the momentum and d^ acts as the momentum rising operator. 

If instead of using spatial translation symmetry we would have used time translational 
symmetry of the FST with the corresponding unitary operator 

f = e~ if[t , (2.59) 

we would have found that infinitesimal time translation gives 

H = j (w k o k d k + cfj dk, (2.60) 

where c is an arbitrary constant. The divergent constant term in the above equation in 
FST can be subtracted by appropriate procedures, but at the moment it must not concern 
us. The analogous arguments which relates P with the total momentum operator, leads 
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to H interpretation as the energy operator. Acting on the d^lO, {it}) state we would 
find that the operator a k raises the energy of the state by one unit and that can be 
interpreted as the total energy of that unit or quantum. 

2.2.2. Quantization in Curved Space-Time 

In the previous section we have described how fields are quantized in FST. This procedure 
is sufficient for the particle physics models, which studies only three fundamental forces 
of nature: electromagnetic, weak and strong. But to give a complete description of the 
Universe we need to study how all four fundamental forces, including gravity, shape and 
influence each other as well as the structure of the Universe. This requires a theory which 
puts all four forces on the same footing. Unfortunately such theory is still absent - the 
gravitational force resists the unification with the other three. In the presence of such 
resistance the only hope is to use a semiclassical description of the Nature, where we 
treat a classical gravitational background on which other quantized fields live. 

This approximation can be justified by noting that the Planck scale is the only scale of 
GR. If we consider small perturbations of the gravitational field and try to quantize them, 
then plays the role of the coupling constant. Hence, perturbation theory should be 
a good approximations for the energies much smaller than mp\. 

We gain confidence in this approach from the early development stages of the quan- 
tum electrodynamics (QED) theory, where the electromagnetic field was considered as 
a classical background on which fully quantized matter lives. And this method is fully 
consistent with a complete QED theory. 

2.2.2.1. From FST to CST 

The quantization of the field living in CST proceeds in the same line as the quantization 
in FST. First we write the action for the field. In CST the analog of Eq. (2.24) would be 



where T> g = det {gnu) is the determinant of the metric and is the covariant derivative. 
With the massive free scalar field Lagrangian, which is written in Eq. (2.26), this action 
becomes 




(2.61) 




(2.62) 
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Taking the variation with respect to the scalar field, dS/S(f> = 0, we arrive at the field 
equation (cf. Eq. (2.27)) 

(□ + to 2 ) <j) = 0, (2.63) 

where the □ operator is defined by 

□0 = gTV^Vv = -7^=8/1 (ypDggrdvA ■ (2.64) 



The scalar product of Eq. (2.28) for the Klein-Gordon equation in CST (Eq. (2.63)) must 
be generalized as well (e.g. Ref. [46]) 



(u m ,u n ) = i <5 tl M m dS' 1 , (2.65) 
Jt, 

where dS' 1 = n^dS and dS is the volume element in a given space-like hypersurface 
while is the orthogonal to this hypersurface time-like unit vector. It may be shown 
that the value of (u m , u n ) is independent on the choice of the space-like hypersurface S, 
i.e. (ii m , u n )-£ i = (u m , u n )-£ 2 . 

As in the FST, functions u m must satisfy the orthonormality conditions in Eq. (2.29) 
and then we can write a general solution of Eq. (2.63) as the superposition of a complete 
set of positive frequency {u m } and negative frequency {u^} solutions 



4> (x) = ^2 [ a mUm (x) + aL^m ( x ) ■ (2.66) 



The quantization of the field proceeds exactly as in the FST: change c-numbers <j) and 
7r into q-numbers <fr and n and impose canonical commutation relations of Eq. (2.33). 
Then operators a m and dm are interpreted as lowering and rising operators in the Fock 
space 

{u}) = Vn m + l|n m + l,M), (2.67) 

with the vacuum defined as 

a m |0,{n})=0. (2.68) 

As in the previous section {u} inside the ket reminds us that we are dealing with the 
vacuum defined by the complete set of orthonormal mode functions {u m }. This emphasis 
on the choice of mode functions becomes very important in CST as will be seen in a 
moment. 
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2.2.2.2. Bogolubov Transformations 

In section 2.2.1.1 it was shown that in FST a natural choice for the complete, orthonormal 
set of mode functions exist, which are plane waves of Eq. (2.42). And it was emphasized 
that this happens because the Poincare group is the symmetry group of the Minkowski 
space-time. Hence, using a global time-like Killing vector, d/dt, of this symmetry group 
we could pick-out positive frequency solutions Uk oc exp (— iu^t). And in all Lorentz 
frames, where t is the time coordinate, these mode functions define the same vacuum 
state. But in CST the Poincare group is no longer a symmetry group and in general there 
will be no global time-like Killing vectors in respect to which one could define positive 
frequency solutions. Therefore, the field expansion in mode functions {u m } in Eq. (2.66) 
is as good as in any other complete set of orthonormal functions: 

(2.69) 

m 

After quantization the vacuum state for this expansion is defined by 

S m |0,M)=0. (2.70) 

The definition of the vacuum state in Eq. (2.68) with mode functions {u m } at least 
formally differs from the definition with the mode functions {v m } in Eq. (2.70). Shortly 
it will be clear that this difference is not only formal but indeed both states |0, {u}) 
and |0, {v}) correspond to a different physical vacuum. Which means that there is no 
way to define uniquely a state without particles: what for one is a vacuum state, for the 
other this state contains particles. In such situation the notion of "the physical particle" 
becomes ambiguous. 

Since both sets of mode functions are complete orthonormal sets of solutions, each 
function in one set can be expanded in terms of the another set, i.e. 

<) . (2.71) 

m n 

These are the so called Bogolubov transformations, and matrices ot nm and (3 nm are called 
Bogolubov coefficients. It can be easily checked that these coefficients satisfy the relations 

/. \ ( a nl^ml - PmlPnl) = S nm, (2.72) 
I 

^ (UmlPnl - PmlOtnl) = 0. (2.73) 
I 



<t> (x) = ^ b m v m (x) + b ] m v* m (x) 
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Comparing Eqs. (2.66) and (2.69) and using Bogolubov transformations in Eq. (2.71) 
we find the relation between creation and annihilation operators of one set of mode 
functions and the other: 

Q>m = ^ ] (&nmbn + PnmPrij or = ^ ] (^nm^m ~ Pnm^rnj ' (2.74) 

n m 

Using these relations we can calculate the expectation value of the number operator 
defined by Nm = flm«m (cf. Eq. (2.37)). Acting with Nm on the vacuum defined by 
the mode functions {v m } we find 

(0, M| JV**> |0, {«}) = (0, {V}\ (<Ji + PnJn) (an'Jn* + Pl>J n ) |0, {v}) 

n,n' 

= (0, W| Yl PnmPn'Jnbl |0, {v}) (2.75) 
n,n' 

= ^ ] |/3mn| • 
n 

This shows that the vacuum defined by the complete set {v m } contains particles of the 
mode functions {u m }. 

The freedom of the choice of mode functions and the related ambiguity of the vacuum 
state constitutes the main problem of quantum field theory in curved space-time. One is 
naturally led to ask, which is "the physical vacuum" and what are the observables of such 
theory. In general, there is no way to pick out one particular set of mode functions. But 
in some space-times, which have a high degree of symmetry, this might be possible. As 
will be seen in the following subsection, this for example happens in a space-time with 
maximal spatial symmetry such as FRW and (quasi) de Sitter universes. The phenomena 
described in Eq. (2.75) are very important in inflationary particle creation. 

2.2.2.3. Quantization in Spatially Homogeneous and Isotropic Backgrounds 

In this section we describe the process of the scalar field particle creation in the expo- 
nentially expanding Universe. Let us consider spatially homogeneous and isotropic FRW 
metric of Eq. (1.2). But instead of using the cosmic time t we rewrite this metric in 
terms of the conformal time r defined as 




(2.76) 
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Then the line element with FRW metric becomes manifestly conformal to Minkowski 
space-time 

ds 2 = a 2 (r) (dr 2 - dx 2 ) . (2.77) 
In the FRW metric the action of the free massive scalar field written Eq. (2.62) becomes 

S=^J a 2 (V 2 - (V0) 2 - a 2 mV) dxdr, (2.78) 

where the prime denotes the derivative with respect to the conformal time, ' = 4= and 
V = di is the spatial gradient. As was already performed several times, taking the 
variation of this action gives the field equation 

(f>" + 2-<p' - V 2 </> + a 2 m 2 </> = 0. (2.79) 
a 

This equation is very similar to the Klein-Gordon equation in FST given in Eq. (2.27), 
except that it has a friction term 2a' /a ■ <p' . Let us transform this equation in such a 
way that it does become like the Klein-Gordon equation in FST. This can be achieved 
using the following mathematical trick, which brings any second order linear differential 
equation to it's normal form. If the equation is given as 



+P{x) — 
ax z ax 



2 +P(x)^ + Q(x)y = 0, (2.80) 



then the transformation 



u = ye y x p(*')^' (2.8i) 

brings it into the form of the harmonic oscillator 

For the equation (2.79) the analogous transformation would be 

x = 0e §/2^dr = a(r )^ (2.83) 

which transforms Eq. (2.79) into the form 

X" - V 2 X + (« 2 m 2 - X = 0. (2.84) 

This equation does look like the Klein-Gordon one in FST except the time varying mass. 
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The quantization of the scalar field x again proceeds as in the previous section: find the 
conjugate momentum of the field, which in conformal FRW space-time is tt = 5C/5x' = 
x' — ^rX) make changes of c-numbers into q-numbers, i.e. x ~ * X an d tt — )• 7r and impose 
canonical commutation relations (cf. Eq. (2.33)) 

[x ( T > x ) ) (r, x')] =iS(x- x) , [x (r, x) , x (r, x')] = [tt (t, x) , vr (r, x')] = 0. (2.85) 

The field operator x expanded into the creation and annihilation operators is written 
as (cf. Eqs. (2.66) or (2.69)) 



) = [ a ™Xm (») + a) m x* m (x) , (2.86) 



where x = (r, x) from the metric in Eq. (2.77). 

Mode functions Xm (x) must satisfy the orthonormality conditions in Eq. (2.29). The 
general scalar product of Eq. (2.65) in the FRW metric becomes 

(Xm,Xn) = i j {XnXm ~ XmXn) dx - ( 2 - 87 ) 

Because the (0, i) and (i, 0) components of the FRW metric in Eq. (2.77) are zero, the 
mode functions Xm (%) can be chosen in such a way that the temporal and spatial parts 
are separated 

Xm (r,x)^(2vr)- 3 / 2 Xk (r)e ik - x , (2.88) 

where k is now the continuous expansion coefficient and the factor (27r) ~ 3 / 2 is pulled out 
in order for the normalization of Xk ( T ) ( see Eq. (2.29)) to give the Wronskian of the 
form 

XkXk ~ Xtx'k = h (2-89) 

which is required from the orthonormality condition and is obtained using the scalar 
product in Eq. (2.87). 

With the ansatz in Eq. (2.88) the expansion of the operator x m Eq. (2.86) becomes 



X(x) = ] [akXk (r) e ik ' x + a\ X l (r) e~^ x j — w (2.90) 
Substituting x ( x ) Eq. (2.84) we find that functions Xk (t) must satisfy the equa- 
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tion of motion 



H 



2 Xk(r)=0 



(2.91) 



where we used a" /a' 



3 _ 



H + 2i7 2 . One can see that the equation of motion for Xk. is 



formally the same as of the harmonic oscillator with the time dependent frequency 



Although Eq. (2.91) constraints the time dependence of functions Xk (t) it does not 
determine the function uniquely. In fact, any function which satisfies Eq. (2.89) will be 
as good a choice as Xk ( r )- As was explained in the previous subsection, this fact deprive 
us of possibility to determine a vacuum state, defined as cik|0,{xk})) which would be 
seen as absent of particles by any observer. On the other hand, the quantum field theory 
in FST is a very successful theory although we do live in the expanding Universe, i.e. 
curved space-time. Therefore, we can expect that it is possible to pick out some special 
definition of the vacuum which would give correct predictions for laboratory experiments. 
The main reason why flat space-time QFT is so successful from this point of view is that 
it describes phenomena which take place in a very weak gravitational field, or in other 
words, very weakly curved space-time, which may be neglected. 

This can be easily seen from Eq. (2.91). If we neglect the effect of gravity, which 
corresponds to taking a = 1 and therefore H = 0, the frequency term in Eq. (2.92) 
becomes constant, u 2 . = m 2 + k 2 = const, the same as in FST in Eq. (2. 43). 3 With 
a constant frequency, functions Xk ( T ) ( or functions u m (t) in Eq. (2.41)) have the time 
dependent part exp(— ioj^r). The vacuum defined in this way will be the same for all 
inertial observers at all times. But the frequency term in the expanding Universe in 
Eq. (2.91) is time dependent. Hence, the vacuum defined at time Tj will contain particles 
as seen by the observer at some later time. This is the main reason why particles get 
produced during inflation, but let us postpone this discussion until a bit later. At the 
moment the important thing is the choice of initial conditions which would fix the form 
of mode functions and therefore the initial vacuum state. 



3 This motivates us to interpret the expansion coefficient k as the comoving momentum of the particle 
(cf. section 2.2.1.1) and k/a as the physical momentum. 




2 



(2.92) 



46 



2. The Origin of the Primordial Curvature Perturbation 



2.2.2.4. The Vacuum State in FRW Background 

Although, as was mentioned earlier, in general, the particle concept in CST is ambiguous, 
in some special cases it is possible to define an approximate particle concept which would 
be as close as possible to the one known from QFT in FST. This is the case, for example, 
in space-times described by the FRW metric or anisotropic Bianchi universes. In the case 
of the present interest we may look for the solution of Eq. (2.91) with the ansatz 

X t W~e-^^^^ (2.93) 
V 2Wk (r) 



The factor l/v2 is chosen so that Xk ( T ) would satisfy Eq. (2.89) and the function (r) 
satisfies 

which can be found by substituting the ansatz in Eq. (2.93) into Eq. (2.91). If the time 
variation of u>k (r) is very slow, it is said that it satisfies the adiabatic condition and the 
vacuum defined when this condition is valid is called the adiabatic vacuum. By "slow" 
we mean that oj^ (r) and all its derivatives change substantially, Aoj^/lo^ ~ O (1), only 
during the time interval T 3> lo^ 1 (Ref. [48]). In the adiabatic case, derivative terms in 
Eq. (2.94) will be small and this equation can be solved using the recursive method. For 
example, to the zeroth order we can take 

wi 0) (T)=u k (T). (2.95) 



Note that for the constant frequency, = const, the mode functions in CST (Eq. (2i 
with Xk (t) given by Eq. (2.93) reduce to the mode functions in FST (cf. Eqs. (2.42) and 
(2.47)). 

In Eq. (2.91) the adiabatic vacuum can be defined for light particles (m/H <C 1) whose 
Compton wavelength is much smaller than the curvature scale H , or in other words 
whose physical momentum is much grater than the Hubble expansion rate, k/a S> H. 
We may say that such particles do not "feel" the gravitational field. Therefore, by sub- 
stituting (t) into Eq. (2.93) and taking that 

w k « k, (2.96) 



4 For the quasi de Sitter expansion H/H 2 « 1 (cf. Eq. (1.42)) and for the FRW Universe H/H 2 ~ O (1), 
so that, when k/a^$> H, terms of order one or less are subdominant in Eq. (2.92). 
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we find the initial condition for the mode function Xk ( r 

1 



X k (2.97) 



which is the same as that of the massless field in FST (cf. Eqs. (2.42) and (2.47)). This 
vacuum state is often called the Bunch-Davies vacuum. 

The result of Eq. (2.97) was obtained by assuming that the space-time can be consid- 
ered flat at the zero order approximation for subhorizon modes. This must be always valid 
for the Einstein gravity in accordance to the equivalence principle. But in many models 
the inflationary energy scale is just couple of orders of magnitude below the Planck scale. 
At such energy scales it might be that Einstein's theory of gravity is not precise enough 
to describe Nature. In this case, one may expect that the equivalence principle does not 
hold anymore. But this failure might be only at the level of (H/mpy) 2 < 10~ 10 , where 
H is the inflationary Hubble parameter [5]. 

2.2.2.5. The Field Perturbation in the Inflationary Universe 

It was already mentioned in section 1.5 that inflation provides a natural mechanism to 
explain the origin of the curvature perturbation in the early Universe. Upon entering the 
horizon this perturbation sets the initial conditions for the tiny density inhomogeneities 
which seeded the subsequent growth of large scale structure such as galaxies and galaxy 
clusters. In this subsection we describe how field perturbations are generated in the 
inflationary Universe and in the next section how they are transformed into the curvature 
perturbation. 

The generation of the field perturbation can be computed from Eq. (2.91) with appro- 
priate initial conditions and assumptions relevant for the inflationary expansion. Dur- 
ing inflation the Universe undergoes quasi de Sitter expansion for which the condition 
H /H 2 <C 1 is satisfied (see Eq. (1.42)). But for our purpose in this section and for 
ater discussions in Chapter 3 it is enough to take the approximation of a quasi de Sitter 
Universe. Hence, we will set H = which is equivalent to considering exact de Sitter 
space-time. 5 

Another assumption we make is that initially the state corresponds to its vacuum, 
i.e. the average occupation number of particles with momentum k is much less than 
1, nfc <C 1. This assumption is easily justified if enough amount of inflation occurred 
before the horizon exit of the scales of interest (see e.g. Ref. [5]). In this case the 



5 More precisely only a part of de Sitter space-time is considered because inflation lasts only for a finite 
time. 
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initial conditions for the mode functions is determined by the Bunch-Davies vacuum in 
Eq. (2.97). 

Perturbations of scalar fields with the mass comparable to the Hubble parameter do 
not grow significantly as will be seen in Eq. (2.105). While in section 2.2.2.6 it will 
be shown that perturbations of heavy fields, with the mass m > |-£f, do not become 
classical after horizon crossing. Therefore, only light scalar fields are considered for the 
generation of the curvature perturbation. Assuming de Sitter expansion the general 
solution of Eq. (2.91) becomes 

X*(T) = ^e^^HP(-kT), (2.98) 

where the initial state was matched to the Bunch-Davies vacuum in Eq. (2.97), and 
{—kr) denotes the Hankel function and we used in the de Sitter space-time r = 
— (aH)~ l . The order of is defined as 



9 /m\ 2 
4 ~ \H, 



(2.99) 



Well after horizon exit, when \kr\ <C 1, this solution approaches to 



Xk(r)^ 6 = P ^ . (2.100) 



2k V 71 " 

For the light field m < the parameter v is real and this solution is not oscillatory, 
therefore interpretation of corresponding states as physical particle states in the Fock 
space is problematic. The reason for this, as can be seen from Eq. (2.92), is that the 
dispersion relation for a light field, with m <C H, becomes imaginary and the mode 
function does not oscillate. 

However, the amplitude of quantum fluctuation in the state \ip) is always well defined 
and we can calculate the expectation value for the vacuum state |0) as 



1 f°° , s , , ..osmkLdk f°° ^ „ , sinfcL d/s 

x(r,x)x(r,y)|0) = _/ o *.| Xk(T) |»__ = y r x (k)— T , 



(2.101) 

where L = |x — y| and V x (k) = (/c 3 /27r 2 ) |xk| 2 is t ne power spectrum. 

The power spectrum of the superhorizon massive scalar field perturbations can be 
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easily calculated using Eq. (2.100) 

p * = -ir-(to) (^h) • (2 - 102) 

However, this expression is derived for the comoving field \- Going back to the physical 
field <j) = \ja (see Eq. (2.83)) and considering the massless limit v = 3/2, this expression 
reduces to 

H\ 2 



^=(-1 . (2.103) 



Or more generally, for a light field, m < we can express Eq. (2.99) as 



3 m 2 



2 3# 2 ' 



(2.104) 



and the power spectrum becomes 



fH\ 2 f k 

(=-) (~) ■ (2-105) 



V2vr7 \2aH 

Although Eq. (2.103) was calculated for the exact de Sitter expansion, one can easily 
include a very slow variation of the Hubble parameter, H ^ 0. Due to this small variation, 
the horizon size changes very slightly during inflation and therefore each mode k exits 
a horizon of slightly different size. But from Eq. (2.103) it is clear that the amplitude 
of the field perturbation is proportional to the horizon size. Therefore, different modes 
will have slightly different amplitudes. This mild dependence of the power spectrum on 
k may be accounted for in Eq. (2.103) by writing 

where the Hubble parameter in this expression has to be evaluated at the horizon 
exit for each mode, i.e. when aHk = k, and is slowly varying with k. 

It is important to note that because we have assumed the initial state for each mode to 
start in the Bunch-Davies vacuum and considered a free field, perturbations of the field 
are Gaussian. This is the result of the equivalence principle valid for Einstein's gravity. 
But if at inflationary energies the equivalence principle does not hold, for example due 
to modified gravity, the perturbations of the field may be significantly non-Gaussian. 
In addition, because we were concerned in this section about scalar fields which are 
rotationally invariant, the perturbations of the field are statistically isotropic. This will 
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not be the case in Chapter 3 where we discuss perturbations of vector fields. 
2.2.2.6. Quantum to Classical Transition 

As it is clear from the discussion so far, the origin of the field perturbation is quan- 
tum mechanical. But as was claimed in section 1.5.1 the greatest success of inflationary 
paradigm is that it can explain how these quantum mechanical perturbations give rise 
to the initial density perturbations in the Universe which are observed as CMB tem- 
perature anisotropics and which seed the formation of galaxies. But CMB temperature 
anisotropics and galaxies are not quantum but classical objects. Hence, there must be 
some transition period where quantum mechanical perturbations are transformed into 
classical ones. This process is analogous to the decoherence in quantum mechanics. The 
result of it is that the quantum mechanical superposition principle is violated and the 
wavefunction collapses to a particular state obeying the classical evolution. The co- 
herence between different states is lost, since after the collapse only one state can be 
observed, although quantum mechanically all states should be allowed. In usual applica- 
tions of quantum mechanics this happens due to the wavefunction interaction with the 
degrees of freedom of the environment. But in cosmological context the transition from 
quantum-to-classical does not require environment, therefore in Ref. [49] it was named 
"decoherence without decoherence". In the exposition of quantum-to-classical transition 
below we will follow Refs. [49, 50, 51, 52]. 

What does it mean, that perturbations in the Universe are classical? As was discussed 
in section 2.1 these perturbations can be described by random fields (3. The classicality 
of perturbations means that /3 is described by classical stochastic variables. But to make 
the discussion easier, instead of treating some general variable (3 let us specialize to a field 
perturbation x ( x ) defined in Eq. (2.83). Then x ( x ) will be classical if it is described as 
classical stochastic variable with the probability distribution function p(\x\ > KD) where 
7r is the canonical conjugate of X- 

In quantum mechanics a classical limit is achieved when the state collapses to the 
definite numerical value. But in the cosmological context we cannot assign the definite 
numerical value to a collapsed state. So we say that the quantum state becomes clas- 
sical if the field modes become equivalent to the classical stochastic functions with the 
probability distribution pdxU 71 "!) [49]. This can be written as 




(2.107) 
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where from Eq. (2.86) 

X(x)= I X m (x)dm = I a m Xm (x) + aJ„Xm ( x ) dm, (2.108) 



m is a continuous index and X m \ = Re (X m ), X m 2 = Im (X m ). Note that in this equation 
operators are denoted by hats and classical fields without hats. 

Of course Eq. (2.107) is not valid in general. But this equality is valid when quantum 
fields can be treated as classical, i.e. when conjugate variables commute. To show this 
let us use a quantum field \- For the moment restoring physical units, the non-zero 
commutation relation for this operator and its conjugate pair in Eq. (2.85) becomes 

[ X (r, x) , vr (r, x')] = iHS (x - x') . (2.109) 

The classical limit of a quantum description must be achieved when the Planck constant 
becomes negligibly small, H — > 0. In this limit the commutator in Eq. (2.109) becomes 
zero. If the operator x is expanded into the complete set of orthonormal mode functions 
as in Eq. (2.86), then in the limit H — >• the orthonormality condition in Eq. (2.29) for 
the mode functions {xm} becomes (Xm,Xm) — > 0. Using Eq. (2.28) we find that this 
condition results in the mode function Xm and its complex conjugate Xm being different 
only by the time independent phase factor 

X*m = CmXm- (2.110) 

But the phase of Xm is completely arbitrary. Therefore we are free to choose it in a way 
that makes Xm real. And because c m is time independent, Xm is rea l at all times and 
X m in Eq. (2.108) can be rewritten as 

X m (x) = Xm{x) a m + cJ m . (2.111) 

One can easily check that with this expression the commutator in Eq. (2.109) is zero. 
In Refs. [49, 50] it was calculated explicitly that an operator of the form in Eq. (2.111) 
satisfies the equivalence equation (2.107) for the quantum field and stochastic classical 
field. Which shows as well that the classical stochastic field can be expressed as 

X m (x) = Xm (x) e m , (2.112) 

where e m are time independent, complex, stochastic c-number functions with zero average 
and unit dispersion: (e m ) = and (e m , e* ) = 5 (m — re). And e m obeys the same statistics 
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as X m (t ,x). 

Note that the time dependent part Xm( T , x ) can be factored out from both: the 
quantum field X m in Eq. (2.111) and from the classical field X m in Eq. (2.112). As 
the result the evolution of a given mode function is completely deterministic after the 
realization of some stochastic amplitude have occurred. In other words, if we measure 
the amplitude of the field perturbation some time after the horizon exit, it will continue 
to have a definite value. 

To show how a quantum field becomes of the form in Eq. (2.111) in the accelerating 
Universe, let us consider such a field in the de Sitter background. In subsection 2.2.2.5 it 
was shown that if we choose mode functions Xm (%) to be Fourier modes (see Eq. (2.88)) 
on the superhorizon scales they will have the solution 



where we have used Eq. (2.100) and v is defined in Eq. (2.99). If v is not imaginary, 
corresponding to m 2 < 9H 2 /A, this mode function does not have a time dependent 
phase and can be made real by choosing a time independent phase rotation Ck (defined 
in Eq. (2.110)). Note that this transformation makes Xm re al at all times on the super- 
horizon scales. Therefore, quantum mechanical operator X m satisfies Eq. (2.111) and 
consequently is equivalent to the classical stochastic field X m . 

The process of quantum-to-classical transition described above suffers from the usual 
interpretational problem of measurement in quantum mechanics. The first question is 
why did Nature choose this particular value for the realization of the field amplitude 
when other infinite possibilities were available? Another is a cosmological variant of the 
Schrodinger's Cat problem related to the question of when the state collapsed into its 
observed value. According to the usual Copenhagen interpretation, this happens at the 
time of measurement. But does that mean that CMB perturbation pattern did not exist 
before we have measured it for the first time? 

2.3. The Primordial Curvature Perturbation 

Quantum field perturbations described in the last section on superhorizon scales give rise 
to the classical cosmological perturbations. These perturbations are most conveniently 
described by the intrinsic spatial curvature £, which is commonly called the curvature 
perturbation. It is defined on the hypersurfaces of constant energy density. We will post- 
pone the discussion how quantum field perturbations may generate £ until section 2.4. In 




(2.113) 
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this section we discuss the properties of £ and the formalism which relates perturbations 
of quantum fields with the curvature perturbation. 

2.3.1. Gauge Freedom in General Relativity 

When discussing the FRW Universe in section 1.1 we implicitly chose the coordinate 
system in which the metric in Eq. (1.2) attained it's elegant form. Although the physical 
results in GR should not depend on the coordinate system, the homogeneity and isotropy 
of the Universe singles out a preferred reference frame in which equations reduces to their 
simplest form. However, in space-times without symmetries such preferred coordinate 
system does not exist. Therefore, the choice of coordinates is purely arbitrary and may 
be selected depending on the problem at hand. 

Fixing the coordinate system in GR specifies how space-time is threaded by the lines 
of constant spatial coordinate x {threading) and how it is divided into the hypersurfaces 
of constant coordinate time t (slicing). And because the coordinate system is arbitrary, 
so are the threading and slicing. 

The arbitrariness of the coordinate system becomes especially problematic when de- 
scribing tiny departures from homogeneity of the actual Universe. Since the real space- 
time with these departures does not posses any symmetry, it would be impossible to 
solve the exact relativistic evolution equations because GR is a non-linear theory. But 
luckily on large enough scales these departures from ideal homogeneity are very small, 
only of the order 10 -5 (see the discussion on the cosmological principle in section 1.1). 
Therefore, a very good approximation to the actual Universe is to treat these inhomo- 
geneities as tiny perturbations of the otherwise homogeneous and isotropic background. 
And using perturbation theory, non-linear equations of GR may be linearized. 

However, separating physical quantities into background value and perturbations are 
not so trivial. Due to the freedom for the choice of the coordinate system this separation is 
not unique. Such complication is due to the fact that by separation we mean that for each 
space-time point in the background or reference manifold we associate a perturbation, 
corresponding to the actual or physical space-time. But because these are two different 
manifolds with different curvature we must specify how the mapping from one manifold 
to the other is performed. This may be done by choosing the specific threading and 
slicing of the space-time. Or in perturbation theory it is called by fixing the gauge. 

By changing the gauge we must redefine what is the background value and what is 
the perturbation. Let us consider an infinitesimal gauge transformation in which new 
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coordinates xf 3, are related to the old ones by 

= x » + fix* ( x ) , (2.114) 

and see how the perturbation of some scalar quantity 5f (x) changes by this transforma- 
tion, where the perturbation is defined as the difference between the actual value / (x) 
and the background value fo (x), Sf (x) = / (x) — fo (x). In the new gauge this perturba- 
tion will be Sf (x) = f (x) — fo (x). But because the physical value of the quantity does 
not change under the gauge transformation, only the separation into the background and 
perturbation does, we can write / (x) = f (x). Further, because / is a scalar, it should 
be invariant under the coordinate change, thus f (x) = f (x). On the other hand, we 
are keeping fixed the point in the background manifold and investigate the change in the 
mapping to the perturbed manifold, therefore fo (x) = fo (x). However, although it is 
the same point on the background manifold, in a new coordinate system it will have a 
different value. With the infinitesimal transformation in Eq. (2.114) this may be written 
as fo (x) = fo (x) + (dfo (x) jdx) ■ 5x. Putting all this discussion together, we may write 

Sf(x) = /(x)-/o(x) 

= f{x)-h{x) (2.115) 

= 5f{x) + fo{x)-^±5x»-fo{x) 

In the later discussion, of special importance will be the transformation of the scalar 
quantity when only the slicing is changed, keeping the threading constant. This corre- 
sponds to a time shift, for which Eq. (2.114) reduces to 

t = t + 5t(x). (2.116) 

And Eq. (2.115) becomes 

6f-6f = -f 6t. (2.117) 

2.3.2. Smoothing and The Separate Universe Assumption 

Calculations of the curvature perturbation in this thesis are performed using the separate 
Universe assumption. And the central concept for this assumption is that of smoothing. 

Let us assume that we are interested in the perturbation of the energy density bp. 
Then the smoothed value 5p(t,x,L) will correspond to the value of the energy density 
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a b 

Figure 2.1.: The schematic illustration of the separate Universe assumption. According 
to this assumption all threadings coincide with the comoving one. Hence, 
curves (a) and (b) represent comoving trajectories of any of these threadings 
(adapted from Ref. [53]) 

at the space-time point (t, x) after averaging over the sphere of comoving size L. If we 
expand 5p(t,x,L) in a Fourier series, smoothing would correspond to dropping off all 
modes bigger than k ~ L . The dynamics of this smoothed quantity is assumed to 
be determined by the averaged Einstein equations. However, GR is a non-linear theory 
which makes the issue of smoothing non trivial. For example, it is not clear how small 
scale fluctuations with k > can influence the evolution of the quantity smoothed 
on scales L . At present there is no satisfactory conclusion to this issue, but we will 
assume that such length scale does exist above which the smoothed Universe is a good 
approximation to the actual one. 

The separate Universe approach assumes that each region, smoothed on scales larger 
than the horizon size, locally evolves as a separate unperturbed Universe. The basic 
idea is presented in Figure 2.1. In this figure L corresponds to the smoothing length 
scale which is larger than the horizon size H -1 but smaller than the largest box size 
i?boxj within which we perform our calculations (see Eq. 2.10). Ideally i?box — * °o, but, 
as was discussed in section 2.1.1, one should keep a box size finite in order to avoid 
unknown physics and keep calculations under control. The lines (a) and (b) represent 
two comoving worldlines for two different space points. 

Another assumption made in the separate Universe approach is that all length scales 
introduced by the energy momentum tensor are much smaller than the smoothing length 
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scale Then is the only relevant superhorizon length scale and all spatial gradi- 
ents of order k/a are negligible. When this assumption is satisfied, the locally measurable 
parts of the metric should reduce to those of the FRW [54]. In other words, every comov- 
ing location smoothed on superhorizon distances, evolves as the unperturbed Universe 
with the FRW metric of Eqs. (1.1) or (1.2). 

2.3.3. Conservation of the Curvature Perturbation 

Every smooth space-time metric can be decomposed into 3+1 components as [55]: 

ds 2 = N 2 dt 2 - jij (dx i + ftdt) (dx j + (3 j dt) , (2.118) 

where M is the lapse function, /3 l is the shift vector and — 7^ is the spatial three metric 
tensor. 

With this decomposition one can define a unit time-like vector, n^, normal to the 
hypersurface of constant coordinate time t. The components of this vector are 

n M = (AA,0); „M = r__LjH. ( 2.n9) 

Then the volume expansion rate of the hypersurface along some integral curve 7 (r) of 
will be given as 

■& = V m tj/\ (2.120) 

where V M is the covariant derivative and r is the proper time, which can be found from 
Eq. (2.118), dr = Mdt. Along each of these integral curves we may define the number 
of e-folds of expansion 

iV(ti,t 2 ;x) = I [ Mr = I $Afdt, (2.121) 

6 Jy(r) 6 Jti 

where the vector x is chosen to be the comoving spatial coordinate. 
The spatial metric of Eq. (2.118) can be further decomposed as 

7 y =a 2 (i,x)7 ij . (2.122) 

With the requirement det (7y) = 1, a (t, x) becomes a local scale factor. Note, that t in 
this equation is not necessarily a proper time, it is just the coordinate time labeling the 
slices. Since we are interested in the non-homogeneity of the scale factor, we may further 
decompose a(t, x) into some global scale factor a(t), which is independent of position, 
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and the local deviation ip (t, x) 

a(t,x) = a(i)e^*' x) . (2.123) 

This decomposition into the global quantity and its perturbation is completely arbitrary. 
Hence, we may choose a (t) in such a way that ip (t, x) vanishes somewhere inside the 
observable Universe. Then if) (t,x) becomes small everywhere inside this Universe [54]. 
A similar decomposition may be done for the 7^ part of the metric 

7y=(le ft ) , (2.124) 

where I is the unit matrix and hij must be a traceless matrix due to the requirement 
det (7) = 1. It can be shown that h%j corresponds to the primordial tensor perturbation, 
i.e. gravitational waves. But, for the time being, we assume that GWs are negligible so 
that we can set h = 0. 

According to the separate Universe assumption, if the metric is smoothed on super- 
horizon scales, at each space-time point we should be able to find such coordinates which 
reduce the metric into the form of FRW: 

ds 2 = dt 2 - a 2 (t) flydxW. (2.125) 

This metric is chosen to be flat in agreement with observations, but as noted in Ref. [54] 
a small homogeneous curvature should not make much difference. 

In accord with this assumption and with the appropriate coordinate choice, the metric 
in Eq. (2.118) should reduce to the form 

ds 2 = M 2 dt 2 - a 2 (t) e 2 ^ >3 %dxM^'. (2.126) 

The separate Universe assumption does not pose any constraints on J\f and ip since they 
are not locally observable quantities. And in view of this assumption we have neglected 
all terms of order O (k/aH), which on superhorizon scales approach zero, k/aH — > 0. In 
particular, in Ref. [54] it was shown that 

Pi = O (k/aH) and % = O \(k/aH) 2 ] . (2.127) 

In the following discussions we will keep in mind that the separate Universe assumption 
is valid up to this order, but will omit terms O (k/aH) from equations. Note, however, 
that the smallness of just corresponds to our choice of the coordinate system. The 
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generalization of the formalism to a threading with non negligible is straightforward 
[54]. 

With the line element of Eq. (2.126) the local expansion rate defined in Eq. (2.120), 
takes the form 

•£>(() 



"^UM+^J- (2 - 128) 

For later convenience we define the local Hubble parameter H (t, x) = ^fi: 

^=K^! + ^ x >)' (2 - 129) 

In what follows, an important quantity is the number of e-folds of the local expansion, 
which is defined in Eq. (2.121). With the line element in Eq. (2.126) it becomes 

iV(t 1 ,t 2 ;x)= f H(t,x)dr= (^ +i,(t,x))dt. (2.130) 

J -fir) Jtl V°W / 

According to the separate Universe assumption each space point evolves as the unper- 
turbed Universe with the locally defined expansion rate in Eq. (2.128) (or equivalently 
local Hubble parameter in Eq. (2.129)). Therefore, at each point we can write the energy- 
momentum conservation law, V U T^ U = 0, from which it follows 

= -3H (t, x) [p (t, X )+p (t, x)] . (2.131) 

It has the same form in the FRW Universe (c.f. Eq. (1.4)). 

This equation is valid independently of the slicing. But let us specialize further to the 
slicing on which energy density is uniform, i.e. independent on space coordinate at each 
given time. Such slicing is called comoving or uniform density slicing and the value of 
ip on this slicing is usually denoted by £. It determines the perturbation in the intrinsic 
curvature of the slices. Then, Eq. (2.131) can be rewritten as 



P(t) 



a(t) 



[p(t)+p(t,x)], (2.132) 



where we have used Eq. (2.129) as well. 

Now let us limit ourselves to the case where pressure is adiabatic, which is equivalent 
to saying that pressure is a unique function of the energy density, i.e. p = p (p). In this 
case, because p is independent of position, the pressure must be independent of position 
as well, p = p(t) only. Therefore, the same must be true for £, i.e. £ = £ (t) only. 
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On the other hand, the decomposition of the spatial part of the metric in Eq. (2.123) 
into the background value a (t) and deviation from that value tp (t, x) was purely arbitrary. 
So we may choose the normalization of a (t, x) such that o (t) corresponds to the scale 
factor at our location (or any other location). In other words, we choose a (t) in such a 
way that £ vanishes at our location at all times. Hence, at this location 

C = 0. (2.133) 

But because £ is independent on position (when p = p(p)) it must vanish everywhere. 

In this way we have found a very important quantity, the curvature perturbation £, 
which determines the intrinsic curvature of constant time spatial hypersurfaces. As was 
shown above, on superhorizon scales £ stays constant whenever pressure is the unique 
function of the energy density. In the history of the Universe this happens when the latter 
is dominated by radiation or matter. More generally, the pressure of the multicomponent 
fluid is adiabatic if each component of the fluid satisfies the relation p a = p a (p), where 
p a and p are the energy densities of each component and of the total fluid respectively. 
Thus, around the matter-radiation equality era, £ is constant too if perturbations are 
adiabatic. By adiabatic perturbations we mean that on uniform total energy density 
slices, perturbations of each component are independent of position. A more rigorous 
proof of the constancy of the curvature perturbation £ can be found for example in 
Refs. [54, 53]. In Ref. [53], the constancy of £ was proved using perturbation theory, 
without the assumption of separate universes. As shown in these works, the change in 
the curvature perturbation to the first order is proportional to the non-adiabatic part of 
the pressure 5p na d as 

C = —Spnad, (2.134) 

P+P 

where <5p na d is defined as the pressure perturbation on the uniform density slicing. 



2.3.4. The SN Formalism 

On superhorizon scales all threadings are equivalent to the comoving threading up to 
the order O (k/aH) 2 [54]. Hence, only the slicing is arbitrary. In this section we show 
how using this fact and the separate Universe assumption one can relate the curvature 
perturbation £ with the energy density perturbation Sp without invoking cosmological 
perturbation theory. 

Let us consider one of the comoving threads drawn in Figure 2.1. Going from the 
coordinate time t\ to ti along this thread we can calculate the change in the value of ip 
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for some slicing using Eq. (2.130) 

^(t2,x)-^(ti,x) = JV(ti,t 2 ;x)-JVo(tl,fe), (2-135) 

where Nq (ti,t 2 ) = m [a (£2) / a (ti)] is the number of e-folds for the background expansion. 
For the flat slicing we have ^g at = 0, and thus the number of e-folds for the local expansion 
coincides with the background one, 2Vfl a t (ti,t 2 ;x) = iVo (ti,t 2 ). 

Let us apply now Eq. (2.135) for two different choices %[)\ and tp-Q corresponding to two 
different slicings, which coincide at time t\. Then at time t 2 the difference between ip A 
and ip-Q will be 

rP A (t 2 , x) - <0b (h, x) = N A (tx , t 2 ; x) - N B (h , t 2 ; x) . (2 .136) 

Let us further specify slicings A and B in the following way. The slicing B will be the 
flat slicing, giving = Nq. And let the slicing A be such that at time t\ it coincides 
with the flat slicing, ipA (t\, x) = 0, while at time t 2 it coincides with the uniform density 
slicing, V'A (^2; x ) = C (*2) x )- Then Eq. (2.136) takes the form 

((t 2 ,*) = N A (t 2 ; X )-N (t 2 ). (2.137) 

Due to our choice of A such that ip A (ti,x) = 0, the number of e-folds of the local 
expansion on this slicing becomes N A (t 2 ,x) = In [a (t 2 ,x) /a (ti)]. This means that 
Eq. (2.137) is independent on the initial time and this is why we omitted the notation 
of t\. In other words, the calculation of £ is independent on the initial epoch, because 
when going from one flat slice to the other the expansion is uniform. 

From this equation it is clear that the curvature perturbation £ (t, x) specifies the 
perturbation in the number of e-folds of the local expansion starting from any flat slice 
and ending on the uniform density slice at time t: 

C(t,x) = SN(t,x). (2.138) 

Until now our discussion didn't require that perturbations should be small. Hence, 
they are valid to any order in the perturbation expansion. To relate SN (t, x) with the 
field perturbation (discussed in section 2.2.2) we will need to specialize further in small 
perturbations. 

Let us assume that the local expansion of the Universe is determined solely by the 
value of a classical scalar field, (t, x) = N {<p (t, x)). By this assumption we neglect the 
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contribution, for example, by the kinetic term of the field, <p (t, x). This is valid in most 
cosmologically interesting cases, for example during inflation with the slowly varying 
field. Then Eq. (2.137) can be written as 

C(t,x) = N(<f>(t,x))-N(<f>{t)). (2.139) 

Taking the field perturbation to be small 8(p (t, x) <C (p (t), where (p (t, x) = (p (t) + S(p (t, x) , 
this equation becomes 

C (t, x) = N^S<t> + ^N H (5(f)) 2 + ..., (2.140) 

where = dN ((p (t)) / d(f> and N H = d 2 N ((p (t)) / d(p 2 . Note that derivatives are 
taken of the unperturbed value of N, and the field perturbation is evaluated on the 
initial flat slice. This equation can be easily generalized to the many field case, when 
N (t,x) = N (tpi (t, x) , (p2 (t, x) ,...), in which case Eq. (2. 140) becomes 

C (t, x) = Ni8<i>i + \Y1 N uHiHj + ■■■■ (2.141) 
/ ij 

In the rest of this thesis it is sufficient to consider the curvature perturbation only to 
the second order in the field perturbations, i.e we will drop out terms denoted by the 
ellipsis. 

2.3.5. The Power Spectrum and Non-Gaussianity of C 

To calculate the power spectrum and bispectrum of the curvature perturbation, £ must 
be transformed to the Fourier space by Eq. (2.2). Then Eq. (2.141) becomes 

Ck = N+Sfa + ijV^ (5(P k ) 2 . (2.142) 

Note, that in this expression ^ is dependent only on the modulus of k = |k|. This is 
because perturbations of the scalar field are rotationally invariant. We will drop this 
assumption in Chapter 3 when discussing perturbations of vector fields. 

The two point correlation function for the curvature perturbation in Eq. (2.142) is 

(C* (t) , Ck' (*)> = N 2 (5<p k , Sfa) + \n 2 h ((5(P k f (5(P k ,f) . (2.143) 

Because we have assumed that field perturbations are Gaussian, the first term of this 
equation is the Gaussian contribution. The second term gives a non-Gaussian contribu- 



62 



2. The Origin of the Primordial Curvature Perturbation 



tion and according to observations this contribution must be subdominant. 

Taking only the dominant part in Eq. (2.143) we find that the power spectrum of the 
curvature perturbation is related to the power spectrum of the field perturbation by 



where the power spectrum of the field perturbations V§ is the one in Eq. (2.103) for de 
Sitter inflation or in Eq. (2.106) for the slow-roll inflation. 

As was discussed in section 2.1.1 if perturbations are Gaussian the two point correlator 
is the only non-zero correlator. The non-Gaussianity manifest itself in the non-vanishing 
higher order correlators. For the aim of the present thesis it is enough to consider only 
the three point correlator, although in some models it might be that, for example the four 
point correlator is even larger than the three point (see e.g. Ref. [56]). The bispectrum 
(defined in Eq. (2.13)) of the curvature perturbation usually is parametrized by the non- 
linearity (or non-Gaussianity) parameter /nl defined in Eq. (2.20). 

If the field perturbation 5(f) is Gaussian, then /nl becomes practically independent of 
k. In Ref. [57] it was calculated that if the first term of Eq. (2.142) is dominant the /nl 
parameter becomes 



2.3.6. Density Perturbations 

In the previous subsection we have shown how to calculate the curvature perturbation 
C, which is conserved on superhorizon scales whenever the pressure of the cosmic fluid 
is adiabatic. However, our primary interest is in small perturbations of the energy den- 
sity, which upon horizon entry form the seeds for the subsequent structure formation in 
the Universe. To make a connection between the curvature perturbation £, which we 
calculated so far, and inhomogeneities of the energy density, we will use a limit of small 
perturbations up to the first order. 

As was discussed earlier, on superhorizon scales the threading is defined uniquely, 
changing the slicing corresponds only to a shift in the coordinate time. So let us consider 
a change from the uniform density slicing to some generic one. At any given position this 
will correspond to a time change 5t (t, x), so that on a new slicing t = t + St (t, x). Then, 
the local scale factor on the new slicing can be found using Eq. (2.117) and considering 
that the background value is the same for both slicings 



V c = Np>4> 



(2.144) 




(2.145) 




(2.146) 
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Separating local scale factors into the background value and the perturbation as was 
done in Eq. (2.123) and considering that the perturbation is small, we find to first order 

rp = (-H6t. (2.147) 

A similar reasoning applies to the energy density giving 

6p4,(t,yL) = -p(t)6t(t,x), (2.148) 

where Sp = on the uniform density slicing. For the time being we use the index '■0' to 
remind ourselves that density perturbation is defined on an arbitrary slicing. While 5p 
without this index will correspond to the density perturbation in a flat slicing. 
Combining the last two equations we arrive at 

p{t) 3 p + p 

where in the last equation the continuity equation for the FRW Universe (Eq. (1.4)) was 
applied. Thus we derived the equation for the transformation going from the uniform 
density slicing to arbitrary slicing. An important choice of the latter is the flat slicing, 
ip = 0, for which we get 

C = ~H 5 4 = \^. (2.150) 
p 3 p + p 

To derive this equation the uniform density slicing was defined with respect to the 
total energy density of the cosmic fluid. But for the fluid with several components, we 
can equally well define the uniform density slicing for each component. Then if there is 
no total energy exchange between these components, Eq. (2.150) can be rewritten as 

Cn = -H 6 -^ = \^, (2.151) 

Pn 3 p n + p n 

where n is the index for a particular component of the fluid and 5p n is the energy density 
perturbation of that component on a flat slicing. Using Sp = ^2 n Sp n we can calculate 
the total curvature perturbation from Eq. (2.150) 

/' • P 

This equation will be important when we consider curvaton models where the primordial 
perturbation can be generated by several fluids. 
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2.4. Mechanisms for the Generation of the Curvature 
Perturbation 

Sections 2.2.2.5 and 2.2.2.6 described how, in the inflationary Universe, quantum fluctua- 
tions are amplified and converted into classical field perturbations 5(f). Then in section 2.3 
we have shown how to calculate the intrinsic curvature perturbation of the space-time 
which is measured after the horizon entry. In this section we will connect those two parts 
and show three mechanisms by which fluctuations of quantum fields during inflation can 
generate the curvature perturbation (. These three models of the generation of the cur- 
vature perturbation by no means are the only possible. However, only these three are 
necessary for our purpose when we discuss vector fields in Chapter 3. 

2.4.1. Single Field Inflation 

Let us assume in this section that the single field which drives the slow-roll inflation, 
as discussed in section 1.5, is the same field which is responsible for the total curvature 
perturbation in the Universe. 

During single field inflation the value of the field (ft (i, x) at any given instant determines 
the energy density p(t,x). Therefore, we can calculate the curvature perturbation £ 
directly from Eq. (2.150) by using the expression for the energy density and pressure 
of the scalar field in Eqs. (1.49) and (1.50). Imposing the slow-roll condition for which 
p ~ V ((f)) and 3H(j) ~ — Va, we find to the first order 

3 cf) 2 

Alternatively we can use the SN formula directly. This method renders the second or- 
der calculations more straightforward and comparison with other methods for generation 
of (" easier. 

In Eq. (2.140) and are derivatives of the number of e- folds of expansion of the 
unperturbed Universe. During single field inflation the local evolution of the Universe 
is determined only by the value of a single scalar field <f>. Therefore, the change of (f) 
corresponds to the shift in time along the same unperturbed trajectory, Nj, and N^s are 
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independent on the final epoch, which makes £ independent on time 6 

C (x) = N+Stj) (x) + N H [5<P (x)] 2 . (2.154) 

The definition of the number of unperturbed e-folds of expansion is given in Eq. (1.57), 
from which it follows that N = —H, or alternatively 

H IV 

AT = - = (2.155) 

(j) m£j Vcj, 

where for the last equality the slow-roll equation of motion in Eq. (1.54) was used. 
Inserting this expression into Eq. (2.154) at the first order one recovers the same equation 
as in Eq. (2.153). 

The field perturbations 8<pk was discussed in section 2.2.2.5, and the power spectrum 
for a light scalar field V$ was calculated in Eq. (2.106): 

M*)- (§) 2 - ( 2 - 156 ) 

A few Hubble times after horizon exit £ becomes constant as shown in Eq. (2.133), hence 
it is enough to evaluate (2.144) at this time, which gives 



2 



Using this expression for the power spectrum of the scalar field and Eq. (2.155) with 
3mpji/| ~ V (<£)|fc gives us the power spectrum of the curvature perturbation 



(2.158) 



where e is a slow-roll parameter defined in Eq. (1.55) and the right hand side has to be 
evaluate at horizon exit, aH^ = k. 

For exponential inflation, when H = 0, is scale independent. But in the slow- 
roll inflation H is only approximately constant, which makes the power spectrum slowly 
varying with k. This variation usually is approximated by a power law (c.f. Eq. (2.14)) 
and is parametrized by the spectral index n such that V( (k) oc k 71 ^ 1 . One can take the 



The situation is different for the multifield inflation, where the change in the field space, <p = 
(4>i, 4>2, ■ ■ ■), does not only correspond to the shift in time along the unperturbed trajectory, but 
also the rotation in this space. Then £ becomes time dependent until the end of inflation or until the 
trajectories of fields 0i, fa, ■ ■ ■ become straight lines. 
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definition of the spectral index to be 



n — 1 



dhxV, 



C 



d In k 

To evaluate this equation it is useful to derive the following relations 



(2.159) 



d In k = d In (afciJfc) — d In a/% = H^dt, 



(2.160) 



from which we can write 



^^i_ 
~N4>d0' 

where Eq. (2.155) was used. Then one can readily calculate 



d _ (^<i_ 
dlnfc ~ Hdcp 



(2.161) 



dlnU _ 1 _ 2 /% 
dlnA; " iV^V" mpi U 



-2e, 



(2.162) 



where we have used Eq. (2.155) and the definition of the slow-roll parameter in Eq. (1.55). 
Now, we also have 



dine (V 



d In k 



V& dmk V V 



V, 



2mp! 



V 



4e - 2ry. 



(2.163) 



Using the last two relations and the power spectrum in Eq. (2.158) from the definition of 
the spectral index in Eq. (2.159) we find that the spectral index for single field slow-roll 
inflation is 

n - 1 = 2r/ - 6e. (2.164) 



From Eq. (2.155) we may write 
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The number of e-folds of inflationary expansion when the observable cosmological scales 
leave the horizon correspond approximately to N ~ 10 [5]. Let us denote the change in 
the field value during this period by Ac^xch Considering that the low-roll parameter e is 
almost constant during this period from Eq. (2.165) we find 



e ~ 10" 



fA(j> 



10 



(2.166) 



which is the value of e when the cosmological scales leaves the horizon. Observational 
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constraints on n where discussed in section 2.1.2.1 Eq. (2.16), i.e. n — 1 ~ 0.04, which 
gives 2r; — 6e ~ 0.04. Therefore, from Eq. (2.166) we see that for small field inflationary 
models, for which A0io <C mpi, the slow-roll parameter e is much smaller than r\ and 
Eq. (2.164) may be written as 

n- 1 = 277. (2.167) 

One can go further and consider the running of the spectral index as well, i.e. the scale 
dependence n = n (k). The running is defined as n' = dn/d In k, from which one finds 

n' = -Wen + 24e 2 + 2f, (2.168) 

where 

Z = m Pl ^^, (2.169) 

where = d 3 V/d(p 3 . 

So far we have discussed the two point correlation function of Eq. (2.143). If perturba- 
tions are exactly Gaussian all higher correlators must vanish. Indeed, the non-Gaussianity 
of the curvature perturbation generated by the inflaton field is very small in single field 
inflation and this can be seen by calculating the three point correlators. To show this let 
us calculate the second derivative of Eq. (2.155) 

^ = 2e-r). (2.170) 

Then using Eq. (2.145) we find that non-Gaussianity is of order the slow roll parame- 
ters: 

^/ NL = r/-2e, (2.171) 
5 

were we have used Eq. (2.170). From this relation it is clear that |/nl| 1. The 
cosmic variance limits the detection of /nl to the values |/nl| > 3 [42]. Therefore, 
non-Gaussianity produced by a single field inflation is too small to ever be observed. 
However, there are other scenarios of the generation of the curvature perturbation for 
which non-Gaussianity can be large enough to be observable in the near future. 

One may consider even higher order correlators of the curvature perturbations. And 
in some models they can be large enough to be observable as well. But this is beyond 
the scope of this thesis. 
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2.4.2. At the End of Inflation 

In the previous section we discussed the scenario when the inflation is driven by a single 
scalar field <j>. In this scenario inflation ends when the inflaton field reaches a critical 
value 4> c , where the slow-roll conditions in Eqs. (1.55) and (1.56) are violated, which is 
solely determined by the inflaton field. Therefore, inflation ends on the uniform energy 
density slice. 

In Ref. [58] a scenario was suggested where the critical value 4> c is modulated by some 
other scalar field a, <j) c = </>c (c)- For single field inflation, the contribution from a 
to the inflaton dynamics must be negligible. <j) c in this case must depend only on the 
perturbation 5a (x). Then the hypersurface of constant (j) c does no longer coincide with 
the uniform density hypersurface. 

For clarity let us assume for the moment that the perturbation 5(f) of the inflaton field 
generated during inflation is negligible. Then, if the inflaton is a free field, the slice of 
constant (j> will coincide with the flat slice and with the constant energy density slice 
even at the end of inflation, i.e there will be no curvature perturbation. But if the end of 
inflation value (f) c is modulated by some other perturbed field, cp c (a (x)), then the uniform 
density slice at the end of inflation no longer coincides with the flat slice. According to 
the section 2.3.4, this produces a perturbation in the amount of expansion between the 
flat and uniform energy density slices which, from Eq. (2.138), is equal to the curvature 
perturbation ( en d- 

If, on the other hand, 6(f) is not negligible, then the same argument holds, but C en( j 
will correspond to the perturbation in the amount of expansion between the uniform 
energy density just before the end of inflation and the one just after the end of inflation. 
However, £ en d may be still large enough to dominate the curvature perturbation which 
is generated during inflation, £ ent } ^> (j n f. 

For the following discussion we will assume that after inflation the Universe undergoes 
prompt reheating, i.e. the inflaton field energy is promptly converted into radiation. 
Then from Eq. (2.140) up to the second order we can write 

Cend = NMc + N cc {5(f> c ) 2 , (2.172) 



where we have defined 



ON d 2 N 
N cS - and «„ = (2.173) 



In the former expression, 54> c is the perturbation of the end-of-inflation field value due 
to the coupling to the field a. Because, as was argued before, cp c depends only on the 
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perturbation 5a, we can expand to the second order 

6<j> c = <lJ c 6<T+±<%(6v) 2 , (2.174) 

where <p' c = d(j) c /da and cj)" = d 2 cj) c /da 2 . Keeping terms only to the second order in 5a 
the curvature perturbation £ enc j becomes 

Ce„d = N c <t>' c 5a + X - [N c ^ + N cc (</>' c ) 2 ] (5a) 2 . (2.175) 

To illustrate this model in a concrete example, let us consider the hybrid inflation 
scenario [59, 60, 61] first. In this scenario the slowly rolling inflaton field is coupled to 
another scalar field, called the waterfall field. The potential for this type of models is 
given by 

V (0, X ) = V (<P) - l -m 2 xX 2 + ^A X 4 + - VV, (2.176) 

where is the slow rolling inflaton field and \ is the waterfall field. From this expression 
it can be seen that the effective mass of % is 

m 2 eS = A 2 - m 2 x . (2.177) 

Initial conditions are such that m 2 s > and the waterfall field is located at x = while 
the inflaton field <j> slowly rolls towards zero. With this configuration the dominant term 
of the potential in Eq. (2.176) is the term V (4>). 

Inflation ends when the inflaton field reaches a critical value (f> c and the effective mass 
in Eq. (2.177) becomes negative. This destabilizes the waterfall field which very rapidly 
rolls down to the minimum of the potential and acquires the vacuum expectation value. 
This also very promptly changes the evolution of the inflaton (f>: instead of slowly rolling 
it is quickly driven towards zero. 

For the end-of-inflation scenario the critical value <f) c is additionally modulated by 
including one more field a, which is coupled to the waterfall field. The potential in 
Eq. (2.176) then becomes 

V (0, X ,o) = V (0) - \m 2 x X 2 + ^A X 4 + \\^ 2 X 2 + \ko 2 x 2 + V(a). (2.178) 
It is clear that the effective mass of x becomes 

m 2 cS = \ <i> 4> 2 + \ a a 2 -m 2 x . (2.179) 
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In this case the waterfall field is destabilized and inflation ends when 



m. 



(2.180) 



As we can see, the critical value <f> c is a function of another field, <p c (a). The first and 
second derivative with respect to this field are 
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(2.181) 
(2.182) 



In Ref. [58] it was considered that the curvature perturbation generated at the end of in- 
flation dominates over the one generated at the horizon exit. This happens if N c <p' c 3> Nj, , 
where N& is given in Eq. (2.155). Using Eq. (2.181) this condition can be rewritten as 



K_o_ 



(2.183) 



where e c is the slow-roll parameter just before the end of inflation. iV c in this expression 
was taken from Eq. (2.155) giving = 1/ (2mp 1 e c ). 

With £ en d dominating over £j n f from Eq. (2.175) we obtain the power spectrum of the 
produced curvature perturbation in this model as 



2mf,je c 
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(2.184) 



where a en d is the value of a at the end of inflation. Since the potential of the field a 
is flat during inflation, we may apply calculations in section 2.4.1 for its perturbation 
spectrum. Therefore, using Eq. (2.167) we find the spectral tilt as 



n 



1 = 2 Va , 



(2.185) 



where r) a is the second slow roll parameter for the field a. 

The non-Gaussianity parameter for this model when £ en d is dominant can be calculated 
as follows. From the expression of /nl m Eq. (2.145) and the curvature perturbation in 
Eq. (2.175) one finds 
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(2.186) 
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Since we already know that N cc /N^ is of the order of slow roll parameters, the first term 
dominates. And using Eqs. (2.181) and (2.182) we arrive at 



Expressions of the power spectrum and non-Gaussianity parameter /nl in Eqs. (2.184) 
and (2.187) involve the homogeneous mode of the field a. This mode is defined as the 
average value of the field, where the averaging is done over the same comoving box in 
which perturbations are defined. The comoving box size L must be larger than the 
observable Universe, but not too large [5, 62, 35]. As discussed in subsection 2.1.1, 
depending on the accuracy required the box should be such that In (LHq) ~ O (1). 

In single field inflation the value of the field at horizon exit may be calculated from the 
number of e-folds of remaining inflation. However, in general it is not possible to calculate 
the unperturbed value of the field and it must be specified as the free parameter of the 
model. In some cases, it can be evaluated using the stochastic formalism and assuming 
that our Universe is the typical realization of the whole ensemble. 

2.4.3. The Curvaton Mechanism 

In section 2.4.1 we have demonstrated the mechanism for the generation of the curvature 
perturbation at the horizon exit during single field inflation. In section 2.4.2 the curvature 
perturbation was generated at the end of inflation. Here we will discuss a mechanism by 
which ( is generated some time after inflation when the Universe is radiation dominated. 
Such model is called the curvaton model and was first introduced in Refs. [63, 64, 65, 
66]. It is possible that the total curvature perturbation is generated by the curvaton 
mechanism, or only a part of it. For simplicity, let us assume first, that the curvature 
perturbation is generated only by the curvaton mechanism. 

In these models the field which is responsible for the curvature perturbation is different 
from the field which drives inflation. In fact, it is not even necessary to assume any 
particular model of inflation and validity of Einstein gravity during that era: it might be 
slow roll inflation due to scalar fields, due to modified gravity or any other mechanism. 
The only assumptions necessary are that inflationary expansion is almost exponential and 
that after inflation the Universe undergoes reheating and becomes radiation dominated. 

The curvaton mechanism liberates inflation models from the need for the inflaton field 
to drive inflationary expansion as well as generate the primordial curvature perturbation. 
Therefore, it substantially increases the available parameter space for viable inflation 




(2.187) 
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models [67]. 

Let us denote the scalar curvaton field by cr. Although it is a different field from the 
one in the previous section, in this case too cr during inflation is subdominant with a 
sufficiently flat potential, \V aa \ <C H 2 , where V aa = d 2 V/da 2 . The unperturbed curvaton 
field satisfies equation of motion 

a + 3H& + V a = 0. (2.188) 

As for all light fields, quantum fluctuations of the curvaton field during inflation are 
promoted to classical perturbations after horizon exit. Then the curvaton perturbations 
5a (x) satisfy 

5a + 3Ha + V aa 5a = 0, (2.189) 

where we have taken into account that on superhorizon scales all gradients vanish since 
k/a —J- 0. With vacuum initial conditions the power spectrum for the curvaton field 
perturbations is given by Eq. (2.106) 

P„=(§) 2 . (2.190) 

At some epoch after inflation when V aa ~ H 2 (t) the curvaton field starts to oscillate 
around its VEV. Let us assume that at this epoch Einstein gravity is already valid and 
that this happens during the radiation domination (the energy density of the dominant 
contribution, i.e. radiation, decreases as p 7 oc a -4 ). Then 3rrip^H 2 = p 7 and the Hubble 
parameter decreases as H oc cT 2 . 

The onset of oscillations of the curvaton field is taken to occur much before the cos- 
mological scales enter the horizon. The potential of the curvaton near its VEV can be 
approximated as 

V{a)= l -m 2 a a 2 . (2.191) 

Then oscillations start when H ~ m a . However, even if the potential is not of that form 
at the start of oscillations, after a few Hubble times, when the amplitude decreases, it 
can be approximated to high accuracy by this quadratic form. Then, from Eqs. (2.188) 
and (2.189), it is clear that the unperturbed and perturbed values of the curvaton field 
satisfy the same equation of motion resulting in da/a = const. But even if this condition 
is not satisfied, it will result only in a scale independent factor which does not spoil scale 
invariance [65]. 
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During oscillations, the curvaton field evolves as the underdamped harmonic oscillator 
with the energy density decreasing as p a (t, x) im 2 (7 A (i, x) oc a -3 (i,x), where <t a is 
the amplitude of oscillations. This decrease is slower than that of radiation. Therefore, 
the relative energy density of the curvaton increases, p a /p-y oc a. If the curvaton decay 
rate is small enough it can dominate the Universe, resulting in the second reheating at its 
decay. Or the curvaton can decay when it is still subdominant. In both cases during the 
period when the relative energy density of the curvaton field increases the total pressure 
of the Universe is not adiabatic. According to Eq. (2.134) this results in a growth of 
the curvature perturbation £ which settles at its constant value soon after the curvaton 
decay. 

To calculate the power spectrum and non-Gaussianity of the curvature perturbation 
we use the 5N formalism (see Ref. [57]). This can be applied with the sudden decay 
approximation. The curvature perturbation without this approximation can be calcu- 
lated using Eq. (2.134) and knowing the decay rate of the curvaton. Such calculation 
can only be done numerically. However, in Refs. [68, 69] it was shown that the sudden 
decay approximation agrees with the numerical results within 10%. 

As we have assumed that the curvature perturbation in the radiation dominated Uni- 
verse is negligible before the curvaton starts oscillating, the number of e-folds from the 
end of inflation until oscillations is unperturbed. Therefore, the initial epoch for the 5N 
formula can be taken just after oscillations commence. Let us denote the value of the 
curvaton field at the onset of oscillations by <7 OS c (x). In general this will depend on the 
value of the curvaton field cr* a few Hubble times after horizon exit of a given scale k, 
i.e. o"osc = fosc (f*)- Then the number of e-folds from the beginning of oscillations until 
the curvaton decay is 

N (pdec,Posc,o-*) = ln°^ = ~ln^^, (2.192) 

Oosc o At, dec 

where 'osc' and 'dec' denotes values at the start of curvaton oscillations and at the decay 
respectively. On the other hand, for the energy density we have p osc (x) ~ ^m^.a^ sc , 
which gives 

1 -m 2 <7 2 

N G°dec, Pose, O"*) = o ln " — (2.193) 

" Pa, dec 

Let us denote the total energy density at the start of oscillations p osc . Because the curva- 
ton energy density at this epoch is negligible p sc corresponds primarily to the radiation 
energy density p osc ~ p 7)OS c- But the curvaton is not negligible just before its decay, mak- 
ing the total energy density at this epoch p& ec = /? 7i dec + Pa dec- From the scaling laws 
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of matter and radiation (see Eq. (1.8)) we find p a ,dec/ Pa,osc = (p-y,dcc/ Vy.osc) ■ Putting 
everything together we find 



2 2 I /'dec /V,dcc 



3/4 

Ar.dec = ^"V^osc ■ (2.194) 

2 V Pose / 

The derivative of N with respect to the curvaton field may be found using the chain 
rule d/da* = a' osc ■ d/da osc and keeping p osc and p^ec fixed (the prime in this relation 
denotes differentiation with respect to a*). Then from Eq. (2.194) we find 

da osc P(T,dec 8 Pdcc ~ Pa,dec ^ ^gj^ 

Pa,dec &osc ^Pdec Per, dec 

Using this relation, from Eq. (2.193) we calculate 

N a , = (2.196) 

" "Vise 



where 



Cl„ = „ 3pa ' de A C • (2.197) 

•iPa,dec "T ^P^^dec 



If the curvaton energy density is subdominant at the decay, then p a ,dec "C P7,dec> giving 
Cl a = where Q a = p a ,dec/ Pdec is the density parameter of a at decay. If, on the 

other hand, the curvaton is dominant at that epoch, then Cl a = Q a = 1. In both cases it 
is a good approximation to write Ct a ~ CT . The error introduced by such approximation 
is not bigger than that of the sudden decay approximation [5]. 

Inserting Eq. (2.196) into the equation (2.140) we find that the curvature perturbation 
generated by the curvaton mechanism to first order is 

C CT = \si a ^6<r*. (2.198) 



3 (T 



The power spectrum becomes 



* - »- (§) = >° it) (§) • (2 - 199) 

where for the curvaton field perturbation the power spectrum was given in Eq. (2.190). 
To find the non-Gaussianity parameter /nl for this model we have to know Q a at least 
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up to second order. Calculating the second derivative of Eq. (2.196) we find 

*/ NL = -^ = -2 + ^ + ^(l + ^). (2,00) 

Eqs. (2.199) and (2.200) were calculated using the SN formalism and they agree very 
well with the calculations performed using the first and second order perturbation theory 
in Refs. [57, 63]. 

In the previous calculations we have cr osc (a*) as a general function. In some models 
it might happen that this is a highly non-trivial function (for example Ref. [70]) but for 
the future reference let us consider the case when cr OS c — Then the power spectrum 
in Eq. (2.199) becomes 

(2 - 201) 

And the non-Gaussianity from Eq. (2.200) is 

i** = h < 2202) 

where we have considered the case when /nl 3> 1- 

The result shows that in the curvaton scenario the non-Gaussianity can be very large, 
/nl 3> 1- This is in contrast to single field inflation, where /nl was of order of the 
slow-roll parameters (see Eq. (2.171)). In the curvaton case the non-Gaussianity can be 
large because N a and N aa have nothing to do with the slow roll parameters. 

There is another notable difference from single field inflation - the power spectrum of 
the curvature perturbation in Eq. (2.201) depends on the homogeneous mode of the field. 
This situation is analogous to the one discussed in subsection 2.4.2. The value of cr* (x) 
should be taken as an average within the observable Universe. Then it can be calculated 
assuming that our Universe is typical. 

Until now we have considered only the case when the curvaton is the only source of 
the curvature perturbation Q. In other words, we have assumed that the perturbation in 
the Universe is negligible prior to the domination (or near domination) of the curvaton. 
But if we drop this assumption then Eq. (2.152) can be used to calculate the resulting 
curvature perturbation from both components: 

c = , Ap T; — c 7 + , 3p Tl — (2.203) 

"At, dec ~r ^Py^dec "At, dec "T ^P-f,dec 

where £ 7 is the curvature perturbation in the radiation dominated background. Using 
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the definition in Eq. (2.197) this equation becomes 

C = (l - C 7 + JkCr- (2-204) 

Such curvaton models with a two component contribution to the total curvature perturba- 
tion were considered in Refs. [71, 62, 72]. While for the negligible curvature perturbation 
from the inflaton this equation reduces to 

C ~ n a ( a . (2.205) 
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Perturbation from Vector Fields 



3.1. Vector Fields in Cosmology 

In the previous chapter we have discussed the generation of the primordial curvature 
perturbation in the Universe. It was shown that during inflation the quantum mechanical 
fluctuations of light scalar fields are transformed into the classical curvature perturbation. 
Long after inflation upon horizon entry this perturbation seeds the formation of large 
scale structure. Until very recently the generation of the curvature perturbation by this 
mechanism was assigned solely to scalar fields. The main reason for this is that scalar 
degrees of freedom are the simplest ones. The preference for scalar fields is further 
supported by the observational fact that the Universe on large scales is predominantly 
isotropic and statistical properties of the temperature perturbation in the CMB sky are 
predominantly isotropic too. In addition, particle physics theories beyond the Standard 
Model are abundant with scalar fields. 

In the rest of this thesis we will show that quantum fluctuations of vector fields may 
influence or even generate the total curvature perturbation in the Universe. But why 
should we consider something else than the scalar field? The motivation comes from 
both sides: theoretical as well as observational. 

From the theoretical side, for the inflationary model building only scalar fields have 
been used to generate the curvature perturbation, even if no fundamental scalar field 
is discovered yet. Although it is widely accepted that all elementary particles possess 
masses due to the Higgs scalar field, this might be explained by other mechanisms, with- 
out invoking fundamental scalar fields (see e.g. the technicolor model in Ref. [73]). It 
is expected that in the near future the Large Hadron Collider (LHC) in CERN will dis- 
cover the Higgs boson and prove the existence of the fundamental scalar fields. But if the 
Higgs boson is not discovered, the alternative models explaining masses of elementary 
particles will become more favorable. On the other hand, in the case of inflation, the 
generation of the curvature perturbation from scalar fields will become much less attrac- 
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tive. Furthermore, even if scalar fields are discovered, particle physics theories, such as 
supersymmetry or supergravity, incorporate many vector bosons and fermions. However, 
despite this, the possible contribution to the curvature perturbation from other kind of 
fields usually has been ignored. 

Observationally there is some indication that the simplest scalar field scenario may not 
be sufficient to explain some features of the CMB sky. These features, which challenge the 
simple homogeneous and isotropic model of the Universe, were first discovered already 
after the release of the first year WMAP satellite data. For example it was found that 
the quadrupole moment of the power spectrum of the temperature perturbation was too 
small compared with predictions of the currently favored ACDM cosmology. The lack of 
power in the quadrupole was still present in the five year data (see Ref. [17]). Another 
discovery was that 2-4-8-16 spherical harmonics of the CMB temperature map seem to 
be aligned, suggesting the presence of a preferred direction in the Universe, so called 
the 'Axis of Evil" [74] . In addition this preferred direction is aligned with the large cold 
spot in the CMB [75], with the large void in the radio galaxy distribution [76] and with 
the galaxy spin directions [77]. Although currently these anomalies are under intense 
debate about their statistical significance, they might be an indication that the Universe 
is mildly anisotropic on large scales. 

If these anomalies are confirmed it will prove the existence of the preferred direction in 
the Universe and this can not be explained solely by scalar fields. On the other hand, for 
vector fields the existence of the preferred direction is natural. But employing vector fields 
for the generation of the curvature perturbation we encounter with two complications: 
conformal invariance and excessive large scale anisotropy of the Universe. 

In this Chapter it will be enough to approximate the inflationary expansion of the 
Universe to be exactly exponential, i.e. with the constant Hubble parameter H = 0. 

3.1.1. Conformal Invariance 

The evolution of the conformally flat space-time, such as de Sitter or matter /radiation 
dominated FRW universes, can be modeled as the conformal rescaling of the Minkowski 
space-time, i.e. by transforming the metric g„v (x) — > a 2 (r) g^ u (x), where r is the 
conformal time. As we have seen in section 2.2.2 for a light, minimally coupled scalar 
field on superhorizon scales this leads to the amplification of vacuum fluctuations. But 
this is not the case for the conformally trivial theories, for which the field equations 
are invariant under the rescaling of the metric. For such theories the form of the field 
equations is time independent in the conformally flat space-times. 
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This is the case, indeed, for the massless U(l) vector field with the Lagrangian 

C = -lF^F^, (3.1) 

where F^ v is the field strength tensor 

Ffxu = d»A v - d v Ap, (3.2) 

and is the four-vector. In contrast to the light scalar field case discussed so far, 
quantum fluctuations of a vector field with this Lagrangian do not undergo amplification 
after the horizon exit. But if the vector field is to generate the non-negligible curvature 
perturbation, its fluctuations have to undergo amplification. Therefore, the conformal 
invariance of the U(l) massless vector field must be broken. 

This problem is well known in the literature on the generation of primordial magnetic 
fields (PMF) during inflation [78, 79, 80, 81, 82]. In this literature there are numerous 
suggested ways of breaking the conformal invariance for vector fields: (i) introducing a 
mass for the vector field, (ii) making the kinetic term time dependent, (iii) introducing an 
anomaly term (iv) coupling a vector field to another field which is not conformally coupled 
to gravity, (v) or using non-Abelian vector fields. In this thesis we use only the first two 
methods. In the context of the curvature perturbation the breaking of the conformal 
invariance of U(l) field by introducing the mass term was first considered in Ref. [83] 
and using the time dependent kinetic term in Ref. [84]. The first attempt to calculate 
the generation of the curvature perturbation by the non-Abelian SU(2) vector fields was 
reported in Refs. [85, 86]. The other methods of breaking the conformal invariance were 
investigated only for the generation of PMF. 

3.1.2. Large Scale Anisotropy 

If the vector field is to influence or generate the curvature perturbation, its energy density 
must dominate or nearly dominate the Universe for the effect to be non-negligible. But 
the energy-momentum tensor of the light vector field has an anisotropic stress. For 
example the energy-momentum tensor of the light Abelian vector field can be written as 
[83] 

T M „ = diag (p, -p ± , -p ± , +p±) , (3.3) 

where we have chosen the coordinate axis in a way that spatial components of the ho- 
mogeneous vector field are equal to A = (0, 0, A) and we use this choice of coordinates 
in the rest of this thesis. From Eq. (3.3) it is clear that if such a vector field is to domi- 
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nate or nearly dominate the Universe, the expansion along the vector field direction will 
be different from the transverse directions. This would induce an excessive large scale 
anisotropy which is ruled out by observations. 

To bypass this problem there are four methods proposed in the literature. The au- 
thor of the earliest one in Ref. [87] considered three orthogonal identical vector fields. 
In this model the total energy- momentum tensor, which is the sum of all three vector 
fields is isotropic. Therefore, these fields can dominate the Universe and even drive the 
inflationary expansion. Another mechanism, with vector fields responsible for the infla- 
tionary expansion, was proposed in Ref. [88]. The authors introduced a large number of 
identical vector fields which are randomly oriented in space with identical initial condi- 
tions. Due to random orientation, the average pressure becomes almost isotropic. The 
residual anisotropy is proportional to iV -1 / 2 , where N is the number of vector fields. If 
this number is sufficiently large, the induced large scale anisotropy can be small enough 
to agree with observational bounds. In another model in Ref. [89] the vector field is 
always sub dominant and therefore does not generate excessive large scale anisotropy. It 
cannot be responsible for the inflationary expansion, but the vector field influences the 
generation of the curvature perturbation by coupling to the scalar field. In particular, 
the authors of this paper consider that the vector field modulates the end of inflation. 
Finally in Ref. [83] the excessive large scale anisotropy is avoided by introducing the 
vector curvaton scenario. In this scenario, a massive vector field is subdominant during 
inflation and afterwards until it becomes massive. As was demonstrated in that work, 
when the vector field becomes massive it starts to oscillate with a frequency much larger 
than the Hubble time. The pressure components in Eq. (3.3) induced by the oscillating 
vector field oscillates rapidly themselves. Therefore, the time averaged value of the pres- 
sure over one Hubble time is zero and the vector field on average acts as the pressureless, 
isotropic matter. It can dominate the Universe without generating excessive large scale 
anisotropy. 

In this thesis we will use the end-of-inflation scenario of Ref. [89] as an example to 
calculate the statistical properties of the curvature perturbation (section 3.6). But mostly 
we will be occupied with the vector curvaton scenario (section 3.3). 

3.1.3. The Physical Vector Field 

Before going into the description of the vector field quantization and the generation of the 
curvature perturbation let us make a comment about the distinction between the field 
which appears in the Lagrangian and the physical vector field. Consider, for example, 
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the Lagrangian of the massive Abelian vector field 1 

C = -i^i^ + \m 2 A»A^ (3.4) 

with the field strength tensor Fu V defined in Eq. (3.2). Using the FRW metric the mass 
term may be expanded as 

C = -\F, V F^ + \m 2 A 2 - imV 2 A 2 , (3.5) 

from which one notices that the spatial term dependents on the scale factor. This might 
be alarming, because the Lagrangian is the physical quantity and cannot depend on the 
arbitrary choice of the normalization of the scale factor. This shows that the vector 
field A in the Lagrangian is the comoving field defined with respect to the comoving, 
Cartesian space coordinates x l . The physical four- vector field is 

W„ = (A), A/a) , (3.6) 

which is defined in in the basis of the physical coordinate system a (t) x l . For the comoving 
and physical vector fields the corresponding upper-index quantities are A % = —a~ 2 Ai and 
W* = -Wi. 

This is the case for the Lagrangian with the canonically normalized field such as in 
Eq. (3.4). In section 3.5 we will consider a vector field with the time dependent kinetic 
term. In this case will denote a canonically normalized physical vector field (see 
Eq. (3.174)). 

3.2. Vector Field Quantization and the Curvature 
Perturbation 

3.2.1. 6N Formula with the Vector Field 

In this section we will generalize the 5N formalism introduced in section 2.3.4 to include 
perturbations of the vector field. For simplicity we will assume that only one perturbed 
vector field affects the local expansion rate in the Universe. And keeping one scalar field 



To be more precise, the Lagrangian of the Abelian vector field whose gauge symmetry is broken by an 
explicit mass term. 
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we can write to the second order (c.f. Eq. (2.141)) 



C(*,x) 



<JiV(^(x),Wi(x),t) 

N+Sf + N l w 5W l + ^N H (5(b) 2 + % w 5cP5Wi + ^iV^aWitfW,- + (3.7) 



where 




dN 




d 2 N 



d 2 N 



Ni 



d 2 N 



(3.8) 



dWi 



dcj) 2 ' 



w — 



dWidWj 



<j>w = 



with i and j denoting spatial indices running from 1 to 3. As with scalar fields, the 
unperturbed vector field values are defined as averages within the chosen box (see the 
discussion below Eq. (2.187)). 

For this expression there is no need to define Wi as components of the vector field. 
Even more, this expression is valid not only for the isotropic background expansion, but 
for anisotropic as well. Although for the aim of this thesis it will be enough to consider 
spatially flat isotropic geometry with the line element in the conformal time 



3.2.2. The Vector Field Quantization 

To quantize the vector field, let us expand perturbations of the field in Fourier modes, 
similarly to the case of the scalar field 



The massive vector field has three degrees of freedom, and the massless vector field 
has two, in contrast to the scalar field case which has only one degree of freedom. In 
Eq. (3.10) we have included only spatial components of the vector field, because the 
temporal component is non-dynamical, i.e. it is not a degree of freedom and is related to 
the spatial components through the equation of motion (for the massive vector field) . The 
perturbation of each degree of freedom may be parametrized using polarization vectors 
as 



A 

where e\ are polarization vectors, k = h/k is the unit vector in the direction of k and 
k is the modulus k = |k|. The most convenient choice is the circular polarization for 




(3.9) 




(3.10) 




(3.11) 
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which two transverse vectors have different handedness. Because both of them transform 
differently under rotations, the rotational invariance of the Lagrangian prevents any 
coupling between them. Choosing the coordinate z axis to point into the direction of k, 
the circular polarization vectors ef take the form 

,L _ n „• n\ /./o „R 



= (1, i, 0) /V2, ef = (l,-*,0)/V2 and el 1 = (0,0,1)- (3.12) 

In these expressions superscripts 'L', 'R' and '||' indicate the left-handed, right-handed 
and longitudinal polarizations respectively. For the massive vector field all three polar- 
izations are present, but for the massless one wu = in Eq. (3.11). These expressions 
define polarization vectors only up to a rotation about the k direction, but this is enough 
for the present purpose. Under the transformation k — > — k one of the axis x or y change 
the sign as well. We will choose that x changes the sign and y stays the same. Then 
e\ ^— ~kj = —e\* and because 5W{ (r, x) is real, imposing the reality condition onto 
Eqs. (3.10) and (3.11) gives w* x (-k) = -w\ (k). 

Later it will be useful to rewrite polarization vectors in Eq. (3.12) when k points to 
an arbitrary direction, not only along the z axis. Using the Cartesian coordinate system 
they become 

e i(k) = k = {k x , ky, k z ) , 

e^(k) = 1 (~ky + ik x k Z , k X + ikyk Z , ~i (kl + kl)), (3.13) 

V 2 ( fc «+*«) 

ef(k) = e^(k). 

Quantization of the vector field proceeds in the same way as for the scalar field: we 
expand each degree of freedom in a complete set of orthonormal mode functions and 
promote the expansion coefficients to operators with appropriate commutation relations. 
Vector field components satisfy the Klein-Gordon equation and the field lives in the 
homogeneous and isotropic FRW space-time, therefore, the complete set of orthonormal 
mode functions was already chosen to be e lk ' x in Eq. (3.10). The quantized vector field 
then takes the form 2 

6W i = J2[ k A (k)^ A (r, k) a x (k) e* k ' x + e}*(k)w* x (r, k) a[ (k) e^H -^L (3.14) 



2 Note that in this chapter we have changed our normalization of Fourier modes from (2ir) 3//2 to 
(27r)~ 3 . 1 
operators 



(27r) 3 . This resulted in the (2-7r) 3 factor in the commutation relations for creation and annihilation 
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where 

[a x (k) , ay (k')] = (2vr) 3 5 (k - k') 5 XX , (3.15) 

and other commutators being zero. As in the scalar field case, after the horizon exit vector 
field perturbations become classical in the sense that the commutator dWi (x) , d T 5Wj (x') 
approaches zero. 

In later sections we will discuss several mechanisms to generate scale invariant pertur- 
bations of the vector field. In all of these mechanisms perturbations will be Gaussian, 
with no correlation between different polarizations A or between perturbations of scalar 
and vector fields. With these conditions it is sufficient to consider only the spectra of 
vector field perturbations V\ = V Wx - They are defined by the analogue to Eqs. (2.3) and 
(2.4) as 



2tt 2 

(w x (k) w* x , (k')> = (2^) 3 S xx ,6 (k - k') - W V X (k) , (3.16) 



(w x (k) w x > (k) > = - (2vr) 3 5 XX ,5 (k + k') ^P A (k) , (3.17) 

where V x (k) = -Pa (k) k 3 / (27r 2 ) as defined in Eq. (2.6), and we have suppressed the 
notation of time, i.e. w x (r, k) = w x (k). 

If the Lagrangian is parity conserving then Pl = Pr, which will be the case in all 
models considered in this thesis. Parity violation is introduced by terms involving the 
dual of Fnu, i.e. F^ u = e^ upcT F pa , where e nvp(J is the totally antisymmetric tensor. 
Examples of such theories may be found in Refs. [90, 91, 92]. The difference between the 
left-handed and right-handed power spectra would indicate a parity violation. Therefore, 
it is convenient to define parity conserving P + and parity violating P_ power spectra by 

Pr ± Pt 

P± = • (3.18) 

We also define two parameters which quantify the anisotropy in the particle production 
of the vector field 

p{k)= v + {k) — and q{k) = T^Wy (3 - 19) 

where we have also assumed that the expansion during inflation is isotropic making V\\ 
and V± dependent only on the modulus of k. By the isotropic particle production it is 
meant that the perturbation spectrum for all three degrees of freedom is the same, i.e. 
p = and (7 = 0. In this case the curvature perturbation generated by the vector field 
does not differ from the scalar field. But if p ^ and/or q ^ the particle production of 
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■ 




Figure 3.1.: Random field maps. The first is statistically isotropic random field. If the 
vector field undergoes isotropic particle production with p = and q = 0, the 
curvature perturbation generated by such a vector resembles the first map. If, 
on the other hand, p ^ and/or q ^ the curvature perturbation generated 
by the vector field is statistically anisotropic and resembles the second or the 
third map. The second map is of the statistically anisotropic random field 
with the preferred direction pointing vertically, while the preferred direction 
of the third map is horizontal. 



the vector field is anisotropic. If such a vector field generates the curvature perturbation, 
its statistical properties are not invariant under the rotations, i.e. it is statistically 
anisotropic (see Figure 3.1) 

Calculating the two-point correlators of 5Wi (k) we find 



(SWi (k) SWj (k')> = (2tt) 3 6 (k + k' 



a 2^ 
k 3 



elej*V h + efef*V R + e l le l jv l 



-i ~ 3 



, (3.20) 



where we have suppressed the notation of k for e^(k) and of k for V\ (k). This expression 
may be rewritten in terms of V± in Eq. (3.18) as 

27T 2 

(SWi (k) SWj (k') > = {2irf S (k + k') -p- x 

x T^(k)V + (k) + zT°f d (k)V- (k) + 7j ng (k)P|| (k)] , (3.21) 

where we have introduced tensors 



T- ven (k) = e, L (k) ej L *(k) + ef (k)ef* (k) = «£(k)ef (k) + ef (k) e J"(k), 
7- dd (k) = 



ef(k)ej*(k)-ef(k)ej* (k) =i ef(k)ef (k) - ef (k) ej L (k) , (3.22) 



lj ng (k) = 4'(k)e]((k). 
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With the circular polarization vectors derived in Eq. (3.13) these tensors take a simple 
form 



rpeven 
ij 



(k) 
t f d (k) 



5 



plong 

ij 



(k) 



y 
k% kj . 



k% kj , 



(3.23) 
(3.24) 
(3.25) 



3.2.3. The Power Spectrum 



Since £ is Gaussian to high accuracy, it seems reasonable to expect that £ will be domi- 
nated by one or more of the linear terms in Eq. (3.7). Keeping only them (corresponding 
to what is called the tree-level contribution) we find 



n(k) 



TfJ^{k)V + (k) + T^V\\{k) 



ATj^ (A;) + (k) + (N W ■ k)" [V\\ - V+) . 



(3.26) 
(3.27) 



Note that the power spectrum of £ is dependent on the direction of k. In the upcom- 
ing discussion we will frequently use the modulus of and the unit vector along its 
direction defined by 



N w = |N 



N^N^ and N w = 



(3.28) 



The curvature perturbation power spectrum (k) may be further separated into isotropic 
and anisotropic parts 



v c (k) = vr (k) 



l + g(k) (N w -k 



(3.29) 



which has the same form as anisotropic power spectrum in Eq. (2.7) keeping only up to 
the quadratic term. Comparing this expression with Eq. (3.27) we find that the isotropic 
part of the spectrum is 



Vf° (k) = Nfa (k) + (k) = NP* (k) 



(3.30) 
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where we have introduced the parameter £ 




(3.31) 



This parameter specifies the relative contribution from the vector field to the statistically 
isotropic part of the curvature perturbation. 

From Eq. (3.27) we can also find that the anisotropy in the curvature perturbation 
power spectrum is equal to 

V\\ (k) - V+ (k) £ 
g(k) = N$y " . , 7 = rm 77T1 p (k) . (3.32) 

If the vector field perturbation dominates £, i.e. £ 3> 1, the anisotropy in the power 
spectrum of the curvature perturbation is equal to the anisotropy in the particle pro- 
duction of the vector field g ~ p. As was mentioned in section 2.1.2.1 the observational 
bound for the anisotropy in the power spectrum of the curvature perturbation is g < 0.3. 
Therefore, if the anisotropy in the particle production of the vector field is larger than this 
bound and there is no other vector field contribution, the produced statistical anisotropy 
would violate observational constraints. To prevent this, the dominant contribution to £ 
must come from one or more statistically isotropic scalar field perturbations. 

Equations (3.26) and (3.27) for the power spectrum of the curvature perturbation were 
calculated only at the tree level. Analogous calculation for the one-loop contribution may 
be found in Ref. [93]. 

3.2.4. The Bispectrum 

Working to the leading order in the quadratic terms of the 5N formula, we calculate the 
tree-level contribution to the bispectrum defined in Eq. (2.13) 

B ( (k l5 k 2 , k 3 ) = (ki, k 2 , k 3 ) + B m (k l5 k 2 , k 3 ) + % (ki, k 2 , k 3 ) , (3.33) 

where we have separated the bispectrum into three parts: one due to perturbations in 
the scalar field, another part due to the vector field perturbations and the mixed term. 
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These terms are given by 

(ki,k 2 ,k 3 ) 
B<f,w (ki,k 2 ,k 3 ) 
B w (ki,k 2 ,k 3 ) 



NlN d . 



^ ^2 
pp V t + C -P- 



4tt 4 



V^Mi (k 2 ) + c.p. 



(3.34) 



4vr 4 

^1^2 



Mi (ki) iV^M,- (k 2 ) + c.p. 



The power spectrum (k) in the above equations depends only on the modulus of 
k because we assumed that the expansion during inflation is isotropic, and the vector 
Mi (k) characterizes perturbations of the vector field: 



Mi (k) = V+N w N l w + h k-% p + i kx% )? 



(3.35) 



Reversal of the three wave- vectors in Eq. (3.34) corresponds to the parity transfor- 
mation, and using the reality condition £ (— k) = Q* (k) we find that it changes each 
correlator into its complex conjugate. This does not affect the power spectrum because 
the reality condition also makes the spectrum real. This may not affect the isotropic 
bispectrum as well, because the reality condition and statistical isotropy make the bis- 
pectrum real. In our case, the bispectrum is anisotropic, and is guaranteed to be real 
only if the theory is parity conserving, i.e. if q = 0. 

The second order contribution of the quadratic terms in the SN formula gives the 
one-loop contribution to the bispectrum. It could be significant or even dominant. It 
has been calculated for the scalar case in Ref . [94] , and has been investigated for the case 
of multifield inflation in Refs. [95, 96] for example. The one-loop contribution from the 
vector perturbation is calculated in Ref. [97]. 



3.2.5. The Non-Linearity Parameter / NL 

In calculating the non-linearity parameter defined in Eq. (2.20) we will be interested in 
two configurations: equilateral, with k\ = /c 2 = k^, and squeezed, with k\ ~ A; 2 ^> k$. In 
the equilateral configuration the bispectra from Eqs. (3.34) become 

-N^wT* [Mi (ki) + Mi (k 2 ) + Mi (k 3 )] , (3.36) 
M(k 1 )iV^(k 2 ) + c.p. J 



fljr(k 1) k a> k 8 ) = 

^ Uil (k!,k 2 ,k 3 ) = 
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where we have defined for the equilateral configuration 
jequil^ u = f J^L\ R ec l uil 



B^ (k 1; k 2 , k 3 ) = -\ 5 ( equu (k l5 k 2 , k 3 ) , (3.37) 



and jg° qml = 6^ qui _|_ g^^il _|_ ^^ ul1 . In this case the non-linearity parameter /nx" 1 i s 
expressed using the power spectrum and the bispectrum as: 

6 , B^fkj.ka.ka) 

u j-equil _ C v " z ' 3; /o qo\ 



3 (^ so 

Observations give a limit on the anisotropy g < 0.3 (see the discussion above Eq. (2.19)). 
Therefore, since the anisotropic contribution to the curvature perturbation is subdom- 
inant compared to the isotropic one, we have included only into the expression of 

i-equil 
•/NL ■ 

In the squeezed configuration we have for the two vectors ki ~ — k 2 , but the third 
vector k 3 is of much smaller modulus than the other two and almost perpendicular to 
them. For this configuration Eqs. (3.34) take the form 

Bj^ki.ka.ka) = 2NlN (j>4> V^(k 1 )V (t> (k 3 ), 

B^Cki.kz.ka) = -JV^ w {^(*i)Mi(k3)+^(*;3)Re[M(ki)]}, (3.39) 
B^i^MM = 2Re[M l (k l )]N t ^Re[M j (k 3 )}, 

where Re [. . .] denotes the real part and ; g^ ocal (ki, k 2 , k 3 ) is defined similarly to Eq. (3.37) 

4 ocal (kl,k 2j k 3 ) = -j-lB^i^MM). (3.40) 

Then, the non-linearity parameter f^L^ m the squeezed configuration becomes 

6 . , ^ ocal (ki,k 2 ,k 3 ) 

^ ^ local _ C v Jy (o 41 \ 

5 ^nl 2Vl so (k 1 )Vf°(k 3 )' { ' ; 

In the treatment so far we have calculated the curvature perturbation generated by 
the vector field as well as the resulting anisotropic power spectrum and the non-linearity 
parameter /nl m equilateral and squeezed configurations. However, we didn't discuss 
how the perturbation of the vector field is generated and which mechanism transformed 
the field perturbation into the curvature perturbation. In the rest of the thesis we con- 
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sider several examples of the conformal invariance breaking for the vector field, which 
generates scale invariant perturbation spectrum and determines the values of p and q pa- 
rameters. We also consider two scenarios in which the vector field perturbation influences 
or generates the curvature perturbation (. 

3.3. The Vector Curvaton Scenario 

In section 3.1.2 it was discussed that a dominant light vector field would generate large 
scale anisotropy in the Universe which violates observational bounds. But to generate 
or influence the curvature perturbation by a vector field it has to dominate or nearly 
dominate. One of the ways to overcome this difficulty is the curvaton scenario. This 
scenario was summarized in section 2.4.3 with the curvaton acted by a scalar field. In 
this section we consider scenarios with the curvaton acted by a vector field, which was 
first introduced in Ref. [83]. 

3.3.1. The Vector Curvaton Dynamics 

Let us consider a massive Abelian vector field in the Universe dominated by matter 
(radiation) with the barotropic parameter w = (w = 1/3). The Lagrangian of the 
vector field is 

C = - l -F^ v + l -m 2 A»A^. (3.42) 

And let us choose a coordinate system in such a way that spatial part of the homogeneous 
mode of the vector field has components Ai = (0,0, A). In Ref. [83] it was shown that 
the energy momentum tensor for this field may be written as 

T l = &&g{pw,-p±.,-p±,+Pi.) , (3.43) 

where 

Pw 

with 

Pkin 

V w 



Pkin +V W , p± = pkin - V W 



(3.44) 



If f^- x -(- 

4 "2a 



l -m 2 A„A» = -m 2 (-\ 2 . 
2 M 2 Ui 



(3.45) 
(3.46) 
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From which we see that in general the energy momentum tensor will have anisotropic 
stress due to the opposite sign of the pressure components in the direction parallel to 
the field and the perpendicular one. Indeed, the equation of motion for the homogeneous 
mode of the vector field with the Lagrangian in Eq. (3.42) is given by [83] 



A + HA + m A = 0. 



(3.47) 



In the matter or radiation dominated Universe the solution of this equation is 



A (t) = f [dJ v {mt) + C 2 J- V {mt)] , 



(3.48) 



and the time derivative of the vector field is given by 



A (t) = — [Ci J„_i {mt) - C 2 Ji- v {mt)} , 
m 



(3.49) 



where C\ and C 2 are constants of integration, J v is the Bessel function of the first kind 
and 

v= ) + 3W , . (3.50) 
6(l + w) v 1 

In the FRW Universe the Hubble parameter is H ~ t~ l and the light field corresponds 
to mt < 1. In this limit Bessel functions can be approximated by power law functions 
and the growing mode of the vector field changes with time as 



A (t) oc t 



2v 



(3.51) 



Inserting this into Eqs. (3.45), (3.46) and using Eq. (3.43) it is clear that the energy- 
momentum tensor of the light vector field has anisotropic stress. If such a vector field 
dominated the Universe, the expansion rate in the direction of the field would be different 
from the transverse directions and the Universe would become predominantly anisotropic. 
Such excessive large scale anisotropy is forbidden by observations, therefore a light vector 
field cannot dominate the Universe. 

In the opposite regime, when the vector field is heavy, mt ^> 1, the Bessel functions 
may be approximated by trigonometric functions giving 



A{t) = t v 



irmt 



,, ■ i l-2v \ _ / l-2v 
L>\ sin ( mt H — ir I + C 2 cos I mt — ir 



A {t) 



-mt 



1 - %v \ _ . / l-2v 

Ci cos mt H — 7r — C 2 sm mt — tt 

irmt V 4 / V 4 



(3.52) 
. (3.53) 
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From these equations it is clear that the evolution of the heavy vector field resembles the 
evolution of the underdamped harmonic oscillator with decreasing amplitude of oscilla- 
tions. To see this let us rewrite Eqs. (3.52) and (3.53) in the form 



Ait) = At v \l — sin (mt + <p) , (3.54) 
irmt 



Alt) = -hmt v \l cos(mt + <p) , (3.55) 

irmt 



where constants A and cp are related to the original constants by 



A = yJCf + Cf + 2dC 2 cos (ttv), (3.56) 

and 



(p = arccos 



d cos (^tt) + C 2 sin (^tt) 



A 

Calculating the energy density and pressure from Eq. (3.44) we find 



(3.57) 



p w = a' 3 — A 2 , (3.58) 
p± = a" 3 — A 2 cos (2mt + 2p) . (3.59) 

TT 

Since the vector field is heavy, i.e. m 3> H, the frequency of oscillating functions in 
Eq. (3.59) is much larger than the Hubble parameter. Therefore, during one Hubble 
time the average pressure of the heavy vector field is zero and we can write 

pw oc a -3 and p± « 0. (3.60) 

Thus the heavy vector field acts as the pressureless isotropic matter and it can dominate 
the Universe without generating excessive large scale anisotropy. 

This property of the vector field is utilized in the vector curvaton scenario (see Fig- 
ure 3.2). During inflation the vector field is light, and although its energy-momentum 
tensor has non- vanishing stress in accordance with the curvaton scenario, it is subdom- 
inant, allowing the expansion of the Universe to be isotropic. During this period the 
vector field with broken conformal invariance undergoes particle production. (In sec- 
tions 3.4 and 3.5 we consider two mechanisms of breaking the conformal invariance of 
the vector field and producing a flat perturbation spectrum.) After inflation, the Hubble 
parameter decreases as H oc t . When it becomes smaller than the mass of the vector 
field, the latter starts to oscillate and, as was shown in Eq. (3.60), acts as the pressureless 
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m< H 



m> H 




log (a) 



Vector HBB 
decay 



Figure 3.2.: The schematic graph of the vector curvaton scenario: I. during inflation the 
vector field is subdominant but it acquires a scale invariant perturbation 
spectrum; II. inflation ends and the inflaton field starts oscillating around 
its minimum, while the energy density of the Universe scales as p oc a -3 (see 
section 1.5.3); III. after reheating the Universe is radiation dominated and the 
vector field is still subdominant until it becomes heavy; IV. the heavy vector 
field oscillates and acts as the pressureless isotropic matter with the energy 
density decreasing as pw oc a -3 ; V. the vector curvaton dominates (solid line) 
or nearly dominates (dashed line) the Universe, imprints its perturbation 
spectrum and decays, recovering the standard HBB cosmology. Note, that 
in this graph we depicted the situation when the vector field becomes heavy 
after reheating. However, this might happen before reheating as well. 
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isotropic matter. After reheating the relative energy density of the vector field increases 
as pw/Py °c a, where p 7 is the energy density of the radiation dominated Universe. As in 
the original curvaton scenario (see section 2.4.3) the curvaton field during the radiation 
dominated period dominates (or nearly dominates) the Universe, imprints its pertur- 
bation spectrum and decays, recovering the standard Hot Big Bang cosmology. If the 
curvaton decays when it is dominant, the Universe undergoes a second reheating. 

3.3.2. The Generic Treatment of / NL 

In this section we obtain analytic expressions for the non-linearity parameter /nl without 
assuming a specific vector curvaton model. In contrast to the original curvaton idea we 
include perturbations already present in the radiation dominated Universe when the 
vector field energy density is still negligible (for a similar study in the scalar curvaton 
case see Refs. [71, 72]). 

Some time after inflation the mass of the vector field becomes bigger than the Hubble 
parameter. Then the field starts to oscillate and acts as the pressureless isotropic matter. 
The total contribution to the curvature perturbation by the vector field before its decay 
can be found using Eqs. (2.151) and (2.204) 

Cvec = &w(w = o^iy-^ 4 -, (3.61) 
3 p w 

where £l\y is defined similarly to the Eq. (2.197) 

a _ 3p w 3JV , , 

= „ : — = a — (3.62) 

and energy densities pw and p 7 are evaluated at the curvaton decay with £lw = pw/p 
and p = pw + p 7 . This expression is valid to the first order in 5pw, which is evaluated 
on a flat slice, where a (t, x) is unperturbed. 

We assume that the curvaton decays instantly (sudden-decay approximation) and eval- 
uate Cw j us t before the curvaton decays, leaving £ constant thereafter. Evaluating 5pw 
to the second order we have [98] 

2 - Wi 1 ~ 1 

Cvec = -n w ^5Wi + -n w —SWiSWi, (3.63) 

where W = |W| is evaluated just before the vector field decays. This is valid only for 
Clw C 1. To calculate the same expression when ~ 1 one could evaluate N and 
hence 5N directly. All of this is the same as for a scalar field contribution, where the 
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calculation of N was done in Ref. [98]. 

Comparing Eq. (3.63) with (3.7) we find and to be equal to 

N l W = l^w^, (3.64) 
< = (3-65) 
Using Eq. (3.64) the isotropic part of the total power spectrum in Eq. (3.30) becomes 

Vf° (k) = Nfa (k) (l + £ , (3.66) 

and from Eq. (3.64) the preferred direction in the power spectrum in Eq. (3.29) is 

N w = W, (3.67) 

where W = W/W. 

Then the vector part of the bispectrum for equilateral configuration in Eq. (3.36) 
reduces to 

B^Oci.ka.ks) = (~) ^V + (h)V + (k 2 ){l+ P (k 1 )W?+p(k 2 )Wi+ 



+w 1 w 2 



q{ki)q(k 2 ) - ^p(h)p{k 2 ) 



+ 



-W 2 p (k 2 ) q (h)] + \q (fci) q (Afc) J + c.p. (3.68) 

In the above we used the notation W\ = W • ki etc. Because the configuration of wave 
vectors k 1; \i 2 and k3 is equilateral, with the angle between any two of them being 2n/3, 
we find ki • k2 = ki • £3 = \t 2 • £3 = —\. Eq.(3.68) simplifies further if we consider a 
scale invariant power spectrum, then the expression for /^i!" 1 becomes: 



6 fCq uil ,2^2 3 (l+ 1 2 q 2 ) + [p+l(p 2 -2q>)]W ± 



where we have taken into account that the non-Gaussianity generated during the single 
field inflation is negligible (see Eq. (2.171)). The quantity W± < 1 is the modulus of the 
projection of the unit vector W onto the plane containing the three vectors ki, k2 and 
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k3. The calculation of W± in the equilateral configuration is explained in more detail in 
the Appendix A. 

Similarly to the definition of the anisotropic power spectrum in Eq. (3.29) we may 
separate the non-linearity parameter /nl in Eq. (3.69) into the isotropic and anisotropic 
parts 

lhL = f$i{l + G-W 2 ± ), (3.70) 

where Q parametrizes the anisotropy in /nl. Comparing this equation with Eq. (3.69) 
we find that in the equilateral configuration 

f equil _ >2-p2 3 i 1 + W) (o 71 \ 

and 

Cf*** = P+ \ { f~ 2 2q2) - (3.72) 

For the squeezed configuration the bispectrum from the vector field perturbation in 
Eqs. (3.39) becomes 

x [1 + p (ki) Wl + p (k 3 ) W%] . (3.73) 

Working as in the equilateral case, we find that the non-linearity parameter for the scale 
invariant power spectra is 



6 .local _ ,2,p2 3 1+pWl 

5 /NL -^ + 27V(^ + W 2 - (3 ' 74) 



Using the parametrization of Eq. (3.70) we write for the squeezed configuration 

/nlIo = f^l - J 1 9 . (3-75) 

and 

g local = p. (3.76) 

As one can see from the above equations /nl in general depends on W± in both 
configurations, i.e. Q ec i ml ^ and Q local ^ 0. This means that /nl is anisotropic, with 
the same preferred direction as in the power spectrum (c.f. Eqs. (3.29) and (3.67)). The 
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isotropic parts of /nl in both configurations may be rewritten as 

,equil __ f local A , 1 _ 2 3 _ 1 + \^ (o 77 \ 

Given the (quasi) exponential expansion of the Universe during inflation, the value of p 
depends only on the Lagrangian of the vector field. As Eq. (3.32) suggests it relates to g 
only indirectly, through parameters determining the generation of anisotropic as well as 
isotropic parts of Q. In other words, specifying the value of p does not determine g. In 
view of this, from Eq. (3.77) we see that the amount of non-Gaussianity is correlated with 
the statistical anisotropy in the spectrum, /nl oc g 2 . If, instead, the particle production 
is isotropic (i.e. V\\ = V+ and V- = 0) Eqs. (3.18) and (3.19) give p = q = and 
therefore ^ ec i ml = g^ ocai = o. In this case /j^l"' an d become isotropic too and both 

reduce to /nl = 5/4f2iy as in the scalar curvaton scenario, where we used Eq. (3.32) 
with the assumption <C V+, i.e. that the dominant contribution to the curvature 
perturbation is due to the vector curvaton field only. 

In addition to the /nl being anisotropic, having the same preferred direction as the 
spectrum and its magnitude being correlated with the anisotropy in the spectrum g 
from Eqs. (3.71), (3.72) and (3.75), (3.76) we find more observational signatures. From 
Eq. (3.77) it is clear that for parity conserving vector fields with q = 0, the isotropic 
parts of /nl are identical, i.e. /nl^so = /nl^so- Any departure from this equality would 
indicate parity violating terms in the Lagrangian of the vector field. But the anisotropy 
in the non-linearity parameter is configuration dependent with g c i ml ^ ^l° ca l j n both - 
parity conserving as well as parity violating - theories. In the squeezed configuration, 
giocal - g s ens itive only to the anisotropy in the parity conserving perturbations of the 
vector field. But in the equilateral configuration, (J et i ul1 is also correlated with the amount 
of parity violation of the field. In both cases the anisotropy in /nl is proportional to 
the anisotropy in the particle production of the vector field, p and q. Therefore, if the 
anisotropy in particle production is of order one or larger, anisotropic parts of /nl in 
both configurations are not subdominant. At the moment, observations do not provide 
any information about the values of Q ec i ml and ^ local . However, as can be seen from 
Eqs. (3.72) and (3.76) the observational detection of Q ec i ml and C? local would allow to 
determine uniquely the values of parameters p and q, therefore, allowing to constraint 
very narrowly the possible range of conformal invariance breaking Lagrangians for the 
vector field. 
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3.3.3. Generation of £ 

To calculate the curvature perturbation generated by the vector curvaton field consider 
an era after reheating when the Universe is radiation dominated, with the energy density 
decreasing as p 7 oc a" 4 . As was discussed at the end of section 3.3.1 the relative energy 
density of the heavy vector field during this era increases as pw/p^ oc a (see Eq. (3.60)). 
When the vector field becomes dominant (or nearly dominant) it imprints its perturbation 
spectrum onto the Universe. 

The curvature perturbation £w generated by the vector field is calculated as follows. 
On the spatially flat slicing of space-time using Eq. (2.151) we can write for the vector 
field 

Am,, 

(3.78) 



dec 



3p w 

where we considered that the decay of the vector field (labeled by 'dec') occurs after the 
onset of its oscillations so that it is pressureless, as shown in Eq. (3.60). Note that, since 
Qw is determined by the fractional perturbation of the field's density, which is a scalar 
quantity, the perturbation (\y is scalar and not vector in nature. 

Since Eq. (3.47) is a linear differential equation, A and its perturbation 5 A satisfy the 
same equation of motion. Therefore, they evolve in the same way, which means that 
5 A I A remains constant, before and after the onset of oscillations. As was discussed in 
section 3.3.1 the massive vector field acts as an underdamped harmonic oscillator. The 
energy of such oscillator is determined by the amplitude of oscillations. Therefore, we 
may write pw = \ m2 ||W^|| 2 > where we used the physical vector field W = A/a defined 
in Eq. (3.6) and ||W|| is the amplitude of the oscillating physical vector field. From the 
above we obtain 



t Sp w 



3p w 



dec 



2 \\SW\\ 


2 SW 


2 5W 




3 \\W\\ 


dee " 3 W 


ose " 3 W 


end 



(3.79) 



where 'osc' denotes the onset of oscillations and 'end' denotes the time at the end of 
inflation. Therefore, from this equation we may write 

In the usual scalar curvaton scenario the curvaton field generates the total curvature 
perturbation in the Universe. This may be realized in the vector curvaton scenario as well 
if the curvature perturbation £w, generated by the vector field, is statistically isotropic, 
or if its statistical anisotropy is within the observationally allowed region (see Eq. (2.19)). 
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However, if £w is predominantly anisotropic it can only be a subdominant contribution 
to the total curvature perturbation f, while the dominant part must be generated by 
some statistically isotropic source. In analogy to the Eq. (2.204) in this case the total 
curvature perturbation in the curvaton scenario may be written as 

C = (1 - jv) c 7 + (iwCw, (3.81) 

where £ 7 is the dominant and statistically isotropic curvature perturbation which is 
present in the radiation dominated Universe before £jy is generated. In this equation 
Cl\Y m < 1 because the vector field must decay before it starts dominating (the 

dashed line in Figure 3.2). Assuming that £ 7 is generated by the scalar field (ft, the 
anisotropy in the power spectrum g from Eq. (3.32) becomes 

9 ~ ^P, (3.82) 

were we also assumed scale invariant power spectra for the scalar and vector field per- 
turbations. 

Because the isotropic part of the curvature perturbation is dominant, from Eq. (3.29) 
we may write (see (2.10)) 



C „ Apiso = g-^N W y/V + (g+p), (3.83) 



C 

were we used Eqs. (3.30) and (3.82). 

Using the definition of p in Eq. (3.19) we find that the typical amplitude of the vector 
field perturbation is 

6W « JV\\ + 2V+ = V^+(p + 3). (3.84) 
Combining the last two equations we obtain 



C«f 1/J %W^, (3.85) 
V 3 + p 

where 5W is evaluated at the vector field decay. For the vector curvaton scenario the pa- 
rameter Nyy can be found from Eq. (3.64) as N\y ~ Slw/ (2W), where W is evaluated at 
the field decay. Therefore, using Eq. (3.79) we find that the total curvature perturbation 
given in Eq. (3.81) is of order 

C ~ g- 1/2 n w (w, (3.86) 
were we have taken y/(g + p) / (3 + p) ~ O (1). 
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After the vector field decays, the curvature perturbation £ stays constant. This hap- 
pens at the time r^, , where T\y ~ h?m is the field decay rate and h is the vector field 
coupling to its decay products. Due to gravitational decay the lower bound for h is 

h>— • (3.87) 
mpi 

However, during its oscillations the vector field is subject to thermal evaporation. Were 
this to occur, all memory of the superhorizon perturbation spectrum would be erased; 
therefore, no Qw would be generated. Considering that the scattering rate of the massive 
vector boson with the thermal bath is Tgc ^ h T we can obtain a bound such that the 
condensate does not evaporate before the vector field decays, i.e. 

Tsc < Fw, (3.88) 

which is evaluated at H ~ The temperature of the Universe at the vector field 

decay is T ~ y/mpiT\y, giving r sc ~ yjmp\m. Substituting this into Eq. (3.88) and 
combining with Eq. (3.87) the range for h becomes 

/ \ 1/6 



< h < . (3.89) 

mpi V m Pi/ 

The lower bound in the above is due to decay through gravitational interactions, while 
the upper bound is relaxed if the vector field dominates the Universe before it decays. 
This happens if Tyy < -ffdomi where H^ om is the Hubble parameter when the vector field 
starts to dominate. Then the energy density of the thermal bath is exponentially smaller 
than pw and the vector field condensate does not evaporate. Thus it is enough to ensure 
that the vector field condensate does not evaporate before it dominates, i.e. T sc < H^om- 
For the dominant vector curvaton this bound can be satisfied even if the one in Eq. (3.88) 
is violated. 

Having discussed the general predictions of the vector curvaton scenario for the power 
spectrum and non-linearity parameter /nl let us turn now to the realization of this 
scenario in two concrete examples. In sections 3.4 and 3.5 we present two mechanisms 
for breaking the conformal invariance of the vector field and find under which conditions 
the field perturbation power spectrum is scale invariant. This allows us to calculate 
parameters p and q as well. Then we implement these models into the vector curvaton 
scenario and compute the parameter space for these models. 
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3.4. Non-minimally Coupled Vector Curvaton 

In Ref. [83] it was shown that a massive vector field may acquire a scale invariant per- 
turbation spectrum if its effective mass during inflation is —2H 2 . In this section we 
consider the realization of this scenario. The negative mass squared can be achieved by 
non-minimally coupling the vector field to gravity through the Ricci scalar term. The 
vector field Lagrangian for this model is written as 

C = -^F flu F^ + ^(m 2 + aR)A^, (3.90) 

where 

Ffiv = d„A u - d u A^ (3.91) 

and R is the Ricci scalar, with a being a real coupling constant. For the further discussion 
let us define the effective mass of the vector field as 

M 2 = m 2 + aR. (3.92) 

Starting from Ref. [79] this action with m = was invoked by many authors for the 
generation of primordial magnetic fields, where is identified with the electromagnetic 
field. Note, that due to the non-minimal coupling, this field is no longer gauge invariant. 
This is in contrast to the electromagnetic field in the Standard Model. But because in 
the present Universe R is very small, it is thought that at present the electromagnetic 
field is approximately gauge invariant. 

In our case we don't have to worry about the gauge invariance because we don't 
associate A^ with the electromagnetic field. Even more so, we don't assume that the 
vector field A^ couples to any scalar field through the covariant derivative of the form 

3.4.1. Equations of Motion 

During inflationary stage the spatial curvature of the Universe is inflated away. And in 
accordance with the curvaton scenario, the vector field during inflation is subdominant 
and does not influence the expansion of the Universe. Therefore, we can assume to a 
good approximation that inflationary expansion is homogeneous and isotropic with the 
flat space-time metric in Cartesian coordinates given by 

d 2 s = d 2 t - a 2 (t) d 2 x. (3.93) 
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In this case the Ricci scalar takes the form 

R = -6 + =-6(H + 2H 2 ^j = 3 (3w - 1) H 2 , (3.94) 

where u>~— 1, w = 1/3 and w = during (quasi) de Sitter inflation, radiation and 
matter dominated epochs respectively. We will further assume that inflationary ex- 
pansion is of the (quasi) de Sitter type with H ~ constant, making the Ricci scalar 
R ~ —12H 2 . With this condition, the effective mass of the vector field during inflation 
becomes M 2 ~ m 2 — 12aH 2 ~ constant. 

Calculating equations of motion for the vector field components we will mainly follow 
Ref. [83]. Using Eq. (3.90) and the variation principle 



8A. B '^ktf = ° <M5) 

we find the field equation for the vector field as 

(d^ln^Vg)] F» v + M 2 A v = 0, (3.96) 



where T> g = det(<w). With the FRW metric in Eq. (3.93) and the field equation in 
Eq. (3.96) the equation of motion for the temporal component [y = 0) is found to be 

V • A - V 2 A t + (aM) 2 A t = 0, (3.97) 

where V is the divergence and V 2 = didi is the Laplacian. In the same way we may find 
the equation of motion for the temporal component (y = i): 

A + HA - a~ 2 [V 2 A - V (V • A)] + M 2 A = V (A t + HA t ) . (3.98) 



A third useful relation is the integrability condition, which is obtained by contracting 
Eq. (3.96) with d v : 

(aM) 2 A t + 2 (aM) 2 ^A t - M 2 V • A + 3H (v 2 A t - V • A) = 0. (3.99) 
Combining the integrability condition with Eq. (3.97) we find 

A t + I 3H + 2^) A t - a~ 2 V ■ A = 0. (3.100) 
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From this equation we can see that the temporal component of the vector field is non- 
dynamical. Taking the gradient of Eq. (3.100) and plugging it into Eq. (3.98) we arrive 
at 

A + HA + M 2 A - <T 2 V 2 A = -2yH+ VA t . (3.101) 

Classical inhomogeneities of the vector field are diluted by inflation. Therefore, we can 
neglect all gradient terms 

d i A lx = V/i. (3.102) 

Using this condition in Eqs. (3.97) and (3.98) we find that for the homogeneous mode 
the temporal and spatial components of the vector field obey 

A t = 0, (3.103) 
A + HA + M 2 A = 0, (3.104) 

where we have used the choice of the coordinates such that the homogeneous vector 
field has components A^ = (^4^, 0, 0, j4). We see that a temporal component of the 
homogeneous massive vector field in the FRW Universe is zero. The equation of motion 
for the spatial component in Eq. (3.104) may be rewritten in terms of the physical vector 
field. For this model it is W = A/ a. In the (quasi) de Sitter space-time {H ~ 0) 
Eq. (3.104) becomes 

W + 3 HW + (M 2 + 2H 2 ) W = 0. (3.105) 
However, to quantize the vector field we need to perturb the field 

, A fix) = Kit) + 8A(t,x) , 
A li (t,x) = A ll (t) + SA ll (ik,x)=>{ V ' ' W V ' ; ' (3.106) 

A t (t,x) = 6A t (t,x), 

where we have used Eq. (3.103). From Eqs. (3.97) and (3.101) we find that the evolution 
of perturbations of the vector field in (quasi)de Sitter space-time [H ~ and M ~ 0) 
follow equations 

V- 5 A- V 2 5A t + {aM) 2 5A t = 0, (3.107) 
5A + H5A + M 2 5A-a" 2 V 2 5A = -2HV5A t . (3.108) 

Going to the Fourier space (see Eq. (3.10)) the first equation for the temporal component 
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becomes 



id t k • 5A k ) 
k 2 + (aM) 2 



(3.109) 



where k 2 = kk and the subscript in 5A kfl denotes the Fourier mode of the vector field 
perturbation. Using the Fourier transform of Eq. (3.108) and plugging it in Eq. (3.109) 
we find 



5A k + H5A k + M 2 5A k + 



SA k + 2H 



kd t (k • 5A k ) 
k 2 + (aM) 2 



0. 



(3.110) 



The massive vector field has three degrees of freedom and all three of them have to be 
quantized. Similarly to Eq. (3.11) we decompose 5A k into three polarizations 



5A kl = ^2et['k)^\ x . 

A 



(3.111) 



and choose ef to denote three vectors of the circular polarization in Eq. (3.12). Two trans- 
verse ones are perpendicular to the wave-vector k giving efki = 0, where '+' stands for 
the left-handed 'L' or right-handed 'R' polarizations. Substituting this into Eq. (3.110) 
we find 



a 2 + Hd t + m 2 + 



5 A+ = 0. 



(3.112) 



For the longitudinal polarization e-ki = k, and taking into account that e" = k = k/A; 
from Eq. (3.110) we find 



df+ 1 



2k y 



k 2 + (aMf 



Hd t + M + 



SA U = 0. 



(3.113) 



In the following sections we quantize the transverse and longitudinal degrees of freedom 
separately. 



3.4.2. Transverse Modes 

Let us rewrite the equation of motion of the transverse polarizations in Eq. (3.112) in 
terms of the physical vector field Wa and conformal time r = J (it /a: 



d 2 + 2-d T + (aM) 2 + k 2 + - 
a a 



w + = 0, 



(3.114) 
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where prime denotes the derivative with respect to the conformal time r and w+ (r, k) 
5A + (r, k) /a (r) is defined in Eq. (3.11). 

To find initial conditions for this field let us make a transformation in Eq. (2.81) 



r s-di 

J a 



and bring the equation of motion into the form of the harmonic oscillator 



d 2 T + (aM) 2 + k 2 



(3.115) 



(3.116) 



In the subhorizon limit, for the modes with k 2 3> |aM| 2 , Eq. (3.116) reduces to the flat 
space-time harmonic oscillator. In other words, if we write the action for functions (p+ (k) 
in the limit k/aH — > oo it would correspond to the collection of harmonic oscillators. 
Choosing the initial state for ip + to correspond to a vacuum (no particles or minimum 
energy), from Eq. (2.97) it becomes (p+ (k) = ex.p(ik/aH) /V2k, or going back to the 
physical field 



lim w + 



_^ e ik/aH 



(3.117) 



We are interested in the power spectrum of classical perturbations of the vector field 
on the superhorizon scales when k aH, which from Eq. (3.16) can be calculated as 



k 3 

aH 



(3.118) 



To find the power spectrum we need to solve the equation of motion in Eq. (3.114) in 
the limit k/aH — > with the vacuum initial conditions in Eq. (3.117). For this purpose 
it is convenient to rewrite this equation in the form 



x 



x- 



(3.119) 



where 



x 



aH 



(3.120) 

This is a Bessel equation with a general solution of the form 

w + = x 3/2 [CxHl (x) + C 2 V.l (x)] , (3.121) 
where C\, C2 are constants of integration and Ti, 2 are the Hankel functions of the 
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first and second kind respectively with 



" s vHf) = vi +i2q -(?) 2 (3 - i22) 

Taking the limit x — > oo and matching to the vacuum value in Eq. (3.117) we get 

w + = a-^e^W(?) l Hl(x). (3.123) 
For the superhorizon perturbations (x — > 0) this solution becomes 

-1 r W Q i(2i/-l)w/4 ( x \^ u 



where r (f) is the Gamma function. Plugging this result into Eq. (3.118) we find the 
power spectrum for the transverse polarization to be equal to 



vr v ' V2vr / \2aH 
Note that for the light vector field, when M — > 0, the power spectrum becomes 

Pr =(^) 2 - (3.126) 

Comparing with Eq. (3.117) we see that it is simply the vacuum value. This is in accord 
with the expectation that the massless vector field is conformally invariant and does not 
undergo particle production. 

The scale dependence of the power spectrum can be parametrized in the usual way 
as V+ oc fe nv_1 , so that n v = 1 corresponds to a flat spectrum. Comparing this with 
Eq. (3.125) we find that the spectral index is n v — 1 = 3 — 2v and the scale invariant 
spectrum of the vector field perturbation is achieved if 

3 

n v = 1 => v = - => M 2 = -2H 2 , (3.127) 

which agrees with the findings of Ref. [83]. With this condition the power spectrum 
becomes 

P + =(|) 2 , (3.128) 



the same as for the massless scalar field. 
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The condition in Eq. (3.127) is satisfied if the coupling constant of the vector field to 
gravity is 

I 1 / rn \ 2 

(3.129) 



1 

a « - 
6 



1 (vn\ 2 



2 

from which it is clear that for scale invariance we need a > 1/6. If m > H then scale 
invariance is attained only when a is tuned according to Eq. (3.129). However, if m <C H 
then scale-invariance simply requires a ~ 1/6. In the latter case m and H do not have 
to balance each other through the condition in Eq. (3.129) and can be treated as free 
parameters. We feel that this is a more natural setup. 
With a = 1/6 the v parameter in Eq. (3.122) becomes 



9 /m\ 2 

which is reminiscent of the scalar field case, where perturbations become classical in the 
superhorizon limit if v is real, corresponding to m? < 9H 2 /A (see the discussion below 
Eq. (2.113)). 

3.4.3. The Longitudinal Mode 

Let us first rewrite the equation for longitudinal perturbations of the vector field in 
Eq. (3.113) in terms of the conformal time r 



5Au = 0. (3.131) 



In the previous discussion on the perturbations of the transverse components, we found 
that the scale invariant spectrum is achieved if the effective mass squared of the field 
is negative and equal to M 2 = —2H 2 . But in this case the second term in the above 
equation becomes singular when (k/a) 2 = 2H 2 . This might indicate that the longitudinal 
vector field perturbation becomes unstable when approaching the horizon exit. But the 
two independent solutions of this equation 

5A± oc (-kr + ^±2i^j e TikT (3.132) 

show that this is not the case. 3 

M 2 = —2H 2 = constant corresponds to the flat perturbation power spectrum of 
transverse modes and, as will be seen later, of the longitudinal mode too. However, the 



3 To find this solution we have used the relation r = — (aH) 1 , which is valid in de Sitter space-time. 
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exactly flat spectrum is excluded by observations (Eq. (2.16)). Therefore, one would 
expect that in the realistic theory the condition M 2 = —2H 2 is violated by a small 
amount to give the correct spectral tilt and the Hubble parameter is not exactly constant 
during inflation. In this case the solution in Eq. (3.132) is not valid and one may be 
worried that for general effective negative mass M 2 < and H / the solution of 
Eq. (3.131) is still singular at (k/a) 2 — > |M 2 |. We can prove that this is not the case 
using the Frobenius method for differential equations with regular singular points (see 
for example Ref. [99]). 

Using Eqs. (3.100) and (3.101) with M/0 the equation of motion for the longitudinal 
mode <L4.ii becomes 



<9 2 - (3w + 1) 



k 2 aH 



k 2 + a 2 M 2 



d T + (k 2 + a 2 M 2 ) 



<L4n = 0, 



(3.133) 



where w is the barotropic parameter of the dominant component of the Universe which 
drives inflation. For de Sitter expansion w = —1 and we recover Eq. (3.131). However, 
for this calculation we do not assume de Sitter inflation and consider a constant w in the 
range — 1 < w < — |, which is necessary for the accelerated expansion of the Universe 
(see the discussion in section 1.5.2). This equation is valid for a general non-minimal 
coupling constant a defined in Eq. (3.92) and we used that aR 3> m 2 during inflation. 
Let us first we make a change of variables 



y 



a\M\ 



1, 



(3.134) 



with y varying in the region —1 < y < oo. Eq. (3.131) with this transformation translates 
into the form 



,2 1 (y + 2) \M 2 \ y 

°y 2y(y + l)° y + H 2 (3u) + l) 2 (y + l) 2 



5A n = 0, 



(3.135) 



with M 2 < and the regular singular point at y — > 0, corresponding to (k/a) 2 — > |M 2 |. 
The general solution of this equation can be found using the ansatz 



<L4, 



71=0 



(3.136) 



where Dq ^ 0. In this case the series in Eq. (3.136) is convergent at least in the region 
— 1 < y < 1 (corresponding to |M 2 | < (k/a) 2 < 2|M 2 |) without a singular point at 
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y = 0. Our aim is to prove that the solution in Eq. (3.136) is not singular even at 
(k/a) 2 = |M 2 |, i.e. y = 0. This will be the case if the power series ansatz in Eq. (3.136) 
has two independent solutions and if s + n > for all n, i.e. there are no negative powers 
of y in the series. To show this let us substitute Eq. (3.136) into Eq. (3.135) giving 



n=0 



4 (s + n) (s + n - 2) y s+n ~ 2 
+4 (s + n) I s + n - 



+ 8 (s + n) ( .s -+ ii 
3 



y 



s+n 



+ 12a 



7 
4 

3w- 1 



(3w + 1) 



s+n+1 



0, (3.137) 



where we also used Eq. (3.94). In order for this equality to be valid, coefficients in front 
of each y with the same power must vanish. The coefficient in front of the term with the 
smallest power, i.e. y s ~ 2 , is 4Dqs (s — 2). Because Dq ^ 0, from the indicial equation 
s (s — 2) = we find 

s = or s = 2. (3.138) 

Because these two solutions differ by an integer, it might be alarming that both series in 
Eq. (3.136) with s = and s = 2 do not provide two independent solutions. In this case 
the second independent solution would involve the term Iny, which indeed diverges at 
y — > 0. However, by closer inspection of Eq. (3.137) we find that the coefficient D2 of the 
series with s = is arbitrary, thus the power series in Eq. (3.136) with s = and s = 2 
do give two independent solutions. In addition they do not involve negative powers of y, 
i.e. s > 0, therefore, the solution with the ansatz in Eq. (3.136) converges at the singular 
point y = 0. This proves that during inflation, when M 2 < 0, the solution of Eq. (3.133) 
is stable when the wavelength of the perturbation approaches (k/a) 2 —¥ \M 2 \. 

Let us turn now to the quantization of the longitudinal mode in (quasi) de Sitter 
space-time. From Eq. (3.113) the equation of motion for the longitudinal physical field 
in the conformal time is 



d 2 + 2- 1 + 



A: 2 



a V k 2 + (aMY 



d T + k 2 + (aM) 2 + 2 



A; 2 



a J k 2 + (aNiy 



+ 



mi 1 = 0. 



(3.139) 

To quantize the longitudinal mode, let us write the Lagrangian corresponding to the 
equation of motion in Eq. (3.131) 



C = M 2 



6 A' (r, k) 



(A;/a) 2 + M 2 



a 2 \SA\\ (r,k)| 



(3.140) 
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This Lagrangian can also be achieved by perturbing the full Lagrangian in Eq. (3.90) and 
it is unique up to the total derivative. To set the initial conditions for the subhorizon 
modes we use the transformation 

^l= 7 - 1 «5A||, (3.141) 

where 7 is the Lorentz boost factor 



E 



^(k/af + \AP\ k/a 



1 = WC~^W ^WY (3 ' 142) 

and the last equality is taken in the limit k/a 3> \M\. With this transformation the 
Lagrangian for subhorizon modes reduces to that of the simple harmonic oscillator 



C = ± 



k 2 \<pu I , (3.143) 



where the sign ± is that of M 2 , hence negative for the case of interest M 2 ~ —2H 2 . 

The Lagrangian in Eq. (3.143) is the same as of the harmonic oscillator. Choosing 
initial conditions to correspond to the vacuum state, we have 



n = -^f ■ ( 3 - 144 ) 



This is similar as for the scalar field case, except that for M 2 < the Lagrangian has a 
negative sign. Because of the wrong sign, initial conditions in Eq. (3.144) are not identical 
to the scalar field case, since for the longitudinal mode occupied initial states would have 
negative energy density and pressure. As the pressure is negative it is not dangerous 
for inflation. Instead, it is the negative energy density that is dangerous. As the total 
energy density is required to be positive, the negative contribution of occupied states has 
to be less than the total at the beginning of inflation. This is satisfied by assuming that 
initially the occupation number is much less than 1 (as in the scalar field case), justifying 
both the choice of initial mode function and the assumption of the vacuum state. 

Matching the solution in Eq. (3.132) to the vacuum initial conditions from Eq. (3.144), 
vac = aw \\ = ~^p^ ex P ( — ikr), and using Eq. (3.118) we find that the power spectrum 
for the longitudinal mode with a = 1/6 is 

V ll = 2^y = 2V + . (3.145) 
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This corresponds to p = 1 in Eq. (3.19), meaning that the particle production of the 
vector field is anisotropic. 

3.4.4. The Stability of the Longitudinal Mode 

As shown above, the possible instability of the longitudinal mode when (k/a) 2 — > \M 2 \ is 
absent. However, in Ref. [100] it has been noted that other instabilities might be present 
for non-minimally coupled vector field. The first concern is that for m 2 >C \R\ the kinetic 
term of the longitudinal mode is negative on the subhorizon scales (see Eqs. (3.92) and 
(3.140)). As a result one might suspect that corresponding particles carry a negative 
energy and they can be created from the vacuum making it unstable. This is indeed 
the case for minimally coupled scalar field with negative kinetic term. The latter field is 
called a ghost and is cosmologically unacceptable as it would create too many photons 
from the present day vacuum (Ref. [101]). However, the flat space-time calculation of 
Ref. [101] can not be directly applied to a vector field with non-minimal coupling to 
gravity and currently no such calculation exits. Moreover, the bound on the photon 
creation from the vacuum at the present Universe is irrelevant for the vector curvaton 
scenario as its bare mass squared m 2 in Eq. (3.92) dominates, making M 2 positive. But 
even for models with negligible m 2 one wouldn't expect a large particle creation, as in 
the present day Universe \R\ ~ 10 _66 eV 2 , which is extremely small compared to other 
energy scales. 

Another concern is about the singularity when M 2 — > 0. After inflation R tx i -2 , and 
when both terms in Eq. (3.92) cancel each other out M 2 vanishes. As shown in Ref. [100] 
this results in a singularity which invalidates linear calculations around this point. To 
evaluate the effects of this, one needs to perform a full non- linear calculation which has 
not been done to the present moment. However, as the period of non-linear evolution is 
very brief, one might expect that linear calculations before and after M 2 = will match. 
Or one can assume that m 2 = and consider more complicated models to generate mass 
term for the vector field, in which case the aforementioned singularity can be avoided. 4 

3.4.5. Statistical Anisotropy and Non-Gaussianity 

Let us calculate the statistical anisotropy and non-Gaussianity for this model. Because 
p = 1, the dominant contribution to the curvature perturbation is assumed to be gener- 
ated by the scalar field. The parity conserving transverse power spectrum of the vector 
field perturbation in Eq. (3.128) and the power spectrum generated during the single 

4 A more thorough discussion of these issues can be found in Ref. [102]. 
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scalar field inflation are equal, i.e. V+ = V^. Thus the isotropic part of the curvature 
perturbation spectrum can be written as 

Vf° = V (j> N 2 (l + 0- (3.146) 
While the anisotropy parameter from Eq. (3.32) becomes 



(3.147) 



This model does not have parity violating terms, and from Eqs. (3.19) and (3.145) we 
find 

p=l and q = 0. (3.148) 

Thus, the anisotropy in the vector field is rather strong, which means that it will have 
to remain subdominant, i.e. Qyy <C 1. Using this and Eq. (3.69), the for the 

non-minimally coupled vector curvaton is found to be 



6 ,equii _ £ 2 f i , 9 

5 /nl ~ 2 n w V 8 



Similarly, /*?£ for the squeezed configuration in Eq. (3.74) is 

^/NL Cal = 2£;(l + TVi) (3.150) 

Since V+ = = = (H/2tt) 2 , for the typical values of the perturbation we 

have 8(f) ~ 5W{ ~ H. This means that, in order for the vector field contribution to 
be subdominant, we require N\y <C Nj, (c.f. Eq. (3.7)), which from Eq. (3.31) gives 
£ <C 1. Using these results and p = 1 (see Eq. (3.148)) from Eq. (3.32) we find g ~ £. 
Thus, in view of Eqs. (3.149) and (3.150), we see that /nl ~ 9 2 /^w- Therefore, we find 
that the non-Gaussianity is determined by the magnitude of the anisotropy in the power 
spectrum. 

This prediction is valid in the regime |<5W/W| <C 1 which corresponds to Q 2 ^ > 
which implies /nl ^ g 3 ^ 2 / ' \ /Pf' ^ or smaner ^W> the contribution of the vector field 



C 

perturbation to £ is of order Qyt/^VF/^VF 2 ) 1 / 2 ]. In other words, it is of order an d 
is the square of a Gaussian quantity. The resulting prediction for its contribution to /nl 
would be given by the one-loop formula which is calculated in Ref. [97]. 
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3.4.6. The Energy-Momentum Tensor 

Let us now study the evolution of the vector field. For the scale invariant perturba- 
tion spectrum with a = 1/6 from Eq. (3.104) we find the equation of motion for the 
homogeneous mode of the vector field 

W + 3HW + m 2 W = 0, (3.151) 

which is identical to the one of a massive scalar field. When m <C H it has the solution 

W = W + CJ ( - W ~ 1 \ (3.152) 

where Wo and C are constants of integration. Because the second term in Eq. (3.152) 
is decaying, as long as m <C H the physical vector field develops a condensate which 
remains constant W ~ Wq. 

We can follow the evolution of the vector field condensate by considering the energy 
momentum tensor, which can be written in the form 

= diag (pw, ~P±, -PL, ~P\\) , (3.153) 

where [88] 

p w = ^W 2 + ^m 2 W 2 (3.154) 
and the transverse and longitudinal pressures are 

P± = g (W 2 - m 2 W 2 ^ + \(^HW + HW + 3H 2 W^j W, (3.155) 
PH = -- (w 2 -m 2 W 2 ^ - -(2HW + HW + 3H 2 W^j W. 

Thus, the energy-momentum tensor for the homogeneous vector field is, in general, 
anisotropic because pu ^ p±. This is why the vector field cannot be taken to drive 
inflation, for if it did it would generate a substantial large-scale anisotropy, which would 
be in conflict with the predominant isotropy in the CMB. Therefore, we have to investi- 
gate whether, after inflation, there is a period in which the vector field becomes isotropic 
(i.e. p± w pn) and can imprint its perturbation spectrum onto the Universe without such 
problems. 

Considering the growing mode in Eq. (3.152) and Eqs. (3.154), (3.155) we see that, 
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during and after inflation, when m <C H , we have 



1 



p w ~ -m 2 WQ and p± 



(l-w)H 2 W$. 



(3.156) 



Hence, the density of the vector field remains roughly constant, while the vector field 
condensate remains anisotropic after inflation. 

The above are valid under the condition m -C H. However, after the end of inflation 
H (t) oc t , so there will be a moment when m ~ H. After this moment, due to 
Eq. (3.94), the curvature coupling becomes negligible and the vector field behaves as a 
massive minimally-coupled Abelian vector field. As shown in Eq. (3.52), when m > H 
a massive vector field undergoes (quasi) harmonic oscillations of frequency ~ m, because 
the friction term in Eq. (3.151) becomes negligible. In this case, on average over many 
oscillations W 2 ~ m 2 W 2 . Hence, Eqs. (3.154) and Eq. (3.155) become 



Pw 
P± 



~ mW 2 , 



(3.157) 



-mH 



3 , . (E 
1 + - (1-w — 

4 \m 



"2^ 3 

The effective barotropic parameters of the vector field are 



W 2 . 



< w± 



"'II 



3 , / H 
1 + - (1-w - 

4 \m 



H 



<1, 



(3.158) 



where w± = p±/pw anci w \\ = P\\/pw- By virtue of the condition m 3> H, we see 
that, after the onset of the oscillations, w±,W\\ — > 0. This means that the oscillating 
massive vector field behaves as a pressureless isotropic matter, which can dominate the 
Universe without generating an excessive large-scale anisotropy. Moreover, as was shown 
in Eq. (3.60) the energy density decreases as pw oc a -3 , i.e. like dust. Thus, if the 
Universe is radiation dominated, pw/p °t a while oscillations occur, so the field has a 
chance to dominate the Universe and imprint its curvature perturbation according to the 
curvaton scenario. 



3.4.7. Curvaton Physics 

As we have seen in section 3.4.5 for the non-minimally coupled vector field with the La- 
grangian in Eq. (3.90) the particle production is anisotropic, and the curvature perturba- 
tion generated by such a field is statistically anisotropic. Therefore, the non-minimally 
coupled vector curvaton may generate only the subdominant contribution to the total 
curvature perturbation £, while the dominant part must be produced by a statistically 
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isotropic source. In such scenario the total curvature perturbation with statistically 
anisotropic contribution fpp was calculated in Eq. (3.86). 

If m <C H during inflation the physical vector field (being non-conformally invariant) 
undergoes particle production and obtains an approximately flat superhorizon spectrum 
of perturbations, as shown. Indeed, if in Eq. (3.130) v ~ | from Eq. (3.145) we find that 
the typical value of the vector field perturbation is (see Eq. (2.10)) 



SW end « yJV\\ + 2V+^^, (3.159) 

where 'end' denotes the typical value of the vector field perturbation at the end of in- 
flation and H* is the Hubble parameter during inflation. The curvature perturbation 
generated by the vector curvaton field was calculated in Eq. (3.80). Using this equation 
and considering that W as constant during inflation (see Eq. (3.152)), i.e. Wend = Wo, 
we may write 

Cw ~ ^. (3.160) 
Thus, from this result and Eq. (3.86) we obtain 

C~<T 1/2 1V~. (3.161) 

Wo 

At the onset of vector field oscillations the density parameter of the vector field is 

Pw ( Wo \ 2 



n w = ~ \ — \ , (3.162) 
P ympij 

where we have used the flat Friedman equation (1.16) p = 3m Pl H 2 . To avoid excessive 
large scale anisotropy the density of the vector field must be subdominant before the 
onset of oscillations, which means that Wo < mpy. 

Let us assume that inflation is driven by some inflaton field, which after inflation ends, 
oscillates around its VEV until its decay into a thermal bath of relativistic particles at 
reheating. In this scenario the Universe is matter dominated (by inflaton particles) until 
reheating. Using the above findings we can estimate the density ratio of the vector field 
at decay 

^~(=^V /2 (^y, (3.63) 

where Ty/ is the vector field decay rate and V is the decay rate of the inflaton field. If 
inflation gives away directly to a thermal bath of particles then we have prompt reheating 
and r — > H*, where H* is the Hubble scale of inflation. 
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Using Eqs. (3.161) and (3.163) and considering that at the vector field decay f^dec = 
we get 

" ( (jl-X< 2 (*^y\ (3 , 64) 



mpi V^dcc/ \min{m,r} 
The Hot Big Bang has to begin before nucleosynthesis (which occurs at the temperature 
Tbbn ~ lMeV). Hence, Ty/ > Tg BN /mpi. Using this and also min{m,r} < we 
obtain the bound 

/ \ 2/5 

H * ~ C " /5 (Odecj ( T BBN-Pl) 1/5 ^ £ 2 1/5 1012 GeV ' ( 3 - 165 ) 

where we used that f^dec i$ 1 an d C = 4.8 x 10~ 5 from COBE observations. For g < 0.3 
this is similar to the case of a scalar field curvaton [103]. 5 

Another bound on the inflation scale is obtained by considering that Tyy ~ h 2 m, 
where h is the vector field coupling to its decay products, for which h > m/mp\ due to 
gravitational decay. Thus, Ty/ > m 3 /mp^ Combining with Eq. (3.164) we obtain the 
bound ^ 

H * ~ c (?w) (mpim)1/2 ^ y * 1/4 ~ 51/4 1011 GeV ' (3 ' 166) 

where we took 1 TeV < m < T. 

Finally, an upper bound on inflation scale can be obtained by combining Eq. (3.161) 
with the requirement Wo < mpi, thereby finding 

H* < 5 1/2 C^dcc m Pi < 9 l/i 10 16 GeV, (3.167) 

where we considered that f^dcc ^ 10~ 3 , in order to avoid excessive non-Gaussianity in the 
CMB. This bound on f^dec may be found considering that /nl ~ <? 2 /^dcc ( see Eq. (3.77)) 
and the observational constraints on |/nl| ^$ 100 (Eq. (2.22)). 

As was discussed in section 3.3.3 we also need to consider the hazardous possibility of 
the thermal evaporation of the vector field condensate. If it evaporates all memory of 
the superhorizon perturbation spectrum is erased and no Qyr is generated. This puts a 
bound on the allowed values of the vector field coupling constant h to its decay products 
which was calculated in Eq. (3.89). 

The above lower bounds on can be substantially relaxed by employing the so- 
called mass increment mechanism according to which, the vector field obtains its bare 
mass at a phase transition (denoted by 'pt') with m/H pt 3> 1. The mechanism was 



5 The cosmological scales re-enter the horizon at temperatures T < IkeV, i.e. much later than nucle- 
osynthesis and well after our vector field condensate decays restoring local Lorentz invariance. 
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firstly introduced for the scalar curvaton in Ref. [104] and has been already implemented 
in the vector curvaton case in Ref. [84]. 

Let us consider now the case when a ^ g. If a = O (1) then, according to Eq. (3.129), 
a scale invariant spectrum is possible only if m ~ H*. Hence, the oscillations begin 
immediately after the end of inflation. With this in mind the previous analysis remains 
valid. In particular, the bound in Eq. (3.165) remains the same. However, the bound in 
Eq. (3.166) becomes much more stringent: 

H* > g 10 10 GeV F* 1/4 > g 1/2 10 14 GeV. (3.168) 

3.4.8. A Concrete Example 

To illustrate our findings let us consider a specific example. Let us choose a ~ g , 
m ~ lOTeV and also T\y ~ 10 _10 GeV such that the temperature at the vector field 
decay is Td ec ~ lOTeV. Such a particle may be potentially observable in the LHC. 
These values suggest h ~ 10 -7 , which lies comfortably within the range in Eq. (3.89). 
For the decay rate of the inflaton let us chose T ~ 10~ 2 GeV so that the reheating 
temperature satisfies the gravitino overproduction constraint T re h ~ ~ 10 8 GeV. 

Then Eq. (3.164) reduces to H*/m vl ~ 10" 4 C (5/^dec) 1/2 - Using this and Eq. (3.161) we 
get Wo/m-p\ ~ 10~ 4 -v/^dec- Hence, with the maximum observationally allowed statistical 
anisotropy g ~ 0.1 the lowest value for the inflationary Hubble scale is > 10 9 GeV. 

3.4.9. Summary of the RA 2 Model 

In Ref. [83] it was demonstrated for the first time that the vector field may influence or 
generate the curvature perturbation in the Universe. It was shown that a massive vector 
field may act as a curvaton field without producing excessive large scale anisotropy. In 
this reference it was also calculated that the perturbation spectrum of a massive Abelian 
vector field is scale invariant if the mass of the field is equal to M 2 = —2H 2 . Section 3.4 of 
this thesis explored the possibility of realizing the negative mass squared by non-minimal 
coupling of the vector field to gravity through the term aRA^A^, where R is the Ricci 
scalar and e is the non-minimal coupling constant. We have calculated the vector field 
perturbation spectrum for the transverse and longitudinal degrees of freedom and found 
that they are scale invariant if a = 1/6. However, the magnitude of the longitudinal 
power spectrum is twice the transverse ones, indicating that the particle production 
of the vector field is anisotropic. If such a vector field generated the total curvature 
perturbation in the Universe, the resulting magnitude of statistical anisotropy in Q would 
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violate observational bounds obtained from CMB measurements. Therefore, the vector 
curvaton considered in this section may generate only a subdominant contribution to £. 

We have also explored the parameter space of the proposed scenario. In this thesis 
calculations of the constraints for the non-minimally coupled vector curvaton model, with 
the statistical anisotropy taken into account, were performed for the first time. We have 
shown that there is an ample parameter space for the model to work by considering all 
relevant constraints in the cosmology. 

Some of recently raised concerns [100, 105] about the stability of the model were 
also addressed. It was shown that although the longitudinal mode is a ghost when it 
is subhorizon, but it may not be dangerous during inflation if we assume no-particle 
(vacuum) initial conditions (as in the scalar field case) and negligible coupling to other 
fields. It was also emphasized that the equation of motion of the longitudinal mode has 
a singular point at (k/a) 2 = |M 2 |, which might indicate that the longitudinal mode 
becomes singular at horizon exit (|M 2 | ~ H 2 ). We have obtained an exact solution for 
non-zero bare mass of the vector field, i.e. m ^ 0, and demonstrated that it is well 
behaved at all time during inflation. However, we have not addressed the instability 
of the longitudinal mode when the effective mass of the vector field becomes zero after 
inflation, i.e. when M — > 0. 

3.5. Vector Curvaton with a Time Varying Kinetic Function 

In this section we consider a vector curvaton scenario with the vector field Lagrangian 
during inflation 

C = -\fF flu F^ + ^m 2 A fl A», (3.169) 

where / = / (t) is the kinetic function and m = m(t) is the mass and both are functions 
of the cosmic time t. = d fJj A u — d v A^ is the field strength tensor. If / is time- 

independent it can be set equal to 1 because any constant value can be absorbed into 
An,. Otherwise, / represents a time-dependent coupling. 

The above Lagrangian density can be of a massive Abelian gauge field, in which case 
/ is the gauge kinetic function. However, we need not restrict ourselves to gauge fields 
only. If no gauge symmetry is considered the argument in support of the above Maxwell 
type kinetic term is that it is one of the few (three) choices [106] which avoids introducing 
instabilities, such as ghosts [105]. 
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3.5.1. Equations of Motion 

We focus, at first, on a period of cosmic inflation, during which we assume that the 
contribution of the vector field to the energy budget of the Universe is negligible. Thus, 
we take the inflationary expansion to be isotropic. As in the previous model in section 3.4 
we also assume that inflation is of (quasi)de Sitter type, i.e. the Hubble parameter is 
H ~ constant. 

Inflation is expected to homogenize the vector field. Following the analogous calcula- 
tions as in section 3.4 and Ref. [84], we find that the temporal component of the homo- 
geneous vector field has to be zero, while the spatial components satisfy the equation of 
motion 

A+ | H+ I I A+ "^A = 0, (3.170) 




where the dot denotes derivative with respect to t. From the above it is evident that the 
effective mass of the vector field is 



M 



m 

77' 



(3.171) 



where we assumed that m, f > 0. 

We perturb the vector field according to Eq. (3.106) and going to the Fourier space we 
calculate equations of motions for the transverse and longitudinal polarizations as 



d 



H + 



f 



2H + 2- 



m 



m 




a + T + W 



5A A 



0, (3.172) 
0, (3.173) 



where 5A+ and <L4|| are defined in Eq. (3.111). 

To continue we need to employ the physical (in contrast to comoving), canonically 
normalized vector field 



W 



/- 



(3.174) 



Note that the definition of W differs from the one in Eq. (3.6) because in this section W 
gets an additional factor of \f~f due to canonical normalization. 
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Expressing Eqs. (3.172) and (3.173) in terms of the physical vector field we obtain 



df + ma* + 



Wf) -7-7 H + iH2 



+ M 2 + 



(3.175) 



and 



m + \2H + 2 



M 



(k/aY 



mi (k/ a y + M 2 



dt + 



1 



+ H 



2/ 



2H + 2 



M 



2 

k\2 



Wi) - L r ! t H+iH2 



M ) (k) 2 + M 2 



+ M 2 + 



w\\ = 0. (3.176) 



Because the theory is parity conserving the Fourier mode w+ of (t, x) perturbations 
denotes both polarizations: the left-handed and right-handed, i.e. w + = yff5A + / a. 

Let us use the following ansatz for the time dependence of the kinetic function and the 
mass 

/ oc a a and m oc a 13 , (3.177) 

where a and f3 are real constants. We will also assume that / — > 1 at the end of inflation 
so that, after inflation, the vector field is canonically normalized. Then Eqs. (3.175) and 
(3.176) become 



w + + 3Hw + + 



l( a + 4)( a - 2)H 2 + M 2 + 



w + = 



(3.178) 



and 



+ 



.. „ , 2-a + 2(3 \ . 
"'ii + I 3 H 1 + r 2 — I Hw\\ + 

2 



2 ,2- Q )^ + 4+ i+r2 



2-c + V iB , + n (1 + r2) 



0, (3.179) 



where r is defined as 



M 
k/a 



(3.180) 
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3.5.2. The Power Spectrum 

To calculate the power spectrum one can proceed as in section 3.4: calculate general 
solutions of Eqs. (3.178) and (3.179), determine integration constants by matching the 
solution to the vacuum at the subhorizon limit, k/a 3> H, and calculating the field 
amplitude at the superhorizon regime when k/a <C H. However, it is difficult to find 
general solutions for these equations, therefore one needs to use approximate methods. 

In Appendix B it is shown that, in analogy to the equation of motion of a scalar 
field during quasi de Sitter inflation and with initial conditions in Eq. (3.117), the scale 
invariant perturbation spectrum for transverse polarizations in Eq. (3.178) is achieved if 

or = -l±3 (3.181) 

(i.e. either / oc a 2 or / oc a -4 ) and 

M* < H, (3.182) 

where the star denotes the time when cosmological scales exit the horizon. The latter 
condition simply requires that the physical vector field Wu is effectively massless at that 
time. 6 

For the longitudinal polarization the initial condition reads 

lim W[l =^^/ aH , (3.183) 



all 



where the Lorentz boost factor is 



In the subhorizon limit r *C 1. After finding the solution of Eq. (3.179) with this 
condition and matching it to the vacuum solution in Eq. (3.183), one can calculate the 
power spectrum of wii in the superhorizon limit. In Ref. [107] it was shown that the 
spectrum is scale invariant if 

= -i(3±5). (3.185) 
As explained in the Appendix B the value f3 = —4 must be disregarded because it implies 



6 Note that this is not the same as having being effectively massless. In the latter case the vector 
field is approximately conformally invariant and does not undergo particle production. However, the 
conformal invariance of the massless physical vector field W u is broken. 
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the massive physical vector field in the subhorizon limit. This contradicts the requirement 
for the scale invariant perturbation spectrum of the transverse modes. 

Having evaluated a and /3 in Eqs. (3.181) and (3.185) to give the scale invariant per- 
turbation spectrum of the transverse and longitudinal modes one can analyze equations 
of motion in Eqs. (3.178) and (3.179) in more detail. In Ref. [107] these equations were 
solved in different approximation regimes as well as solved numerically. Below we provide 
the summary of results. 

Case: / oc a -4 & rn oc a 

For a = —4 and (3 = 1, the equation of motion for the transverse mode functions in 
Eq. (3.178) become 

/ k\ 2 

w + + 3Hw + + f-J (l + r 2 )w+ = 0. (3.186) 

When the kinetic function of the vector field scales as / oc a~ 4 , from Eqs. (3.180) and 
(3.171) we find that r oc a 4 . In addition we assume that cosmological scales exit the 
horizon when the vector field is light. Therefore, for subhorizon perturbations when 
x > 1, where x was defined in Eq. (3.120) as x = kj (aH), the r parameter is very small, 
i.e. r <C 1. When the mode leaves the horizon x < 1 and for cosmological scales r < 1. 
However, because r is a growing function, at some later time it may become large, i.e. 
r > 1. Assuming the Bunch-Davies vacuum initial conditions, when i> 1, the solution 
of Eq. (3.186) in these three different regimes are given by [107] 



w + = a 



~ V2 \[m^i 2{x) - iJ -*i 2{x) ~\ 



w + 



1 



V2k V k 
where z is defined as 



. 2k 

H\ I ZTT 

2 



1 + tt£ 



H 

2k V k 



for x > 1 , 



for xCK 



-J-i/ 2 {z)+iz 1 J 1 / 2 (z) 



for 



<1, 



M 
3H 



(3.187) 
(3.188) 
(3.189) 

(3.190) 



and r = z/ (3x). The solution of Eq. (3.186) was calculated using numerical methods as 
well and it was found that they agree with Eqs. (3.187)-(3.189) remarkably well. 

The equation of motion for the longitudinal component with the same scaling of / and 
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m is 



w\\ + 3 + 



1 + r z 



Hw\\ + 



24 



1 + r : 



r# 2 + 



W\\ 



0. 



(3.191) 



And the solution of this equation in the same three regimes was found to be 



Wii 



6 



l Q 9l2 \lh X 2Z 1 [ J 5/2 (^) - ^-5/2 (3 



1 



i/'ll 



if 

2fc V k 



H 



Z 1 J-l/2 (z) ~ (z) 



for x > 1 , 



(3.192) 



1 



2A: V fc 



if 



-l 



1 



for x < 1< - , (3.193) 



for 



1 



<1, 



(3.194) 



which agrees with the numerical solution of Eq. (3.191) very well too. 

As is seen from Eqs. (3.188), (3.189) and (3.193), (3.194) on the superhorizon scales 
modes w+ and w\\ evolves differently if the vector field is light, M < H , or heavy, M > H. 
When the field is light, u> + is constant and W\\ oc a -1 . Therefore using Eq. (3.118) we 
find 



H\ 2 , „ 1 { H\ 2 



2vr J 



and V\\ 



oc a 



for M < H. 



z 2 \2-k ) 

Thus the typical value of the vector field perturbation is (see Eq. (2.10)) 

3H H 



5W 



M 2n 



oc a 



-3 



(3.195) 



(3.196) 



where 5W = |<5W| and we used V\\ 3> V+- 

On the other hand, when the mass of the vector field becomes comparable with the 
inflationary Hubble parameter, i.e. M ~ H, the transverse and longitudinal mode func- 
tions on the superhorizon scales, with x <C 1, become 



w+ 

W\\ = 



i (H\ /W J i/2( z ) 

u ; V 2 

L (K.\ /W J -1/2( Z ) 
2fc U ; V 2 



\W-i 



Wll 



z>l- 



(3.197) 



When the vector field becomes heavy M S> H from Eqs. (3.189) and (3.194) we find 

H\ sin (z) 



w+ 



it'u 



2k \k / z 
1 f H\ cos (z) 



2k V k 



(3.198) 
(3.199) 
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i.e. they oscillate with the same amplitude, ||iu+|| = | \w\\ 1 1, but with the phase difference 
of it/2. The frequency of oscillations is much larger than the Hubble parameter because 
z> 1, therefore it makes sense to use the average values of the power spectra over many 
oscillations. Using Eq. (3.118) we find 

^ = n = ^{^) for M > H. (3.200) 
Thus, the typical value for the vector field perturbations in this regime is 

W^^-^- — oca" 3 , (3.201) 
y/2 M 2vr v ' 

where for the scale invariant perturbations |<5W| ~ ( see Eq. (2-10)). Thus from 

Eq. (3.196) we see that the typical value of the vector field perturbation is roughly the 
same if the field is light or heavy. 

Case: / oc a 2 & m oc a 

When the vector field kinetic function / is increasing with time, i.e. a = 2, and 
/3 = 1 the effective mass of the field is constant, M = constant. The requirement that 
the field is effectively massless when cosmological scales exit the horizon in Eq. (3.182) 
suggests that M/H <C 1 at all times when the scaling above holds. Using this condition 
and scaling we can calculate the power spectra for all components of the superhorizon 
vector field perturbations generated by the particle production process. 
The equation of motion for the transverse mode functions 

w+ + 3Hw + + f - J (1 + r 2 )w + = 0. (3.202) 

is the same as for the a = —4 case except that now M = constant. This condition 
simplifies Eq. (3.202), making it possible to obtain the exact solution. After matching this 
solution to the initial Bunch-Davies vacuum state, the power spectrum on superhorizon 
scales becomes [107] 

P + =(|) 2 . (3.203) 

The equation of motion for the longitudinal component from Eq. (3.179) and a = —4 
becomes 



^ (l + r 2 ) 
a ' 



wu = 0. (3.204) 
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Again, it is impossible to find an exact solution of this equation. But using the vacuum 
initial conditions in Eq. (3.183) and solving it in two regimes, r< 1 and r > 1, and 
matching those solutions at r = 1 we find 

for x > 1, (3.205) 

for x < 1< -. (3.206) 

z 

From Eq. (3.206) we calculate the power spectrum 

^ -»(£)'(£)'■ <3 - 2o7) 

same as in Eq. (3.195). Since in this case, we have M/H = constant <C 1, the longitudinal 
power spectrum is constant, in contrast to the a = —4 case. Also, we see that V\\ 3> V+. 

3.5.3. Statistical Anisotropy and Non-Gaussianity 

The theory studied in this section has two clear advantages. First we can obtain a com- 
pletely isotropic perturbation spectrum for the vector field, which has previously never 
been achieved. As we discuss below, this means that we may consider vector fields as 
dominating the total energy density of the Universe when the curvature perturbation is 
formed. The second advantage is that we can also account for a small amount of statisti- 
cal anisotropy in the curvature perturbation spectrum depending on when inflation ends, 
again by considering the vector field alone. We also demonstrate this in what follows. 
Finally, statistical anisotropy can also be present in a correlated manner in the bispec- 
trum as well, which characterizes the non-Gaussian features of the CMB temperature 
perturbations. In view of the forthcoming observations of the recently launched Planck 
satellite mission this is a particularly promising and timely result. 

Let us first consider the case with a = —4. As we have seen in the previous section 
in this case the effective mass of the vector field is time dependent during inflation, 
M (x a 3 . When M < H the evolution of the vector field perturbations follows the power 
law on superhorizon scales (see Eqs. (3.188) and (3.193)) and the power spectra for the 
transverse and longitudinal modes are given in Eq. (3.195). Using the definition of p in 
Eq. (3.19) we find that in this regime the anisotropy in the particle production is equal 
to 

p = z~ 2 - 1, (3.208) 
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Figure 3.3.: If the vector field is light at the end of inflation, the particle production will 
be highly anisotropic. If it is heavy, the transverse and longitudinal power 
spectra are equal and p = 0, i.e particle production is isotropic. In the regime 
where rh ~ H the anisotropy in the particle production might be non-zero 
but very small, p < 1. 



where z < 1. 

In the opposite regime, when M > H, the vector field perturbations are oscillating 
(see Eqs. (3.198) and (3.199)) and the average power spectra in Eq. (3.200) are equal 
giving 

p = 0. (3.209) 

At the end of inflation the kinetic function and the mass of the vector field are stabi- 
lized: / en( j = 1 and M en( j = constant = rh. At this epoch the vector field perturbation 
power spectra become constant. Therefore, although in the curvaton scenario £ becomes 
constant only after the curvaton decay, it is enough to evaluate amplitudes of the vector 
field perturbations at the end of inflation. Thus, the value of p is frozen at the end of 
inflation and it depends on the ratio rh/H (see Figure 3.3). If the vector field is light at 
the end of inflation (or if its mass is of the order of inflationary Hubble parameter), the 
particle production is anisotropic and given in Eq. (3.208) with z = z cn d = rh/3H. If the 
field is heavy, the transverse and longitudinal power spectra are equal, i.e. the particle 
production is isotropic giving Eq. (3.209). 

If rh > H, the particle production of the vector field is isotropic, p = 0, and the vector 
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field generated curvature perturbation is statistically isotropic 

g = 0. (3.210) 

Therefore, if the vector field is heavy at the end of inflation the generated curvature 
perturbation is indistinguishable from the scalar field case. Indeed, if we plug Eq. (3.209) 
into the expression of the /nl in Eqs. (3.69) and (3.74) we find 

6 »equil _ 6 „i oca i _ 3 . , 

exactly the same as in the scalar curvaton case. In this expression we considered that 
the only contribution to Q comes from the vector field, i.e. = 0. 

If, on the other hand, the vector field mass at the end of inflation is rh ~ H , from 
Eq. (3.208) we find < p < 1. Using this and Eqs. (3.72) and (3.76) the anisotropy in 
/nl becomes 

gequil _ glocal ~ p < ^ ( 3 . 2 12) 

where again we have considered that only the vector field generates the curvature per- 
turbation. This regime is possible because the observational bound on anisotropy in the 
spectrum (defined in Eq. (3.32)) is not violated 

g&p <0.3, (3.213) 

where g < 0.3 is the observational constraint from CMB on the statistical anisotropy 
(see the discussion above Eq. (2.19)). Using Eq. (3.213) and requiring that z en( j < 1 from 
Eq. (3.208) we find the allowed range of the mass values rh for this case 

0.3 < — < 0.4. (3.214) 
m 

Unfortunately this range is very narrow and initial conditions must be tuned accurately 
to achieve this possibility. 

In the vector curvaton model with the light vector field at the end of inflation the 
dominant part of Q must be generated by the scalar field, while the vector field can gen- 
erate only a subdominant contribution. This is because if rh <C H, the anisotropy in 
the vector field particle production is large, i.e. p> 1, and this would violate observa- 
tional constraints on g. Assuming a light scalar field with perturbation power spectrum 
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= (H/2it) 2 = V+ from Eq. (3.32) we find 

9 = T^P-T^(^r) ' 3 ' 215 
1 + 4 l + £V m / 

where in the last expression we have used zCl Similarly from Eq. (3.77) we find that 
the isotropic part of /nl is 

A% = /NtL = 9 ^(f) 4 - (3-216) 

While amplitudes of the angular modulation of /nl in the equilateral and squeezed 
configurations with zCl are 




(3.217) 



and 

which are much larger than one in both cases. 

The values of the non-linearity parameter /nl calculated in this section correspond to 
the scaling of the kinetic function with a = —4. But in the limit z <C 1 Eqs. (3.215)- 
(3.217) are also applicable for a = 2. In the latter case z = constant <C 1 therefore, if 
/ oc a 2 the vector field may only generate a subdominant contribution to the curvature 
perturbation without violating observational bounds on statistical anisotropy, where the 
dominant part is produced by a scalar field. But for / oc a -4 , as we have seen in 
Eq. (3.210), the vector field can also produce the total curvature perturbation in the 
Universe. If it is heavy at the end of inflation, i.e. m 3> H, the generated £ will be 
statistically isotropic and indistinguishable from the scalar curvaton case. If, on the 
other hand, rh ~ H, it may still generate the total £ which is approximately statistically 
anisotropic within the observational bounds. 

3.5.4. Evolution of the Zero Mode 

In order to calculate the curvature perturbation associated with the vector field one needs 
to study also the evolution of the homogeneous zero mode W. Combining Eqs. (3.170) 
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and (3.174) and using Eq. (3.177), we obtain 

n • i 



W + 3HW + 



l--a\H--(a + A)(a- 2)H 2 + M 2 



W = 0, (3.219) 



where we also used the definition of the effective mass in Eq. (3.171). 
3.5.4.1. During Inflation 

As shown in Appendix B, to obtain a scale invariant spectrum for the transverse compo- 
nents of the vector field perturbation we require /(a) to scale according to Eq. (3.181), 
i.e. q = — 1 ± 3. Using this and considering the (quasi)de Sitter inflation (with H ~ 0) 
the above becomes 

W + 3tfW + M 2 W = 0. (3.220) 

We show below that, when M <C H (true at early times when a = —4; always true when 
a = 2), the solution of the above is well approximated by 

W ~ Ci + C 2 a- 3 , (3.221) 

where C{ are constants. The dominant term to the solution of Eq. (3.221) is determined 
by the initial conditions. We choose initial conditions for the vector field zero-mode 
based on energy equipartition grounds. As is demonstrated in what follows, if the energy 
equipartition is assumed at the onset of inflation, the dominant term turns out to be the 
decaying mode W oc a -3 when a = —4, and the "growing" mode W = constant when 
a = 2. 

To apply energy equipartition in the initial conditions we need to consider the energy- 
momentum tensor for this theory, which, from Eq. (3.169) is given by [84] 

1> = / (jg^FprF*" - F w fA + m 2 (a^A„ - l -g^A p A^ . (3.222) 

If we assume that the homogenized vector field lies along the z-direction, we can write 
the above as [84] 

T ^ = di&g(p w , ~P±, ~P±,+P±) , (3.223) 

where 

PW = Pkin + Vw , P± = Pkin - Vw , (3.224) 
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with 



Pkii 



Vw 



-fF^F^ = -a- 2 fA 2 = - 
4 2 2 



1 2 



W - - (a - 2) flW 



-m 2 A^ = -a" 2 m 2 ^ 2 = -M 2 W 2 , 
2^2 2 



(3.225) 



(3.226) 



where A = |A|, we used Eqs. (3.174) and (3.177), and we assumed a negative signature 
for the metric. 

Energy equipartition corresponds to 



(Pkin)o — (Vw)o 



(3.227) 



where the subscript '0' indicates the values at some initial time, e.g. near the onset of 
inflation. 

Case: / oc a~ 4 

In this case M oc a 3 and the solution to Eq. (3.220) is 



W = a 



C 3 sin cos (§) 



(3.228) 



When M > H the above shows that the amplitude of the oscillating zero mode is 
decreasing as ||W|| oc a -3 . In the opposite regime, when M <C H the solution above is 
well approximated by Eq. (3.221) with C\ = C30,q 3 Mo/3H, where we considered that 
a~ 3 M = a^" 3 Mo = constant. Using this, the constants Ci and C3 in Eq. (3.228) can be 
expressed in terms of initial values of the field amplitude Wq and it's velocity Wq: 



C 2 



Wo 3 1 A 
— a and C 3 



W + 3HW 
Mo 



a 3 , 



(3.229) 



Assuming initial equipartition of energy we can relate Wq with Wq. From Eqs. (3.225) 
and (3.226), setting a = —4, we readily obtain 



Pkin = \( W + 3HW) 2 and V w = ^M 2 W 2 . 



(3.230) 



Then, using Eq. (3.227), we get 



Wo ~ Wo(-3H±M ) 



(3.231) 
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Substituting this relation into Eq. (3.229) we find that the evolution of the vector field 
W in Eq. (3.228) takes the simple form: 



W = W, 



do 



-3 



^ cos (an ± I 



(3.232) 



Note that this equation is valid for any value of M. However, we can see that when 
M <C H the zero mode of the vector field is decreasing as W oc a -3 , but when M S> H 
it oscillates rapidly with a decreasing amplitude proportional to a~ 3 . On this basis we 
can assume that the typical value of the zero mode during inflation always scales as 



W oca' 3 . 



(3.233) 



With the assumption of initial equipartition of energy for the vector field at the on- 
set of inflation we can calculate the kinetic and potential energy densities. 7 Inserting 
Eq. (3.232) and its derivative into Eqs. (3.225) and (3.226) we find 



Pki: 



M 7T 

WW sin( 



and V\y 

Hence, the total energy density is constant 



Pw = M 2 Wo 2 . 



M TT 

W M cos [ 



(3.234) 



(3.235) 



Because this relation is independent of the vector field mass M it is valid in both regimes: 
when M <C H and W follows a power law evolution, and when M 3> H and W oscillates. 
This is valid as long as f(a) and m(a) are varying with time. 

In the vector curvaton scenario the vector field must be subdominant during inflation. 
From Eq. (3.235) we see that assuming this to be the case at the onset of inflation, it 
will stay so until the end of inflation irrespective if the field is light or heavy. 

Case: / oc a 2 

In this case, M = constant, which means that the solution of Eq. (3.220) is 



W = a 



-3/2 



C 



\aS* {h > +C 2 a V4 U) 



(3.236) 



Since in this case M <S H, the above solution is always well approximated by Eq. (3.221) 
and there is no oscillating regime. 



r By "potential" we refer to the energy density stored in the mass-term V~w = —km^A^A 11 . 
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Now, Eqs. (3.225) and (3.226) take the form 

Pkin = ^W 2 and V w = ^M 2 W 2 . (3.237) 
Combining Eqs. (3.221) and (3.237), we find 

Pkin = \n 2 C 2 a'\ (3.238) 
Thus, at the onset of inflation assuming energy equipartition in Eq, (3.227) gives 

where we used that Mq = M <C H . Inserting the above into Eq. (3.221) we find 
W = a 3 C 2 

because, after the onset of inflation, (ao/a) 3 <1< 3H/Mq. 

Therefore, we have found that W remains constant. Since M = constant, this means 
that Vw also remains constant. On the other hand, Eq. (3.238) suggests that p\^ in oc a -6 . 
Thus, since we assumed energy equipartition at the onset of inflation, we find that, during 
inflation, pkin *C Vw- Hence, 

pw^Vw — M 2 Wl (3.241) 

where M = constant = Mq. This result is the same as in the case / oc a -4 in Eq. (3.235) 
and the vector curvaton field is ensured to be subdominant during inflation (as required 
by the curvaton mechanism) if it is subdominant at the onset of inflation. 

3.5.4.2. After Inflation 

At the end of inflation we assume that the scaling of / and m has ended and we have 

/ = 1 and m = m. (3.242) 

Hence, Eqs. (3.235) and (3.241) no longer apply. The evolution of pw is determined as 
follows. 

As mentioned already, after the end of scaling, a = and M = rh. Then, Eqs. (3.225) 



a) M 



constant ~ Wq 



(3.240) 
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and (3.226) become 



Pkin = \(W + HWf and V w = ^rh 2 W 2 . 



(3.243) 



The behavior of p^ n and Vw depends on whether the vector field is light or not. To 
see this let us calculate the evolution of the field after inflation. With the conditions 
in Eq. (3.242) the physical vector field of Eq. (3.174) is W = A/a, while Eq. (3.219) 
becomes 

W + 3HW +(H + 2H 2 + rh 2 ) W = , (3.244) 



where the Hubble parameter after inflation decreases as H(t) = 3 ( 1 _^ w ) f ) with w = p/p 
being the barotropic parameter of the Universe. Solving Eq. (3.244) we find 



IT" 

W + HW 



1 w — l 
t 2 W+ 1 



1 w-l 

rh t 2 w + 1 



C\J V (rht) + C2J-V {pit) 
C\J v -\ (rht) - C2J1-V {rht) 



(3.245) 
(3.246) 



where v 



l+3w 



6(i+ w ) • ^ ne can easu y see that the vector field behaves differently if it is 
light, rht <C 1, or heavy, rht S> 1. 

Let us first see what happens if the vector field is light. Then, Eqs. (3.245) and (3.246) 
can be approximated as 
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W + HW 



1 w-l 
£2 w+1 



Ci 



rht 



1 w-l 

rht 2 w + 1 



T(l + v) \ 2 



T{l + v) V 2 J 



C 2 (rht 
+ W^v) \~2 

v-l 



rht\ 



1 



C 2 fmt\ 



l-v 



l-vT(l-v)\2j 



(3.247) 
(3.248) 



Although the solution has one decaying and one growing mode, it might happen that the 
decaying mode stays larger than the growing mode. To check this we calculate constants 
C\ and C2 by matching the above equations to the values W en d anci Wend at the end of 
inflation (denoted by 'end'). Thus, we find that 



W 



W + HW 



2 / a \ ' 
3w + l \a end 

-2 



i(3to-l) 



H* 



^cnd 



w end + 



w end + 



Wend 

H* 



Wend 

H* 



(3.249) 
(3.250) 



where H* is the inflationary Hubble scale. Plugging these solutions into Eq. (3.243) (and 
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^ = _^L_(*Y7 * ) ^(mtr«l, (3.251) 



using that a 3 ^ 1+w ^ oc t 2 ) we obtain 

2 ( t x2 
Pkin (3w + l) 2 V^*/ \*end 

which implies that the total energy density of the light vector field is 

PW - Pkin = \ ( Wend + Wend-ff*) (—— ) => OC a" 4 . (3.252) 
2 V / V a cnd/ 

Therefore, we see that the energy density of the light vector field scales as that of rela- 
tivistic particles. This is in striking difference to the scalar field case, in which when the 
field is light its density remains constant even after inflation. 

On the other hand, if the vector field is heavy, rht 3> 1, the Bessel functions in 
Eqs. (3.245) and (3.246) are oscillating. Hence, as was discussed in section 3.3.1, the 
heavy vector field oscillates with a frequency much larger than the Hubble parameter 
and with the amplitude decreasing as oc a~ 3 / 2 . In Eq. (3.60) it was shown 

that the energy density of such field decreases as p\y oc a -3 and the average pressure 
is zero, i.e. p± ~ 0. Therefore, on average, the oscillating vector field behaves as 
pressureless isotropic matter and can dominate the Universe without generating excessive 
large scale anisotropy. This is crucial for the vector curvaton mechanism because, to 
produce the curvature perturbation, the field must dominate (or nearly dominate) the 
Universe without inducing excessive anisotropic expansion. 

3.5.5. Curvaton Physics 

In this section we calculate constraints for our vector curvaton model assuming that the 
scaling behavior of f(t) and m{t) ends when inflation is terminated. This implies that 
the scaling is controlled by some degree of freedom which varies during inflation, e.g. the 
inflaton field. 

In the curvaton scenario the total curvature perturbation can be calculated as the sum 
of individual curvature perturbations from the constituent components of the Universe 
multiplied by the appropriate weighting factor. In the current scenario this is written as 
follows 

C = (1 - &w)(~f + &w(w, (3.253) 

where £l\v is defined in Eq. (2.197). As in the scalar curvaton paradigm, the above is to 
be evaluated at the time of decay of the curvaton field. 
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As was discussed in section 3.5.3, if m > H* at the end of inflation, then the vector 
field perturbation spectrum is isotropic and may generate the total curvature perturba- 
tion in the Universe without violating observational bounds on the statistical anisotropy 
of the curvature perturbation. If this is the case, we can assume that £ 7 = 0. On the 
other hand, when rh <C H*, the amplitude of the spectrum of the longitudinal compo- 
nent of the vector field perturbations is substantially larger than the one of the transverse 
perturbations. Hence, the curvature perturbation due to the vector field is excessively 
anisotropic. To avoid conflict with observational bounds, the contribution of the vector 
field to the curvature perturbation has to remain subdominant. Therefore, for this sce- 
nario, we have to consider £ 7 ^ and the curvature perturbation already present in the 
radiation dominated Universe must dominate the one produced by the vector curvaton 
field. 

In Eqs. (3.196) and (3.201) it was shown that the typical value of the field perturbation 
is 5W ~ (3iT*/-M)(iT*/27r). If M <C H* this is because the longitudinal component is 
dominant over the transverse ones (see Eq. (3.195)). If M > H*, then the transverse 
and longitudinal components are oscillating with same amplitudes (see Eq. (3.200)). 

For this reason, at the end of inflation, we can write 

<We„d~^f^^, (3.254) 

where we have taken M = rh and / = 1 at the end of inflation. Wend can be found from 
Eq. (3.235) by using (p w )end ^ W M ~ W end rh (see Eqs. (3.235) and (3.241)). Thus, 



^ end „ V(PwU d (3 255) 

m 

Hence, from Eq. (3.80) we calculate the curvature perturbation of the vector field 

mpi 

where O en( j = (pw / p)end 1S the density parameter of the vector field at the end of inflation, 
p en d is the total energy density dominated by the inflaton field, and we have used the 
Friedman equation: 3m|>i?^ = p cn d- Since the vector field must be subdominant during 
inflation we have C nd ^ 1- 

Eq. (3.256) is valid in both a = — 1±3 cases. The only difference is that, in the / oc a 2 
case, statistically isotropic curvature perturbations cannot be generated. Hence, only 
considerations for statistically anisotropic perturbations in Sec. 3.5.3 are relevant. 
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To calculate the parameter space for this model we note that at the end of inflation the 
inflaton field starts oscillating and w ^ —1. Therefore the Hubble parameter decreases 
as H(t) ~ t . In general, the inflaton potential is approximately quadratic around 
its VEV. Thus, the coherently oscillating inflaton field corresponds to a collection of 
massive particles (inflatons) whose energy density decreases as a" 3 . When the Hubble 
parameter falls bellow the inflaton decay rate V, the inflaton particles decay into much 
lighter relativistic particles reheating the Universe. After reheating, the Universe becomes 
radiation dominated with the energy density scaling as p 7 oc a -4 . 

On the other hand, the evolution of the energy density of the vector field, depends on 
its mass rh. As discussed in Sec. 3.5.4, if rh <C the energy density scales as pw oc a~ 4 
until the vector field becomes heavy and starts oscillating. If rh 3> H*, however, the 
vector field has already started oscillating during inflation and py/ oc a -3 . 

To avoid causing an excessive anisotropic expansion period the vector field must be 
oscillating before it dominates the Universe and decays. This requirement implies that 

r, rh > F w , H dom , (3.257) 

where Ty/ is the decay rate of the vector field and H^om is the value of the Hubble 
parameter when the vector field dominates the Universe if it has not decayed already. 
Working as in Ref. [84], we can estimate -Hdom as 

( - \ 2 /3 (- - ^ -1/6 

#dom ~ ft C ndr 1/2 min j 1; _ J min 1 1; ~ I . (3.258) 

Similarly, if the vector field decays before it dominates, the density parameter just before 
the decay is given by 

r, \ 1/2 ( ^ N 2/3 ( ~ N -1/6 

1 \ ' ml' m 1 ' 



^doc ~ ^end (frrj min 1 1;— | minjl;— \ . (3.259) 

where Od ec = (^lw)dec- Combining the last two equations and using Eq. (3.256) we can 
express the inflationary Hubble scale as 

— ^dcc^mmjl;-j mm|l;-| ^ j . (3.260) 

The bound on the inflationary scale can be obtained by considering that the decay rate 
of the vector field is Tw ~ h 2 rh, where h is the coupling to the decay products. Then 
we can write max{T\Y] H dom } > h 2 rh. Furthermore, we must consider the possibility 
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of thermal evaporation of the vector field condensate during the radiation dominated 
phase. If this were to occur, all the memory of the superhorizon perturbation spectrum 
would be erased. The bound on h, such that the condensate does not evaporate before 
its decay, is given in Eq. (3.89). 

From Eq. (3.260) one can see that the parameter space is maximized if the Universe 
undergoes prompt reheating after inflation, i.e. if T — >■ H*. To find the parameter space 
we investigate two separate cases: when rh S> H* and when rh <C H*. 

3.5.5.1. The Statistically Isotropic Perturbation 

The statistically isotropic perturbation can be realized only in the case when a = —4. 
As mentioned before, if the mass of the vector field at the end of inflation is larger than 
the Hubble parameter, rh > .H*, then the field has started oscillating already during 
inflation. In this case amplitudes of the longitudinal and transverse perturbations are 
equal and therefore the curvature perturbation induced by the vector field is statistically 
isotropic. We can assume, in this case, that the vector field alone is responsible for 
the total curvature perturbation in the Universe without the need to invoke additional 
perturbations from other fields. Thus, we can set £ 7 = in Eq. (3.253) and write 



Using this and the lower bound on h we find from Eq. (3.89) the lower bound for the 
inflationary Hubble parameter 



where we have taken into account that the parameter space is maximised when the 
Universe undergoes prompt reheating, i.e. T — > H*. From this expression it is clear 
that the lowest bound is attained when the vector field dominates the Universe before 
its decay, Odec —> U an d when the oscillations of the vector field commence at the very 
end of inflation, i.e. rh — > H*. With these values we find the bounds 



where V* denotes the inflationary energy scale and we used that (~5x 10~ 5 from the 
observations of the Cosmic Background Explorer. 

In view of the above, we can obtain a lower bound for the decay rate of the vector 



C ~ ^dccCw ■ 



(3.261) 




(3.262) 



> 10 9 GeV & V} /A > 10 14 GeV 



(3.263) 
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field. Indeed, using Eqs. (3.89) and (3.263) we find 

,3 H 3 



m 



rv> ^r>^>io- 9 GeV. 



(3.264) 



From the above we find that the temperature of the Universe after the decay of the 
vector field is Td ec ~ V m Pl^W ^ 10 4 GeV, which is comfortably higher than the temper- 
ature at BBN Tbbn ~ 1 MeV (i.e. the decay occurs much earlier than BBN), and also 
higher than the electroweak phase transition, i.e. the decay precedes possible electroweak 
baryogenesis processes. 

Since m > H*, Eq. (3.263) corresponds to a lower bound on m. An upper bound on 
m can be obtained as follows. Because, m > > T, Eq. (3.259) becomes 



n 



dec ~ ^end 



W 



(3.265) 



From Eq. (3.89) we have T\y > m 3 /m^. Combining this with the above we obtain 



3 < ( ^cnd\ p 
~ V^dcc/ 



PI- 



(3.266) 



Now, when a = —4 we have M oc a 3 during inflation. Since the end of scaling occurs 
when inflation is terminated, for a < a en d we can write 



m 



( 



«end \ ' 



V a 



M ~ e' 



3N os 



(3.267) 



where we considered that the field begins oscillating when M ~ and N osc is the 
number of remaining e-folds of inflation when the oscillations begin. Inserting the above 
into Eq. (3.266) we find 



N < /v max 

iv OSC rO iv osc 



9 



( ^dcc ) V-H* 



J3* 



2, 
< - In 



mpi 



9 V ^I^qqH^ 



(3.268) 



where in the last inequality we used that ^ C nd < 1 an d T < H*. Now, considering that 
m>H*, Eq. (3.262) gives 

( ) , 11 

> C 2 . (3.269) 



mpi 



Hence, combining Eqs. (3.268) and (3.269) we obtain 



N. 



< __ w = 4 4 

osc - 9 



(3.270) 
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Thus, in view of Eq. (3.267), we obtain the bound rh < e 3 ^"ff», which results in the 
following parameter space for rh: 

1 < rh/H* < 10 6 , (3.271) 

where we used Eq. (3.270). The above range is reduced if the decay of the curvaton occurs 
more efficiently than through gravitational couplings, i.e. if h > m/mp\. Nevertheless, 
we see that the parameter space in which the vector field undergoes isotropic particle 
production and can alone account for the curvature perturbation, is not small but may 
well be exponentially large. Indeed, repeating the above calculation with Tw ~ Th (i.e. 
h ~ 1) it is easy to find that 




(3.272) 



Hence, using that fi e nd < 1 an d T < H* we obtain 

^c X = -^lnO dcc <3.1 l<rh/H*<W\ (3.273) 

where we used fidec > 10~ 2 . This is because, in the case considered, /nl is given by 
Eq. (3.211), so a smaller fidcc would violate the current observational bounds on the non- 
Gaussianity in the CMB temperature perturbations (see the discussion in section 2.1.2.2). 

Still, it seems that, to obtain an exponentially large parameter space for m, we need 
pw not to be too much smaller that during inflation and also inflationary reheating 
to be efficient. In the case of gravitational decay (Ty/ ~ rh^/rnp^) Eq. (3.268) has a weak 
dependence on both fi e nd and T: m oc (^cnd^) 1 ^ 3 ' wn i cn means that the allowed range of 
values for rh remains large even when fiend an d T are substantially reduced. This is not 
necessarily so when T\y ~ h 2 rh, with h 3> rh/mp\. Indeed, in this case it can be easily 
shown that rh oc h~ 2 Q 2 nd T. Therefore, if T is very small it may eliminate the available 
range for rh. Fortunately, the decay coupling h can counteract this effect without being 
too small. 

3.5.5.2. Statistically Anisotropic Perturbations 

If the vector field is not responsible for the total curvature perturbation in the Universe, 
the parameter space is more relaxed. In this case, the vector field may start oscillating 
after inflation and hence its mass is m C H*. However, this means that the curvature 
perturbation due to the vector field is strongly statistically anisotropic. For this reason 
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we can no longer set £ 7 to zero in Eq. (3.253) because the curvature perturbation present 
in the radiation dominated Universe must be dominant. In other words, the parameter 
£ defined in Eq. (3.31) needs to be very small, 1. 

In this case the total curvature perturbation is given in Eq. (3.86). Inserting this into 
Eq. (3.260) and considering again that the lowest decay rate of the vector field is through 
the gravitational decay, max{T\y; ifdom} > m 3 /mpj we find 

H* (ge\ m ( m A 5/8 ( r V 3/8 . f ml 1/8 



mpi 



> (^X ( A) ( JLV rnin(l--V (3 274) 

V^dcc/ \mpi) \mpij \ 'rj 



The above suggests that the lower bound on is minimised for prompt reheating with 
r — > H*. Also, from observations we know that the statistically anisotropic contribution 
to the curvature perturbation must be subdominant. Thus, the vector field should not 
dominate the Universe before its decay. Hence, using T — > H* and < 1 we obtain 



— >y/9C\—- (3-275) 
mpi y mpi 

From this expression it is clear that the parameter space for H* is maximised for the 
lowest mass value. The minimum mass of the vector field can be estimated from the 
requirement that the field decays before BBN. Because the lowest decay rate is the 
gravitational decay, this condition reads m 3 /m P1 > Tg BN /mpi, with Tbbn ~ lMeV, 
which corresponds to rh > 10 4 GeV. Using this, we find that the parameter space for 
the vector curvaton model with the statistically anisotropic curvature perturbations is 

> g 1 ' 2 10 7 GeV ^ U* 1/4 > g 1 ^ 10 13 GeV , (3.276) 

i.e. it is somewhat relaxed compared to the statistically isotropic case (c.f. Eq. (3.263)) 
depending on the magnitude of the statistical anisotropy in the spectrum, for which 
g < 0.3 (see the discussion above Eq. (2.19)). This result is valid for both a = — 1 ± 3 
cases. From the above it is evident that there is ample parameter space for the mass of 
the vector field 

lOTeV <m<Z.H*. (3.277) 
3.5.6. Summary for the Massive fF 2 Model 

In section 3.5 we studied a particularly promising vector curvaton model consisting of a 
massive Abelian vector field, with a Maxwell type kinetic term and with varying kinetic 
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function / and mass m during inflation. The model is rather generic, it does not suffer 
from instabilities such as ghosts and may be realized in the context of theories beyond 
the standard model such as supergravity and superstrings (see two tentative examples in 
Ref. [107]). 

We have parametrised the time dependence of the kinetic function as / cx a Q , where 
a = a(t) is the scale factor. Our model offers two distinct possibilities. If rh < H* 
(possible for a = — 1 ± 3) the vector field can only produce a subdominant contribution 
to the curvature perturbation £, but it can be the source of statistical anisotropy in 
the spectrum and bispectrum. In fact, non-Gaussianity in this case is predominantly 
anisotropic, which means that, if a non-zero /nl is observed without angular modulation, 
then our model is falsified in the rh < case. The second possibility (possible for 
a = —4 only) corresponds to rh > H*. In this case the vector field can alone generate 
the curvature perturbation £ without any contribution from other sources such as scalar 
fields. If ra > H*, particle production is isotropic and the model does not generate any 
statistical anisotropy. The vector field begins oscillating a few e-folds before the end of 
inflation but its density remains constant until inflation ends. The parameter space for 
this case can be exponentially large, i.e. 1 <C rh/H* < 10 6 . Significant non-Gaussianity 
can be generated, provided the vector field decays before it dominates the Universe, in 
which case /nl is found to be identical to the scalar curvaton scenario. In other words, if 
rh 3> , our vector curvaton can reproduce the results of the scalar curvaton paradigm. 
Finally, if rh ~ if* the vector field can alone generate the curvature perturbation £ but it 
can also generate statistical anisotropy in the spectrum and bispectrum. In this case, the 
anisotropy in /nl is subdominant and equal to the statistical anisotropy in the spectrum, 
which is a characteristic signature of this possibility. However, the allowed range for rh 
values in this case is very narrow, as shown in Eq. (3.214), requiring accurate tuning 
of the initial conditions. We have also found that inflation has to occur at energies of 
V^ A > 10 14 GeV in the (almost) isotropic and V*^ > g 1//4 10 13 GeV in the anisotropic 
case. 

3.6. The End-of-lnflation Scenario 

3.6.1. Vector Field Perturbations and £ 

In this section we consider another model in which a vector field influences the generation 
of the curvature perturbation. The model avoids excessive large scale anisotropy in 
the Universe by a different mechanism than the vector curvaton scenario described in 
section 3.3. The idea is based on Ref. [58] which was summarized in section 2.4.2, where 
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it was shown that in hybrid inflation models the generation of the curvature perturbation 
can be realized due to the inhomogeneous end of inflation. Yokoyama and Soda [89] used 
this idea to generate the anisotropic contribution to the total curvature perturbation. 
In their model the anisotropy is generated at the end of inflation due to the vector field 
coupling with the waterfall field. In other words the scalar field a of section 2.4.2 is 
changed by the vector field A„. In this section we calculate the non-Gaussianity of the 
model in Ref. [89] using the formalism developed in section 3.2. 

This scenario uses the conformal invariance breaking of the U (1) vector field through 
the non-canonical kinetic function of the form / (t) F^F^: 

S = J ^pD g (~f (t) F^ d 4 x, (3.278) 

where F^ u = d^A u - d v A. 

This action is only written for the conformal invariance breaking term, the dots rep- 
resent other terms which give inflation with practically constant H, and generate / (i) 
without having any other effect on the evolution of the gauge field during inflation. For 
the vector field to be gauge invariant any scalar field coupled to must have zero expec- 
tation value (no spontaneous symmetry breaking) with negligible quantum fluctuation 
around that value. 

This form of the conformal invariance breaking was considered in many papers. Start- 
ing from Ref. [108] such action was often considered for the generation of the pri- 
mordial magnetic fields (see for example Refs. [109, 110, 111, 112]) and recently in 
Refs. [84, 107, 113] it was considered for the generation of £ by the vector field (see 
section 3.5). In these papers it was discovered that a scale invariant perturbation spec- 
trum of the physical, canonically normalized vector field is obtained if / oc a 2 (we 
showed in section 3.5 that this is the case for / oc a~ 4 as well): 

V + ={^j. (3.279) 

In the end-of-inflation scenario of Soda and Yokoyama there are two components of the 
curvature perturbation: one generated during inflation and an anisotropic one, generated 
by a vector field at the end of inflation: 

C = Cinf + Cend- (3.280) 

The first component is due to the perturbation of the light scalar field, while the second 
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one is due to the perturbation of the vector field with the kinetic term in Eq. (3.278). 
Without parity violating terms the power spectra for left handed and right handed polar- 
izations are equal, while the longitudinal polarization is absent for a massless vector field. 
In this situation we find that parameters p (k) and q (k) defined in Eq. (3.19) become 

p=-l and q = 0. (3.281) 

Cinf in this scenario is the statistically isotropic contribution to the total curvature 
perturbation. In the slow roll inflation the spectrum of the scalar field perturbation is 
= {H/2-k) , so that we get 

V+ = V4,, (3.282) 
And the total isotropic part of the curvature perturbation from Eq. (3.30) becomes 

Vf° = V+Nl (1 + , (3.283) 

with £ given by 

{ ^) 2 . (3.284) 

Using the expression for the anisotropy parameter g in Eq. (3.32) we find that in this 
scenario 

g = - T ^. (3.285) 

Taking into account Eq. (3.281), the vector Mi (k) defined in Eq. (3.35) reduces to the 
simple form 

M (k) = NwVtp N w - k (n w ■ k) . (3.286) 
3.6.2. Hybrid Inflation Model 

To calculate /nl we consider a specific example of the hybrid inflation with the potential 
V (4>, X, A^) = V + l -ml4> 2 - l -m 2 xX 2 + \^ + + \\aX 2 A»A^ (3.287) 

which contributes to terms in Eq. (3.278) denoted by dots. Here 4> is the inflaton and \ 
is the waterfall field (compare this with the scalar field case in Eq. (2.178)). The effective 
mass of the waterfall field for this potential is 



^eff = -ml + A^> 2 - —WW, (3.288) 

J 
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where we chose the Coulomb gauge with W% = and d{W l = was chosen. Inflation ends 
when the inflaton reaches a critical value (fi c where the effective mass of the waterfall field 
becomes tachyonic. But one can see from Eq. (3.288) that the critical value is a function 
of the vector field cp c = <j) c (W). With this in mind Nfc and TV™- can be readily calculated: 



dN 



w 



dW,. 



fX 4 



(3.289) 



and 



dN d 2 (t) c 
d6 r dWidW, 



+ 



d 2 N 



N 2 



dWi dWj <t> c N c 



C 5j, 



(3.290) 



where we have defined 



dN 



and C 



f K 



X A W' 



(3.291) 



where W = |Wj| and / are evaluated at the end of inflation and we used the fact that 
Ncc/N 2 ~ N^fi/N^ ~ 0(e) under the slow roll approximation, where e is the slow 
roll parameter defined (Vf/Vy (see Eq. (1.55)), with the prime denoting 

derivatives with respect to the inflaton. As mentioned earlier the total of perturbations 
consists of two components: perturbations of the scalar and vector fields. This gives the 
following bispectrum in the equilateral configuration 



B^Ckx.k^ka 



B 



cquil (k 1 ,k 2 ,k 3 )+^ uil (k 1 ,k 2 ,k 3 



WlN 2 N AA + 



Mi (ki) N&Mj (k 2 ) + c.p. 



VlNi 



e 

N c <f> c 



(C 2 



C 2 - 1 )w 2 



16 



(3.292) 



The mixed term B^y^ 1 is absent from Eq.(3.292) because in this model N^ w = 0. By 
using the expression for the isotropic power spectrum in Eq. (3.283) and the bispectrum 
in Eq. (3.292) from the definition of f^ 11 in Eq. (3.38) we obtain 



6 /.equil 

5 /nl 



(c 2 - 1) 



-c 2 -i)w 2 - —wi 

8 j x 16 x 



(3.293) 



where the slow parameter rj is equal to rj = mpjV^/y = m^mpj/Vo and mp\N c = l/^/2e c , 
with e c being the e parameter evaluated at the end of inflation. Similarly, for the squeezed 
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configuration we find 




W 2 (C 2 -l)-(C 2 -l)^i-^( S in^)M 



(3.294) 



In this equation <p is the angle between the vectors ki and Wj_ (see Figure A.l). 

We find that /ni!" 1 ana - /nl^ are f unc ti° ns of W±, i.e. they are anisotropic and corre- 
lated with the statistical anisotropy. Also the level of non-Gaussianity is proportional to 
the anisotropy parameter squared, /nl oc g 2 , as in the vector curvaton model. However, 
the angular modulation of /nl in this scenario is different from the curvaton scenario. 
From Eqs. (3.293) and (3.294) we see the additional modulation term proportional to 
Wj_. This term is absent in the vector curvaton scenario. 

As was mentioned earlier, in this model the vector field during inflation is massless 
and, therefore, gauge invariant. The homogeneous value of such vector field can be set 
to zero by an appropriate gauge choice. However, as seen from Eqs. (3.291) and (3.293), 
(3.294) calculated predictions do depend on the homogeneous value of the vector field 
W. Therefore, for this model as it stands, the interpretation of the results are not 
clear. Although terms proportional to Wj_ in /nl expressions do not depend on C and 
consequently on the gauge choice. 

3.6.3. Summary of the End-of-lnflation Scenario 

In section 3.6 we have considered a model proposed in Ref. [89]. In this model the 
energy density of the vector field is subdominant throughout the history of the Universe. 
However, it influences the generation of £ by modulating the end of inflation through 
the coupling to the waterfall field. The conformal invariance of the massless vector field 
is broken by the time dependent kinetic function as in section 3.5. We consider a scale 
invariant perturbation spectrum of the vector field with the kinetic function scaling as 
/ oc a 2 . In this model the vector field is gauge invariant, therefore its particle production 
is anisotropic and the curvature perturbation generated due to this field is statistically 
anisotropic. 

We have calculated the non-linearity parameter /nl in this model and found that it is 
correlated with anisotropy in the power spectrum, as in the vector curvaton scenario. In 
this model too /nl has an angular modulation with the amplitude of the same order as 
the isotropic part. In addition it has the modulation term, proportional to Wf, which is 
absent in the vector curvaton model. 
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The successes of the standard Hot Big Bang theory in explaining the structure and 
evolution of the Universe since the very first second until today is very impressive. The 
predictions for abundances of the light elements are in a very good agreement with 
observations. The origin and the process of formation of galaxies and galaxy clusters 
are now well understood. However, to reproduce the observable Universe, the initial 
conditions of the HBB model must be finely tuned. The spatial curvature of the Universe 
must have been incredibly close to zero near its birth for the Universe to have time to 
evolve to the present state, and it must have started being extraordinary smooth even in 
regions which were never in causal contact with each other. In addition, in the framework 
of the standard HBB model, there are no mechanisms to explain the origin of tiny 
primordial density perturbations which are almost Gaussian, adiabatic and correlated on 
superhorizon scales. Such perturbations are observed as temperature fluctuations in the 
CMB sky and they seed the growth of large scale structure. 

The fine tuning problems may be substantially alleviated by postulating a period of 
accelerated expansion at the earliest stages of the evolution of the Universe. This period 
is called inflation. In addition to solving the flatness and horizon problems, the greatest 
achievement of the inflationary paradigm is the explanation of the origin of the primordial 
density perturbation which has the properties observed in the CMB sky. According to 
this paradigm, the primordial density perturbation originated as quantum fluctuations 
during the inflationary period. In Chapter 2 we have demonstrated how the application 
of quantum field theory on a curved space-time background may lead to the amplification 
of quantum fluctuations and their conversion into the classical field perturbation. This 
perturbation, consequently, causes the perturbation in the curvature of space-time. Much 
later, after inflation, when the wavelengths of the perturbation become smaller than the 
horizon size, it seeds the formation of structure in the Universe due to the process of 
gravitational instability. 

To describe the formation and evolution of the cosmological perturbation we have 
used a very important quantity: the curvature perturbation £. This quantity is constant 
throughout the history of the Universe, except during those periods when the total pres- 



149 



4. Summary and Conclusions 



sure of the Universe is not a unique function of the energy density. In other words, when 
pressure is not adiabatic. To show how the classical field perturbation, originating from 
quantum fluctuations, is related to the curvature perturbation £, we used the separate 
universes approach. In this approach the evolution of the Universe on superhorizon scales 
at each space point is treated as that of the separate, unperturbed Universe with the 
locally defined expansion rate. The latter is determined by the average energy density on 
the flat hypersurface at that point, where the averaging is performed on a superhorizon 
scale of interest. 

The statistical properties of £ provide one of the main tools in cosmology for observa- 
tional tests of models of the very early Universe. We have shown how these properties 
may be calculated using the 5N formalism. It was applied to calculate the power spec- 
trum and the bispectrum at tree level for three models, namely: the single field inflation, 
the end-of-inflation and the curvaton scenarios. In the treatment of these three models, 
we have assumed that £ is generated solely by quantum fluctuations of scalar fields. In 
Chapter 3 we showed that quantum fluctuations of vector fields may contribute or even 
generate the total curvature perturbation in the Universe as well. 

However, a massless, canonically normalized U (1) vector field cannot produce £ be- 
cause, being conformally invariant, its quantum fluctuations are not amplified during 
inflation. And even if they were amplified, such a field cannot dominate the Universe 
without producing excessive large scale anisotropy, i.e. excessive anisotropic expansion of 
the Universe, although in most scenarios the vector field must dominate or nearly dom- 
inate the Universe to generate £. In section 3.1 we discuss possibilities of breaking the 
conformal invariance of vector fields and introduce four mechanisms for the generation 
of (" by vector fields without producing an excessive large scale anisotropy. 

In section 3.2 we have extended the 5N formalism to include perturbations of vector 
fields. In contrast to the scalar field, which has one degree of freedom (DoF), the massive 
vector field has three DoF. Therefore, in a theory with a massive vector field we must 
consider quantum fluctuations for all three of them. To calculate the evolution of each 
DoF they were decomposed into the longitudinal and two circular polarization vectors. 
This choice is advantageous because each polarization vector transforms differently under 
the Lorentz group. Therefore, we can be sure that they do not mix in the course of 
evolution. We found that in general the amplification of quantum fluctuations is not the 
same for all three DoF. In other words, the particle production of a vector field is not in 
general isotropic. This results in different values of n-point correlation functions for each 
polarization. 

To quantify the anisotropy in the particle production we introduced two parameters 
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p(k) and q(k) in Eq. (3.19), where k is the wavevector. The q parameter quantifies the 
difference in the power spectra of two transverse polarization modes. It is non-zero only 
in parity violating theories. The parameter p quantifies the difference in the longitudinal 
power spectrum and the average of the transverse ones. If both parameters are equal to 
zero, the particle production of the vector field is isotropic. However, if any of these are 
non-zero, the particle production is anisotropic. The values of p and q parameters are 
determined by the mechanism which brakes the conformal invariance. 

If the vector field with anisotropic particle production generates or affects the curvature 
perturbation, the latter is statistically anisotropic, i.e. statistical properties of £ are not 
invariant under rotations. The power spectrum of such perturbation will have an angular 
modulation. To the lowest order we can express it as [36] 



where 7^° is the isotropic part of the spectrum, n is the unit vector along the preferred 
direction and g parametrizes the amount of modulation. In the vector field models n is 
in the direction of the homogeneous vector field. V 1 ^ in these models may be solely due 
to the vector field, if g satisfies the observational bounds, or it may be dominated by 
some other, statistically isotropic source of Q. The present observational bound on the 
anisotropy in the spectrum of £ is g < 0.3 (see the discussion above Eq. (2.19)). The value 
of g is determined by the mechanism which generates the curvature perturbation and by 
the value of p. In this thesis we consider two such mechanisms: the vector curvaton and 
the end-of-inflation scenarios. 

The vector curvaton scenario, first proposed in Ref. [83], uses the fact that a heavy 
vector field oscillates rapidly with the frequency much larger than the Hubble parameter. 
The time averaged pressure of such field is zero and the energy density decreases with 
the scale factor as o -3 . Thus, the heavy vector field acts as pressureless, isotropic matter 
and can dominate the Universe without producing excessive large scale anisotropy. In 
accord with the curvaton scenario, the vector field dominates (or nearly dominates) the 
Universe after reheating, when the latter is radiation dominated. The curvaton imprints 
its perturbation spectrum and decays before the BBN. The perturbation spectrum of the 
vector curvaton field is acquired during inflation, when the field is light and its energy 
density is negligible compared to the inflaton one. During this period the values of 
parameters p and q are determined, depending on the mechanism of conformal invariance 
breaking. 

In section 3.3.2 the general predictions for the non-linearity parameter /nl are derived 




(4.1) 
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in the vector curvaton scenario with p / and q ^ 0. First, we find that /nl has an 
angular modulation, similarly to the power spectrum. The amplitude of the modulation 
is parametrized by Q given by 

hh = Ml{l + QWl), (4.2) 

where W± is the projection of the unit vector of the preferred direction onto the plane 
of vectors ki, \s.2 and k3 which were used to calculate the bispectrum. The preferred 
direction is determined by the direction of the homogeneous vector field. Therefore, 
we find that both, the power spectrum and /nl have the same direction of angular 
modulation. Another important prediction of the vector curvaton scenario is that the 
magnitude of /nl is correlated with the anisotropy in the power spectrum, i.e. oc g 2 . 

We calculated /nl in the squeezed and equilateral configurations (Eqs. (3.71), (3.72) 
and (3.75), (3.76)) and found that only the equilateral configuration is sensitive to the 
parity violating terms in the Lagrangian of the theory. If the theory is parity conserving, 
isotropic parts of /nl are equal in both configurations. Therefore, the detection of 
different values of fi|? in the squeezed and equilateral configurations would indicate 
parity violation. However, the amplitude of the angular modulation Q is not equal 
in the squeezed and equilateral configurations for both parity violating and conserving 
theories. In addition, the anisotropic part of /nl dominates over the isotropic part 
if p > 1. Although presently there are no observational constraints on the values of 
Q in the squeezed and equilateral configurations, the detection of them would allow a 
unique determination of p and q, and therefore, would constraint very tightly the possible 
conformal invariance breaking mechanisms for the vector field during inflation. 

If, on the other hand, the particle production is isotropic, i.e. p = and q = 0, the 
predictions of the vector curvaton scenario do not differ from the standard scalar curvaton 
case. However, this offers a possibility to generate the total curvature perturbation in 
the Universe solely by the vector field, without directly invoking scalar fields at all. But 
even if the particle production is anisotropic with \p\ < 0.3 and any value of q, the 
vector field can still generate the total £ with the amount of statistical anisotropy satisfy 
observational bounds. 

To find the values of p and q, we consider two mechanisms of breaking the conformal 
invariance. In the first one a massive Abelian vector field is non-minimally coupled to 
gravity through the Ricci scalar, see Eq. (3.90). We calculate the perturbation power 
spectra for all three polarizations and find that the scale invariance is achieved if the 
non-minimal coupling constant is equal to 1/6 and the bare mass of the vector field 
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is much smaller than the Hubble parameter. Because the given Lagrangian is parity 
conserving, the parity violation parameter is q = 0. The other anisotropy parameter is 
p = 1 in this model. As was discussed after Eq. (3.32), because p > 0.3 such a vector 
field cannot produce the total curvature perturbation in the Universe, without violating 
observational bounds on statistical anisotropy. Therefore, the dominant contribution to 
£ must come from some other, statistically isotropic source. In the context of the vector 
curvaton scenario this means that the vector field must decay before dominating, while 
the dominant contribution to £ must be present in the radiation dominated background 
before the curvaton decay. 

We find that, in non-minimally coupled vector curvaton model, the isotropic part of 
/nl is equal to 2g 2 where Q\y < 1 is the density parameter of the vector field just 
before its decay. The amplitudes of anisotropic parts are 1 and 9/8 times the isotropic 
part in the squeezed and equilateral configurations respectively. After taking into account 
all cosmologically relevant bounds we find that the parameter space for this model is 

H > g 10 10 GeV ^ V 1/4 > g 1/2 10 14 GeV, (4.3) 

where H is the inflationary Hubble parameter and V 1 ^ is the energy scale of the inflation. 
From this result it is clear that the parameter space is large enough for a successful 
realization of this scenario in particle physics models. 

Another model considered in this thesis is of the vector curvaton with time dependent 
kinetic function and mass, see Eq. (3.169). As in the previous model we calculate the 
superhorizon perturbation spectra for all three polarizations and find that they are scale 
invariant if the mass varies with the scale factor as m oc a and is smaller than the Hubble 
parameter when cosmological scales exit the horizon, while the kinetic function scales as 
/ oc a~ 1±3 . 

We assume that degrees of freedom, which modulate the time dependence of the kinetic 
function and mass, are stabilized at the end of inflation. Therefore, the vector field mass 
becomes constant at that moment, i.e. m = constant = rh. Since any constant value in 
front of the kinetic function may be absorbed into the definition of the vector field, we 
may set / = 1 at the end of inflation, and the field becomes canonically normalized. As 
we saw, the scale invariant perturbation spectra may be achieved if the kinetic function 
/ is increasing as well as decreasing. If the vector field is a gauge field, then / is the 
gauge kinetic coupling. In this case it is inversely proportional to the gauge coupling e 
as / oc e~ 2 . Therefore, an increasing / oc a 2 (small during inflation) would correspond 
to the strongly coupled regime, while / oc a" 4 would correspond to the weak coupling. 
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Therefore, only the second case may be realized in the particle physics models. 

First we calculate the anisotropy in the particle production for / oc a~ 4 and find that it 
depends on the mass of the vector field at the end of inflation m. Since the Lagrangian of 
this model has no parity violating terms, q = 0. But the value of p depends on in. If the 
vector field is light at the end of inflation p ^ and is given by p = (3H/m) 2 — 1. If, on 
the other hand, the vector field is heavy, p = 0. Therefore, the light vector field generates 
the statistically anisotropic curvature perturbation, while the heavy field generates the 
statistically isotropic one. 

The isotropic part of /nl in the p » 1 case is equal to /^£ = (2g 2 /Qw) ■ (3H/m) A . The 
amplitude of the /nl angular modulation in this regime is {3H/mf and | {3H/mf times 
larger than the isotropic part in the squeezed and equilateral configurations respectively. 

If the vector field is heavy at the end of inflation, p = and the generated curvature 
perturbation is statistically isotropic. Such a vector field may generate the total curvature 
perturbation in the Universe without the need of scalar field contribution. In this regime 
the standard curvaton scenario predictions for the non-Gaussianity are valid, i.e. if the 
curvaton decays before domination /nl ~ 3/ (2Qw). In the opposite case, when it decays 
being dominant, the generated £ is Gaussian. 

For this model, when the vector curvaton is light at the end of inflation, the allowed 
range of inflationary Hubble parameter and energy scale is 

H > g 1 ' 2 10 7 GeV V 1/4 > g 1/4 10 13 GeV, (4.4) 

while the allowed region for the vector field mass at the end of inflation is 

10TeV<m<#. (4.5) 

For the heavy field, and consequently statistically isotropic curvature perturbation, the 
analogous bounds are 

H > 10 9 GeV ^ V 1 / 4 > 10 14 GeV (4.6) 

and 

1 < m/H < 10 6 , (4.7) 

were we considered that the vector field produces the total curvature perturbation. Al- 
though for the statistically isotropic case the parameter space is somewhat reduced, in 
both cases it is large enough for a successful implementation in realistic particle physics 
models. 
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So far we have discussed only the case / oc a~ 4 . In the case of increasing kinetic 
function with / oc a 2 , the same results apply, but the vector field has to be light and can 
produce only statistically anisotropic £. 

In the final section 3.6 of this thesis we calculate the non-Gaussianity in the end-of- 
inflation scenario introduced in Ref. [89]. The conformal invariance of the vector field in 
this model is broken by the time varying kinetic function, similarly as in the vector curva- 
ton case discussed above. However, in this model the vector field is massless; therefore, it 
has only two degrees of freedom and particle production is necessarily anisotropic. In the 
end-of-inflation scenario the vector field is always subdominant. However, it influences 
the generation of £ through a coupling to the waterfall field of hybrid inflation. In this 
way the end of inflation is spatially modulated by the vector field, i.e. the hypersurface 
of the synchronous end of inflation does not coincide with the uniform density hypersur- 
face. We calculated the non-Gaussianity for this model and found that the amplitude 
of angular modulation of /nl is larger than the value of the isotropic part, as in the 
curvaton scenario. However, in contrast to the curvaton scenario, in this model /nl has 
an additional modulation, proportional to Wj_, where the latter is the projection of the 
preferred direction onto the plane of three k vectors, used to calculate the bispectrum. 

In summary, we have shown that a vector field can influence or even generate the total 
curvature perturbation in the Universe. If the particle production of the vector field 
is isotropic, the generated curvature perturbation by such field is statistically isotropic. 
Then the vector field may generate the total £ in the Universe without the direct involve- 
ment of scalar fields. In this case observational predictions for the curvature perturbation 
are the same as for models with scalar fields. If, on the other hand, the particle pro- 
duction of the vector field is anisotropic, the generated contribution to £ by such field 
is statistically anisotropic. In this case observational signatures, very distinct from the 
scalar field case, will be present: anisotropic power spectrum and /nl, where the mag- 
nitude and the preferred direction of the latter is correlated with the anisotropy in the 
spectrum. 

Until recently CMB analyses were performed assuming statistical isotropy of the curva- 
ture perturbation. Our results suggest a new observable: statistical anisotropy. There- 
fore, it is desirable to reanalyze CMB maps without imposing rotational invariance a 
priory. Although current measurements might not be sensitive enough to constraint the 
statistical anisotropy in the primordial curvature perturbation (see Refs. [38, 43]), with 
an advent of the Planck data the situation will improve considerably. For example, ac- 
cording to Ref. [114] the lowest detectable value of g from WMAP data is \g\ > 0.1. With 
an expected performance of the Planck satellite this bound will be reduced to \g\ > 0.02. 
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Planck measurements will be much more sensitive to non-Gaussianity as well. Current 
WMAP bound is |/nl| ^ 100, in case of no detection with Planck data it will be reduced 
to |/nl| ^ 5, very close to the cosmic variance limit [42]. Presently there are no obser- 
vational constraints on the angular modulation of /nl- With such increase in sensitivity 
of measurements in a very near future one expects that anisotropy in the spectrum and 
bispectrum will be discovered or constrained very tightly. In case of the discovery, with 
the magnitude and anisotropy of /nl proposed above, it will be a smoking gun for a 
vector field contribution to the primordial curvature perturbation. 

Another important advancement in this direction will be the confirmation or falsifica- 
tion of the presence of the "Axis of Evil", which suggests that low multipoles of the CMB 
are aligned along one direction [74] . Presently the statistical significance of the discovery 
of the "Axis of Evil" is still debatable. However, its confirmation will have profound 
implications: this would prove the existence of the preferred direction in the Universe. 
Such direction cannot be accounted for by scalar fields, but for vector fields it is natural. 

It is necessary for vector field models to be confronted with observations, in addition, 
the treatment presented in this thesis should be extended in several directions. We 
have investigated only two mechanisms of breaking the conformal invariance of massless 
Abelian vector fields: non-minimal coupling to gravity and the time varying kinetic 
function. In the literature on primordial magnetic fields there are many more mechanisms 
proposed to brake this invariance. It would be desirable to explore which of them may 
give scale invariant perturbation spectra from quantum fluctuations of vector fields. 

Even more so, one would also like to understand the generation of perturbations from 
vacuum fluctuations in the anisotropically inflating Universe. In the scalar field dom- 
inated Universe it is natural to assume isotropic expansion, provided inflation lasted 
long enough before cosmological scales exit the horizon, so that according to the no-hair 
theorem, initial anisotropy was inflated away. In the presence of light vector fields, the 
backreaction on the expansion of the Universe might not be negligible, generating the 
large scale anisotropy. We neglected such backreaction in vector curvaton models because 
the vector field energy density is negligible during inflation. However, the anisotropic ex- 
pansion can be easily accommodated within these models or it might be obligatory in 
others. 

We considered two scenarios for the generation of the curvature perturbation: vector 
curvaton and end-of-inflation. However, the developed formalism may be easily extended 
to include other scenarios that have already been explored for the contribution of the 
scalar field perturbation. 

Three toy models were presented in this thesis to generate the curvature perturbation 
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by vector fields. Ultimately any model of the early Universe must be firmly rooted in 
realistic particle physics theories. In the context of inflationary model building 'particle 
physics theories' can mean two things. It may be the realization of the inflationary 
expansion of the Universe in some fundamental theory, currently the best developed 
of which is string theory. For the cosmological aspects of string theory one may see 
Refs. [115, 116] and references therein. Or particle theory may mean an effective field 
theory, which accurately describes the Nature at the energies when cosmological scales 
exit the horizon. In this direction an extensive effort exists in explaining the inflationary 
epoch in the context of supersymmetry and supergravity (for reviews see Refs. [117, 118]). 

From this point of view a particularly attractive is the vector curvaton model with time 
varying kinetic function / (t) and mass m (t), presented in section 3.5. These functions, 
/ and m, cannot have an explicit time dependence but must be modulated by some 
dynamical DoF during inflation. It might be an inflaton itself, or some other field. A 
vector field with such Lagrangian is very natural in string theory, where parameters such 
as masses and kinetic functions are modulated by scalar fields called moduli. Moduli are 
not fundamental scalar fields, they parametrize the size and shape of the manifold on 
which extra dimensions are compactified. But from our four dimensional perspective they 
act as scalar fields. In Ref. [107] it was shown that the modulus field with an exponential 
potential (which is reasonable for a modulus field) can play a role of a single DoF driving 
inflation as well as modulating the time dependence of the vector field kinetic function 
and mass. 

From the effective field theory point of view, the time varying kinetic function and 
mass is very general in supergravity theories. In this case / is the gauge kinetic function, 
which is a holomorphic function of the scalar fields of the theory. In supergravity the 
potential of these scalar fields receive a correction from the Kahler potential such that 
their mass become m^, ~ H 2 [119, 120, 121]. In the inflationary model building this 
is known as the r\ problem. Therefore, scalar fields fast-roll during inflation and one 
expects a considerable evolution of the gauge kinetic function, which is modulated by 
these fields. Indeed, in Ref. [107] it was shown that an expectation of /// ~ H is quite 
generic. The time dependence of the mass in these theories may be modulated by the 
same or additional DoF through the Higgs mechanism. In the same work it was shown 
that the required scaling of the gauge field mass, i.e. mota, can be achieved if the mass 
of the Higgs field is mu ~ H. This, again, is very reasonable due to corrections from the 
Kahler potential. However, to be a gauge field, the gauge coupling constant of the vector 
field must be small. As was discussed above, this means that only the kinetic function 
with / oc a~ 4 scaling is applicable in this case. Fortunately, this is a case which have the 
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richest phenomenology. 

In the context of implementing vector curvaton scenario in the particle physics theories 
it is important to note that models considered in this thesis involve only Abelian vector 
fields. However, most of gauge bosons in simple extensions of the Standard Model (SM) 
are non- Abelian. Therefore, an investigation of the particle production and the gener- 
ation of the curvature perturbation by non-Abelian vector fields is desirable (a related 
work can be found in Refs. [85, 86]). 

The investigation of the very early Universe is exciting for two reasons. First, it offers 
a possibility to understand the origin and history of the observable structure in the Uni- 
verse. Secondly, it serves as the giant laboratory to constrain theories of the fundamental 
physics. From the second point of view the research in cosmology is complimentary to 
the research in particle physics which can be tested by large experiments, such as LHC. 
At the time of writing LHC just started operating and everyone is looking forward with 
excitement for new discoveries. First of all, the detection of the Higgs boson is expected. 
This would prove the existence of fundamental scalar fields in Nature. If it is not discov- 
ered, the particle physics models without a fundamental Higgs field will become favorable, 
such as technicolor. But for inflationary model building until very recently only scalar 
fields were considered for the generation of the curvature perturbation. If such a field is 
not discovered, alternative models will become more attractive. However, currently the 
only alternatives being explored in the literature are vector fields. 

Another exciting possibility is the discovery of signatures of physics beyond SM. This 
will have a profound significance for particle physics as well as early Universe theories. If 
these signatures will be compatible with the supersymmetry, it will be a strong assurance 
that investigation of supersymmetric or supergravity models of inflation is the fruitful 
direction. 

From the astronomy side a large contribution towards the particle physics theories 
will be provided by the observations of recently launched Planck satellite. The most 
relevant questions for this thesis which Planck is expected to answer are: does the pri- 
mordial curvature perturbation have a detectable level of non-Gaussianity and statistical 
anisotropy? If non-Gaussianity and anisotropy is detected and if it is of the form sug- 
gested in this thesis, it will prove the non-negligible contribution of vector fields to the 
primordial curvature perturbation. This will provide a new observable allowing to probe 
the gauge field content of the effective field theory which governs the physics at energies 
when cosmological scales exit the horizon. 
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First note that in the equilateral configuration the unit vectors k a satisfy ki +k2 = — k3, 
where a = 1,2,3. If we define scalar products of the unit 3-vector, A, with each k a as 
W a = W • k a , then in the equilateral configuration W\ + W2 = — W3 and 

Wf + Wf + Wf = 2 {Wl + W1W2 + Wo 2 ) ; 

W1W2 + W2W3 + W3W1 = - (W? + W1W2 + Wl) ; (A.l) 
W%W% + VF 2 2 W 3 2 + WlWl = {Wl + W1W2 + w|) 2 . 

Let us define a vector Wj_ which is the projection of W to the plane containing vectors 
ki, k 2 and k 3 (see Figure A.l). Then the scalar product of these vectors and W is the 
same as the product with Wj_: 

W • k a = W ± • k a . (A.2) 

Without loss of generality we can assume that the angle between Wj_ and ki is ip: 

Wi = W • ki = W± ■ ki = W± cos tp, (A.3) 

where W± = |Wj_|. In equilateral configuration the angle between vectors ki and k 2 is 
27r/3, and W2 becomes 

W2 = Wi • k 2 = W± cos ^ip + = -W± cos ip + ^ sin cp^j . (A. 4) 

From the last two equations we get 

Wl + WiW 2 + Wl = -Wl. (A.5) 
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Figure A.I.: W± is the projection of the unit vector W into the plane of three vectors ki, 
k2 and k3. <p is the angle between W± and ki. In this figure the equilateral 
configuration of k a is shown. 

Putting this result back into Eq. (A.l) we find 

Wf + Wf + Wf = \W\\ 

W X W 2 + W 2 W 3 + W Z W X = -\Wl\ (A.6) 

w?w% + w%wl + wlwl = ^wf_ . 
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B. Scale Invariant Perturbation 

Spectrum of the Vector Field with 
Time Varying Kinetic Function 



In section 3.5.2 it was stated that the vector field with the time varying kinetic and mass 
terms in Eq. (3.169) acquires a scale invariant spectrum if the kinetic function scales as 
/ ot a°, where a = — 1 ± 3, and the mass scales as m oc a?, where (3 = 1. Here we will 
prove this result following Ref. [107], where it was derived by Dr. K. Dimopoulos. 

The equation of motion of the transverse modes is calculated in Eq. (3.178). For 
convenience let us rewrite it here using the conformal time 



w''+2-w' + 
a 



■-(a + 4)(a - 2) (aH) 2 + (aM) 2 + k 2 



w + = 0, (B.l) 



where primes denote derivatives with respect to the conformal time r. This equation is 
simpler that the one of the longitudinal mode. Thus we will find the value of a first, and 
then consider the equation for the longitudinal mode to determine f3. 

The value of a can be readily deduced by noting that Eq. (B.l) reduces to the equation 
of motion of a scalar field in Eq. (2.79) if 1 

a = -l±3. (B.2) 

Then in subsection 2.2.2.5 it was calculated that the scalar field acquires a scale invariant 
perturbation spectrum if the field is effectively massless and is initially in the Bunch- 
Davies vacuum state. Therefore, by analogy we conclude that transverse modes of the 
vector field acquire the scale invariant perturbation spectrum if it is effectively massless, 
has Bunch-Davies initial conditions and the kinetic function scales as / oc a 2 or / oc a~ 4 . 

However, the perturbation spectrum of the longitudinal mode depends not only on the 
scaling of the kinetic function but on the scaling of the mass as well, i.e. on the value of 



In Fourier space Eq. (2.79) becomes <j>'l + 2^(f>' k + [(am) 2 + k 2 ] </> k = 0. 
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(3. To determine this parameter let us rewrite Eq. (3.179) as 

1 



„ 4 - a + 2/3 , 

10 1 1 Wu + 



(a -2) (2 - a + 2/3) t~ 2 + kf 



wu = 0, 



(B.3) 



where we have also taken into account that the field has to be light for the perturba- 
tion spectrum of transverse modes to be scale invariant and we used the substitution 
r = — (aH)^ 1 valid in the de Sitter space-time. This equation can be solved using the 
vacuum initial conditions. For the longitudinal mode they are 



lim = 7 ( 

fer— >-oo \j2k 



-ikr 



(B.4) 



where 7 is the Lorentz boost factor defined in Eq. (3.184). For the light vector field it is 
equal to 

k/a 



M ' 



(B.5) 



where M is the mass of the physical vector field and is defined in Eq. (3.171). Solving 
Eq. (B.3) with initial conditions in Eq. (B.4) we find 



Wu 



k \J — T7r e 



aM 2a sin (iris) 



[J v (-fcr) - (-fcr)] 



(B.6) 



where J v denotes Bessel function of the first kind, and 



v = -^2(a - 2)(2 - a + 2/3) + (5 - a + 2/3) 2 . 



(B.7) 



At late times, when the mode exits the horizon, the dominant term of the above solution 
approaches 



lim w\ 



\J — T7T e - ' 



!(«+!) 



k 



V{\-v) a sin(vr^) \M J \2aH 
With this solution we find that the power spectrum is given by 



Vu 



k 3 
2tt2 



lim w\ 

fcr^-o 1 



167T 



(HV ( HY 



sin 2 (7n/)[r(l-i/)] 2 \M ) \2ttJ \2aH 



5-2u 



The above expression becomes scale invariant if v = 5/2, and V\\ becomes 

P "= 9 (S) 2 (S 2 - 



(B.« 



(B.9) 



(B.10) 
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Using Eqs. (B.7) and (B.2) we find that v = 5/2 is achieved if 

= -i(3±5). (B.ll) 

However, the value f3 = —4 must be disregarded. This can be seen using the definition 
of the mass M in Eq. (3.171) 

-^oca- 3 -*/ 2 . (B.12) 
k/a 

The above expression is a decreasing function of a with any value of a in Eq. (B.2). 
Thus, with j3 = —4 the vector field is massive at early times, which is contradictory to 
the requirement for scale invariance of the perturbation spectrum of transverse modes. 
Therefore, only the value /3 = 1 is allowed. 
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